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IN  THE  CRITICAL  CASE 


SINGULARLY  PERTURBED  EQUATIONS 

*  * 

A.  B.  Vasil 1 eva  and  V.  F.  Butuzov 

Introduction 

The  present  monograph,  essentially  the  sequel  to  the  book  "Asynptoti  : 
Expansions  of  Solutions  of  Singularly  Perturbed  Equations"  (13j,  is  devoted 
to  the  asymptotic  theory  of  differential  equations  with  a  small  parameter 
p  before  the  derivative,  such  as 

U  =  F(Z,y,t),  =  f(Z,v,t)  ,  U> 

and  to  other  related  problems  involving  asymptotic  behavior  for  small 
Such  problems  are  said  to  be  singularly  perturbed. 

The  difficulty  with  the  construction  of  the  asymptotic  expansion  of 
the  solution  of  system  (1)  arises  from  the  fact  that  for  ..  =  0  tin- 
order  of  the  system  decreases;  as  a  result,  the  solution  of  the  deqenerati 
system 

°  =  F(z,y,t),  ^  =  ft^^t)  ,  (2) 

cannot  in  general  satisfy  all  of  the  supplementary  conditions  prescribed 
for  (1).  On  account  of  this  singularity  the  asymptotic  expansion  of 
the  solution  of  system  (1)  cannot  be  constructed  solely  in  the  form  of 
an  "ordinary"  series  in  powers  of  (reqular  series)  but  a  boundary 

series  must  be  added  whose  terms  are  important  only  in  a  neighborhood  of 
those  points  at  which  the  given  supplementary’  conditions  for  (1)  are  not 
satisf ied  by  ( 2)  . 

* 

Originally  published  in  1978  by  Moscow  State  University.  Translated  from  the 
Russian  by  F.  A.  Howes  with  the  editorial  assistance  of  R.  E.  O'Malley,  Jr. 
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The  techniques  for  constructing  regular  and  boundary  series  are  des¬ 
cribed  in  detail  in  [13];  in  addition,  methods  for  estimating  the  remainder 
terms  are  also  discussed  there. 

All  of  the  problems  discussed  in  [13 ]  were  characterized  by  the 
fact  that  the  equation  F("z,y,t)  =  0  had  one  or  several  isolated  solutions 
~z  .  However,  in  applications  one  frequently  encounters  cases  where  this 
equation  has  a  family  of  solutions  which  depends  on  several  arbitrary 
functions.  We  shall  call  such  cases  critical  and  we  shall  consider  them 
■in  the  present  monograph.  Many  of  these  results  were,  obtained  in 
[3»k,9,10,ll,12,l5 ]  •  The  critical  case  can  be  distinguished  analytically 
by  definite  signs:  some  eigenvalues  of  a  special  matrix  are  identically 
zero.  We  also  discuss  several  problems  of  concrete  physical  importance 
in  a  number  of  fields:  problems  in  kinetics,  problems  in  the  theory 
of  semiconductors,  numerical  difference  schemes,  etc. 

The  techniques  for  constructing  the  asymptotic  expansion  in  this 
critical  case  are  basically  the  same  as  those  in  [13 ] -  For  a  thorough 
understanding  it  is  advisable  for  the  reader  to  become  acquainted  with 
the  first  three  chapters  of  that  book.  However,  in  the  present  monograph, 
for  an  understanding  of  at  least  the  formal  aspects  of  the  construction 
of  the  expansion,  it  is  not  necessary  to  consult  [13].  Concerning  the 
estimates  for  the  remainder  terms  of  the  asymptotic  expansions,  the 
reader  should  consult  [13]*  The  only  exception  is  Chapter  4  where  the 
proof  is  somewhat  different  and  so  it  is  given  in  detail. 

The  presentation  is  sufficiently  elementary  that  it  can  be  under¬ 
stood  completely  by  those  concerned  with  applied  questions. 
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Chapter  1 

i 

Weakly  Nonlinear  Singularly  Perturbed 
Equations  in  the  Critical  Case: 

Initial  Value  Problems 

§1  Singularly  Perturbed  Differential  Equations 
1.  Statement  of  the  Problem.  We  consider  the  differential  equation 

dx 

^dt  =  A(t)x  +  M.f(x,t,n)  ,  (1) 

where  u  >  0  is  a  small  parameter,  x  and  f  are  m- dimensional  vector 
functions,  A(t)  is  an  (m  x  m) -matrix  and  0  <  t  <  T  .  A  solution  of 

» 

equation  (1)  should  satisfy  the  initial  condition 

*  .  0 

x(0,*i)  *  x  (2) 

If  we  formally  set  u  =  0  in  (1)  then  we  obtain  the  reduced  equation 

A(t)x  =  0  .  (3) 

*7  I*  det  A(t)  /  0  for  0  <  t  <  T  ,  equation  (3)  has  the  unique 

solution  x  =  0  .  In  [13]  it  was  shown  that  if  the  eigenvalues  \^(t) 

*  of  A(t)  satisfy  for  0  <  t  <  T  the  inequalities 


I 


ReXi(t)  <0  (i  -  l,...,m)  , 

% 


then  the  solution  x(t,u)  of  the  problem  (1),  (  ?)  converges  as  n  0 
j  to  x  e  0  for  0  <  t  <  T  . 


Suppose,  however,  that  det  A(t)  =  0  for  0  <  t  <  T  .  Then 
equation  (3)  has  infinitely  many  solutions  and  the  question  arises: 
under  what  conditions  will  the  solution  x(t,|i)  of  the  problem  (1),  (2) 
converge  as  u  -*  0  to  one  of  these  solutions,  and  in  particular,  to 
which  one?  The  present  section  is  concerned  with  this  question  as  well 
as  with  the  question  of  the  construction  of  the  asymptotic  expansion  of 
x(t,u)  with  respect  to  u  . 

We  impose  several  additional  conditions  on  equation  (l).  All  of 
these  conditions  will  not  be  formulated  at  the  same  time  but  as  they  are 
stated  in  the  text.  They  will  be  denoted  by  the  numerals  I,  II,  ... 

The  first  condition  concerns  the  smoothness  of  A(t)  and  f(x,t,p.)  . 
We  require  sufficient  smoothness  in  order  to  construct  the  desired 
asymptotic  expansion.  A  more  precise  formulation  of  Condition  I  will  be 
given  in  Section  3  after  we  describe  the  algorithm’s  construction;  until 
then  we  formulate  this  condition  as  follows. 

I.  Suppose  that  A(t)  and  f(x,t,p.)  are  sufficiently  smooth  for 
0  <  t  <  T  and  for  (x,t,u)  in  the  domain  P(x,t,u)  =  D(x,t)  X  F°>M-q]  , 
where  D(x,t)  is  a  domain  in  (x,t) -space  and  u0  is  a  positive  constant. 

The  following  two  conditions  are  concerned  with  the  eigenvalues 
X^(t)  (i  *  l,...,m)  of  A(t)  .  Note  that  the  assumption  that  detA(t)^-0 
for  0  <  t  <  T  implies  that  at  least  one  of  the  X^(t)  is  identically 
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II.  Suppose  that  for  0  <  t  <  T  the  following  conditions  hold: 

X  (t)  =0  (i  =  l,...,k  ;  k  <  m)  ,  (b) 

ReX^t)  <0  (i  =  k+l,...,m)  .  (5) 

Remark.  In  [13 ]  the  initial  value  problem  was  studied  under  the 
assumption  that  condition  (5)  was  satisfied  for  all  i  =  l,...,m  (that  is, 
the  "noncritical"  case).  If  at  least  one  of  the  X^(t)  has  a  positive 
real  part  then,  generally  speaking,  the  solution  of  the  initial  value 
problem  is  unbounded  as  p  -  0  . 

III.  Suppose  that  there  are  k  linearly  independent  eigenvectors 
e^(t)  (i  =  1,  ...,k)  of  A(t)  corresponding  to  the  k  identically  zero 
eigenvalues  for  each  t  in  [0,T|  . 

Thus  we  are  considering  cases  where  the  number  of  linearly  independent 
eigenvectors  corresponding  to  X  s  0  is  equal  to  the  multiplicity  of 
X  =  0  .  For  the  remaining  eigenvalues  for  which  ReX^t)  <  0  ,  neither 
their  multiplicity  nor  the  number  of  eigenvectors  corresponding  to  them 
is  of  importance;  indeed,  both  of  these  quantities  can  change  as  t 
varies. 

2.  Algorithm  for  the  Construction  of  the  Asymptotic  Expansion  of  the 
Solution.  As  we  already  stated,  our  goal  is  the  construction  of  the 
asymptotic  expansion  of  the  solution  of  problem  (1),  (2).  In  order  to 


j 


oi 
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achieve  this  we  will  follow  the  same  procedure  as  that  adopted  previously 
in  the  book  [13].  First  we  develop  an  algorithm  (rule)  for  the  construc¬ 
tion  of  certain  formal  series  which  determine  the  structure  of  the  solu¬ 
tion.  In  a  neighborhood  of  the  leading  terms  of  these  series,  there 
exists  a  solution  of  the  problem,  and  the  series  is  itself  an  asymptotic 
expansion  of  this  solution.  The  proof  will  be  given  in  Section  3  . 

Thus  we  construct  a  series  formally  satisfying  equation  (1)  and 
condition  (2),  and  having  the  form 

x(t,p.)  =  x(t,u)  +  ^x(t,u)  ,  (6) 

where 

x(t,u)  =  xQ(t)  +  ux^t)  +  ...  +  M.nxn(t)  +  ...  (7) 

is  called  the  regular  series,  while 

rrx(T,ij.)  =  *0x(t)  +  httiX(t)  +  ...  +  ^nTTnx(T)  +  ...  (8) 

is  called  the  boundary  series  for 
T  =  t/u  . 

The  coefficients  in  the  series  (7),  (8)  are  determined  by  formally 
substituting  (6)  into  (1),  (2)  and  equating  terms  with  like  powers  of  u 
according  to  a  definite  rule  which  we  state  below.  First  note  that  the 
asymptotic  expansions  in  all  of  the  singularly  perturbed  initial  value 
problems  considered  in  [13  ]  are  constructed  in  the  form  of  series  having 
the  structure  (6).  This  structure  already  occurs  in  simple  examples. 
Consider,  for  example,  the  problem 
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>*£  *  “  * 1  •  x(0-‘i)  =  x°  • 

where  a  is  a  negative  constant.  The  exact  solution  is  x(t,u)  =  -t/a 
-  y./a2  +  (x°  +  p./a2)exp(  aT)  ,  and  it  consists  of  terms  of  the  type  (7) 
and  (8).  Both  terms  of  the  latter  type  converge  to  zero  exponentially 
as  T  -•  co  >  and  in  a  neighborhood  of  t  =  0  they  serve  as  a  correction 
to  the  regular  part  -t/a  -  \i/a.  which  does  not  satisfy  the  given  initial 
condition  x(0,|i)  =  x°  .  This  structure  of  (6)  reflects  such  behavior: 
TTx(t,vj.)  serves  as  a  correction  to  x(t,p,)  in  a  neighborhood  of  t  =  0  ; 
moreover,  it  converges  to  zero  exponentially  with  increasing  t  . 

The  coefficients  TT  x(-r)  of  the  series  (8)  will  be  called  boundary 
functions,  and  we  will  require  that  the  boundary  functions  converge  to 
zero  as  r  —  <*>  .  Thus  the  formal  algorithm  for  the  construction  of  the 
series  (7)  and  (8)  requires  that 

TT^t)  -  0  as  t  -  ®  .  (9) 

We  pass  now  to  the  procedure  for  determining  the  coefficients  in 

(7)>  (8).  To  this  end  we  first  represent  f(x,t,u)  in  the  form 

f(x(t,*i)  +  7x(t,4)  ,t,n)  =  f(x(t,u)  ,t,n) 

+  [f(x(T|i,y)  +  TTx(t,h)  ,Tp,,u)  -  f(x(Tp.,u)  JTU>U)  1  s  f  +  ^f  • 
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Here,  by  f  we  mean  the  expansion  of  f(x(t,n)  ,t,q.)  in  a  series  of  the 
type  (7),  while  by  TTf  we  mean  an  expansion  of  the  term  in  square 
brackets  in  a  series  of  the  type  (8)-,  namely, 

f  =  f^(t)  +  jif^(t)  +  ...  +  u  f^(t)  +  •••  , 

TTf  =  rr  f(T)  +  +  ...  +  Anf(T)  4  ••• 

We  perform  this  same  operation  on  A(t)x: 

A(t)  (x(t,)j,)  +  tTx(t,vj,))  =  A(t)x(t,|i)  +  A(tu.)TTx(t,u)  ^  Ax  +  TT(Ax)  . 

We  now  substitute  (6)  into  (1)  and  (2),  taking  account  of  the 
transformations  performed  on  f  and  Ax  : 

4  “V  •••>  *  £<V 4  “V 4  •••>  ‘ 

Ax  +  TT(Ax)  +  y.(f  +  h'f)  ,  (10) 

x(0,p.)  +  TTx(0,|jl)  =  x°  .  (11) 

Next  we  equate  coefficients  of  like  powers  of  u  on  both  sides  of  equations 

(10)  and  (U),  and  separating  those  terms  depending  on  t  and  those  depending 

on  t  ,  we  obtain  equations  and  initial,  conditions  for  determining  the 
coefficients  x^t)  and  TT^t)  of  the  series  (7)  and  (8)  . 

For  Xg(t)  we  obtain  a  linear  homogeneous  system  of  algebraic 


equations 
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A(t)x0(t)  =  0  ,  (12) 

which  coincides  with  the  reduced  equation  (3).  By  virtue  of  Condition  III 
the  general  solution  of  (12)  can  be  written  in  the  form 
_  k 

x  (t)  =  Z  a.  (t)e  (t)  ,  (13) 

u  i=l  1  1 

where  e^(t)  (i  =  l,...,k)  are  the  linearly  independent  eigenvectors 
corresponding  to  the  zero  eigenvalues  of  A(t)  ,  and  cn(t)  are  arbitrary 
scalar  functions. 

Remark.  By  virtue  of  Condition  III  the  rank  of  the  matrix  A(t)  is 
equal  to  m-k  for  each  t  in  [0,T],  that  is,  there  is  a  minor  of  order 
m-k  (in  general,  not  the  same  for  all  t  )  which  is  nonzero,  and 
consequently,  the  system  (12)  has  k  linearly  independent  solutions 
(eigenvectors)  e^t)  for  i  =  l,...,k  .  If  this  nonzero  minor  can  be 
found,  it  is  easy  to  construct  eigenvectors  e^t)  having  the  same  degree 
of  smoothness  as  the  matrix  A(t)  .  If  there  is  no  such  minor  of  A(t) 
then  the  question  of  the  degree  of  smoothness  of  the  e ..  (t)  becomes 
more  involved.  From  [27]  it  follows  that  it  is  possible  to  construct 
eigenvectors  e^(t)  (i  =  l,...,k)  having  the  same  degree  of  smoothness 
as  the  matrix  A(t)  .  Such  eigenvectors  are  used  in  (13)  and  below. 

If  we  introduce  the  (mxk)  -matrix  e(t)  whose  columns  are  these 
e^(t)  and  the  k-dimensional  vector  function  a(t)  whose  components 
are  a^Ct)  ,  then  it  is  possible  to  write  (13)  in  the  form 


-  +*9*< 
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xQ(t)  =  e(t)a(t)  .  (11) 

For  TT^x(t)  we  obtain  a  linear  homogeneous  constant  coefficient 
system  of  differential  equations,  namely 

sV '  A(0,V  •  <«> 

The  general  solution  of  this  system  can  be  written  in  the  form  (cf.  for 
example  [ 26 ] ) 

k  m 

TTqx(t)  =  £  c.e  (0)  +  I  CjV  (T)exp(X  (o)t)  ,  (16) 

i=l  i=k4 1 

where  c^i  =  l,«»->m)  are  arbitrary  constants,  e^O)  (i  =  l,...,k)  are 
the  eigenvectors  of  A(0)  corresponding  to  the  zero  eigenvalues,  and 
w^(t)  (i  =  k+l,...,m)  are  known  vector  functions  whose  components  are 
polynomials  in  t  .  [If  h^,...,hn  is  a  Jordan  chain  of  vectors  cor¬ 
responding  to  an  eigenvalue  X  of  A  such  that 

Ah,  =  Xh,  Ah  =  Xh  +  h, Ah  =  Xh  4  h  ,  , 

1  1*  2  2  1*  *  n  n  n-1  * 

then  there  are  n  linearly  independent  solutions  xr(r)  =  wr(T)exp(Xt) 

(r  =  l,...,n)  of  the  system  ^  =  Ax  ,  where 
Tr-1  Tr-2 

VT)  '  TIoJT  hi 4  TT^TT  V  —  4  Vr  *  1 . ">  •  l 

By  virtue  of  condition  (5)  the  second  term  in  the  right-hand  side  of 
(16)  converges  to  zero  as  T  -*  •  .  Therefore,  in  order  that  condition  (9) 
hold,  it  is  necessary  to  set  c^  *  0  (i  =  l,...,k)  . 

The  initial  condition  for  ^(t)  is  obtained  by  equating  the 
coefficients  of  the  zeroth  power  of  4  in  (11)  ,  namely 
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TT0x(0)  =  x°  -  xQ(0)  =  x°  -  t^a.(0)ei(0)  . 

Substituting  now  into  (16)  and  noting  that  =  0  (i  =  we 

obtain 

1c  in 

lQ.(0)e.(0)+  E  c.w.(O)  =  x°  .  (18) 

i=l  1  1  i=k+l  1  1 

The  system  (l8)  is  a  linear  algebraic  system  of  m  equations  in  the  m 
unknowns  ct.(O)  (i  =  l,...,k)  and  c^(i  =  k+  1,  •  •  •  ,m)  .  By  virtue  of  the 
linear  independence  of  the  column  vectors  e^(O)  (i  =  l,...,k)  and 
wi(0)  (i  =  k+l,...,m)  the  system  (18)  has  a  unique  solution. 

Thus  TTqX(t)  is  completely  determined.  By  virtue  of  (|3)  it  is 
clear  that  there  exist  constants  c  >  0  and  k  >  0  such  that  !I^qX(t)H  < 
cexp(-KT)  for  T  >  0  .  (The  symbol  !lxJ!  denotes  the  norm  of  a  vector 
(matrix)  x  which  is  defined,  for  example,  as  the  sum  of  the  absolute 
values  of  its  components  (elements).  ]  The  function  xQ(t)  is  not  defined 
until  the  functions  on^t)  (i  =  l,...,k)  whose  initial  values  0^(0) 
are  found  from  (18)  are  first  defined.  Let  us  set  a(0)  =  a°  . 

For  xx(t)  we  obtain  the  linear  nonhomogeneous  system  of  algebraic 
equations 

_  _  dxn(t) 

A(t)xx(t)  =  -f(xQ(t),t,0)  +  — -  =  cp(t)  .  (19) 

Since  det  A(t)  =  0  for  0  <  t  <  T  a  necessary  and  sufficient  condition 
for  the  solvability  of  system  (19)  is  that  its  right-hand  side  be  ortho¬ 
gonal  to  each  of  the  eigenvectors  gj(t)  (j  =  l,...,k)  of  the  adjoint 

matrix  A  (t)  corresponding  to  the  zero  eigenvalues.  [From  linear  algebra 

* 

it  is  known  that  the  matrix  A  (t),  which  in  the  present  case  is  simply 
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the  transpose  of  A(t)  ,  has  a  zero  eigenvalue  to  which  correspond  k 

linearly  independent  eigenvectors.  We  note  that  the  degree 

of  smoothness  of  the  vectors  g.(t)  is  the  same  as  that  of  the  matrix 

J 

A(t)  .  ] 

We  denote  by  <a,b)  the  scalar  product  of  two  m- dimensional 
vectors  a  and  b  ,  that  is,  the  sum  of  products  of  corresponding  compon¬ 
ents.  Thus  (taking  note  of  (lU))the  solvability  condition  for  (19) 
be  written  as 

<g^(t),  -f(e(t)a(t) ,t,0)  +  -^(e(t)a(t) ))  =  0  ,  j  =  l,..»,k  • 

We  will  obtain  a  system  of  k  nonlinear  differential  equations  for 

the  k  unknown  functions  cn(t)  .  If  we  call  g(t)  the  (k x m) -matrix 

whose  rows  are  the  vectors  g.(t)  (j  =  l,...,k),  then  it  is  possible  to 

J 

obtain  the  corresponding  system  in  matrix  form,  namely 

(g(t)e(t))^=  g(t)  (f (e(t)a(t)  ,t,0)  -  e'  (t)a(t) ) .  (20) 

The  initial  values  a(0)  =  a°  ,  as  noted  above,  are  found  from  (l8)  . 

From  the  fact  that  the  number  of  eigenvectors  e^t)  and  g^(t1  is 

equal  to  the  multiplicity  of  the  zero  eigenvalue,  it  follows  that  the 

determinant  of  the  (kxk) -matrix  (g(t)e(t))  is  nonzero.  For  if  this 

determinant  were  zero  for  some  t  ,  then  a  certain  nontrivial  linear 

k 

combination  of  its  columns  gives  a  zero  column,  that  is,  Z  =  0 

k  i*l  1  J  1 

(j  *  l,...,k)  .  It  follows  that  <g  ,  t  Y.e.)  =  O  ,  that  is,  the 
k  J  i-1 

eigenvector  e  =  C  Y.e.  of  A(t)  ,  corresponding  to  X  *  0  ,  is 
i*l  1  1 
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orthogonal  to  each  g.(j  =  l,...,k)  so  the  solvability  condition  is 

J 

satisfied  for  the  linear  system  Ax  =  e  .  But  this  implies  that  the 
zero  eigenvalue  of  A(t)  has  adjoint  vectors,  which  contradicts  the 
fact  that  the  number  of  eigenvectors  e^(t)  is  equal  to  the  multiplicity 
of  X  =  0  .  Thus,  the  aforementioned  determinant  is  nonzero,  and  there¬ 
fore  the  system  (20)  can  be  solved  for  ,  that  is. 


§  -  ro(a’t'  • 


(21) 


with  the  form  of  FQ  being  clear  from  comparing  (20)  and  (21)  . 


IV.  Suppose  that  the  equation  (21)  with  the  initial  condition 
a(o)  =  has  a  solution  a  =  a(t)  for  0  <  t  <  T  . 

Now  that  a(t)  is  determined,  the  solution  (Ih)  of  the  reduced 
system  (12)  is  complete. 

Let  us  introduce  in  the  space  of  variables  (x,t)  a  curve  L  con¬ 
sisting  of  two  components  (this  curve  is  the  graph  of  the  zeroth  pproxi- 
mation) : 


L1  =  { (x,t) :  x  =  xQ(  0)  +  V(T)  (T  >  0)i  t  =  0}  , 
Lg  "  {(x,t):  x  =  xQ(t);  0  <  t  <  T)  . 


V.  Suppose  that  the  curve  L  lies  in  the  domain  T)(x,t)  appearing 
in  Condition  I. 

Thus  the  zeroth  order  terms  in  the  series  (7)  and  (8)  have  been 
determined. 


The  general  solution  of  the  system  (19)  can  be  written  as 


k 

xn(t)  =  £  0.(t1e.(t)  4  x  (t)  =  e(t)0(t)  4  x  (t)  ,  (22) 

1  i=l  1  1 

where  x^(t)  is  a  particular  solution  of  (19)  and  0(t)  is  an  arbitrary 
k- dimens ional  vector  function. 

For  TT^x(t)  we  obtain  a  linear  nonhomogene ous  system  ~>f  differential 
equations 

dJT  x  _ 

=  A(0)TTlX4  TA'(0)tT0x(t)  4  f(x0(0)  4  TTox(t),0,0) 


-  f(xQ(0),0,0)  . 


(23) 


The  initial  condition  for  IT^t)  is  found  by  equating  the  coefficients 
of  the  first  power  of  p,  in  (11),  that  is. 


1^(0)  *  -x^O)  . 

The  general  solution  of  (23)  can  be  written  as 
km 

rr  x(t)  =  I  d  e  (0)  4  t  d  w  (T)exp(X  .(0)t)  4  TT  x(t)  ,  (24) 

x  i=l  1  1  i=k+l  11  1  x 

where  the  d^  are  arbitrary  constants,  w^(t)  are  the  same  vector  functions 

as  those  in  (16),  and  TT^x(t)  is  a  particular  solution  of  (23)  which,  it  is 

mm 

not  difficult  to  see,  can  be  chosen  so  that  ||TTx1(t)||  satisfies  the  same 
inequality  as  ||TTx0(t)||  ,  that  is,  ||TT^x(t)||  <  cexp(-nT)  for  T  >  0  . 

Remark.  The  positive  constants  c  and  k  in  the  estimate  for 

m * 

TT^Ct)  are,  in  general.,  not  the  same  as  those  in  the  estimate  for  tt^t)  . 
However,  with  a  view  to  simplifying  the  notation,  the  same  notation  will 
be  used  for  analogous  constants.  Constants  of  the  form  c  will  always 
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denote  an  upper  bound  on  such  values,  while  constants  of  the  forir.  >•. 
will  always  denote  a  lower  bound. 

As  in  the  case  of  TTqX(t)  we  require  that  =  0  ( i  =  l,...,k) 
and  then  from  the  initial  condition  for  fT  x(t) ,  we  obtain  a  linear 
algebraic  system  of  m  equations  in  the  m  unknowns  3.(0)  (i  =  l,...,k' 
and  d^(i  =  k+l,...,m)  : 

km  _ 

t  p  e  (0)  +  t  d.v  (0)  =  -5L(o)  -  rr  x(0)  ,  (&) 

i=l  1  1  i=k+l  1  1  L  1  ' 

which,  like  system  (18),  has  a  unique  solution. 

Thus  ThjX(T)  is  completely  determined,  and  it  is  obvious  that 

||TTjX(t)||  <  cexp(-KT)  for  T  >  0  .  The  expression  for  x  (t)  is  not 

determined  until  we  determine  p(t)  ,  whose  initial  value  p(0)  is 

obtained  from  (25)  .  It  can  be  determined  from  a  solvability  condition 

for  the  linear  system  of  algebraic  equations  relative  to  x„( t)  in  a 

manner  analogous  to  that  for  a(t)  .  We  obtain  the  linear  differential 

equation 

H  *  B(t)p  4  Fx(t)  ,  (26) 

where  B(t)  ,t)  is  a  known  matrix  and  F  (t)  is  a  known 

function.  By  linearity  the  system  (26)  with  the  initial  condition  p(0) 

has  a  unique  solution. 

The  determination  of  the  remaining  terms  in  the  series  (7)  and  (8) 
proceeds  analogously.  At  the  i-th  stage  an  arbitrary  vector  function 
(call  it  Y(t))  enters  the  expression  for  x^( tl  .  First  we  determine 


16 


\(0)  from  an  equation  like  (25)  and  then  from  the  solvability  condition 

for  x.  ,(t)  we  obtain  for  \(t)  a  linear  differential  equation  of  the 
l+l 

form  of  (26):  =  B(t)v  +  F. (t)  ,  from  which  y(t)  is  finally 

ClU  1 

determined. 

The  boundary  functions  Tf  x(«r)  are  constructed  like  TT^x(t)  and 
also  satisfy  the  exponential  estimate 

||TT^x(t)|!  <  cexp(-*f)  for  t  >  0  .  (27) 


3.  An  Estimate  of  the  Remainder  Term.  Let  us  denote  by  Xn(t,p.)  the 
n-th  partial  sum  of  the  series  (6),  that  is, 

n  i  - 

X_(t,u)  =  £n  (x  (t)  +  tt  x(T))  . 

We  note  that  for  the  determination  of  the  terms  appearing  in  Xn(t,p.) 
it  is  sufficient  that  A(t)  and  f(x,t,n)  have  continuous  partial 
derivatives  with  respect  to  all  arguments  up  to  order  n  inclusive.  How¬ 
ever,  for  the  proof  of  Theorem  1.1  below  it  will  be  necessary  that  A(t) 
and  f(x,t,u)  have  a  higher  degree  of  smoothness.  Tt  is  now  possible 
to  determine  more  precisely  the  extent  of  the  domain  D(x,t)  appearing 
in  Condition  I  .  Let  us  introduce  a  5- tube  of  the  curve  L  ,  that  is,  the 
set  of  points  (x,t)  whose  distance  from  the  curve  L  does  not  exceed  &  . 

For  the  domain  D(x,t)  we  cam  take  an  arbitrary  6-tube  of  the 
curve  L  where  6  is  independent  of  4  .  Condition  I  can  now  be 


formulated  as: 
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I.  Suppose  that  A(t)  in  0  <  t  <  T  and  f(x,t,u.)  in  the  domain 
D(x,t,u)  =  D(x,t)  X  [0,n  )  have  continuous  partial  derivatives  of  order 
(n  +  2)  inclusive  (with  respect  to  each  argument) . 

Theorem  1.1.  Under  Conditions  I  -  V  there  exist  positive 
constants  and  c  such  that  for  0  <  u  <  Uq  the  solution  x(t,|i) 

of  the  problem  (1),  (2)  exists  in  the  interval  [0,T],  is  unique  and 
satisfies  the  inequality 

||x(t,n)  -  Xn(t,n)||  <  cuml(o  <  t  <  T)  .  (28) 

(In  place  of  an  inequality  of  the  type  (28)  we  will  also  use  the  notation 
x(t,u)  -  Xn(t,u)  =  &(M.m1).) 


Proof.  Let  us  set  f(t,u)  =  x(t,u)  -  Xml(t,u)  .  Substituting 
X  «=  xn+3_(t,M')  +  £  into  (l)  and  (2)  we  obtain  for  £;(t,p.)  the  initial 


value  problem 

=  A(t)|  4nfx(t,u)F  4  G(?,t,n)  , 

(29) 

e(o,i*)  =  o  , 

(30) 

where  fx(t,p)  =  fx(xQ(t)  +  ^qx(t) ,t,u)  and 

G(?,t,n)  =  A(t)Xnfl(t,n)  +  Hf(Xml(t,u)  4  £,t,u) 
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The  function  G(?,t,p)  has  the  following  two  important  properties 
which  are  established  in  the  same  way  as  in  Subsection  h  of  Section  10 

of  [13]  : 

1.  G(o,t,n)  =  &(nm2)  ; 

2.  If  l|?1(t,M.)!l  <  c-^p  and  II  ?2(t ,u)i!  <  c^p  for  0  <  t  <  T  and 

0  <  4  <  (c^  and  p^  axe  certain  constants),  then  there  exist  constants 

cQ  and  pQ  <  p^  such  that  for  0  <  t  <  T  and  0  <  p  <  pQ 

||G(5,,t,p)  -  G(SL,t,p)||  <  c_u2  max  ||?  -  . 

1  2  ^  [0,T }  1 

When  Property  2  holds  we  will  say  that  G(f,t,p)  is  a  contraction 

2 

operator  with  contraction  coefficient  of  order  &(p  )  for  P  =  6(p)  . 

Remark.  The  constant  pQ  appearing  in  Property  2  is,  generally 
speaking,  not  the  same  as  that  appearing  in  the  statement  of  Theorem  1.1. 
However,  for  the  sake  of  simplicity  of  notation  we  will  use  the  same 
symbol  j,Q  in  all  bounds  involving  conditions  of  sufficient  smallness 
for  p  .  It  is  clear  a' priori  that  among  all  such  constants  p^  the 
least  one  will  furnish  a  positive  bound. 

For  the  problem  (29),  (30)  we  now  transform  the  unknown  function 
*>y  ?(t,p)  **  T(t)w(t,p)  ,  where  T(t)  is  a  certain  differentiable 

(m  x  m)  -matrix.  Thus  we  obtain  the  initial  value  problem 

*  T-1(t)A(t)T(t)w  +  pC(t,p)w  +  GjJwjtjp)  ,  (31) 


v(0,p)  =  0 


9 


(32) 


where 


c(t,n)  =  T-1(t)[fx(t,n)T(t)  -  T* (t) ]  and 
G1(w,t,u)  =  T-1(t)G(T(t)wJt,n)  . 


We  choose  for  T(t)  a  matrix  which  transforms  A(t)  in  the  interval 
[0,T]  into  the  block-diagonal  form,  that  is, 


T”  (t)A(t)T(t)  = 


’A]_(t)  0 


0  A2(t). 


where  A^(t)  is  an  (m-k)  x  (m-k) -matrix  whose  eigenvalues  satisfy 
ReX^(t)<0(i  =  k+l,...,m)  and  A,,(t)  is  a  (k  x  k)-matrix  with  k 
zero  eigenvalues  X^(t)  =  0  (i  =  1, . . . ,k)  .  Such  a  nonsingular  matrix 
T(t)  exists  and  it  is  as  differentiable  as  A(t)  (cf.  [P7]).  Indeed, 
in  our  case,  A2(t)  =  0  for  0  <  t  <  T  . 


We  set 


G2(u,v,t.u,)' 


G3(u,v,t,u), 


(where  u  and  G5  have  (m-k)  -  (and  v  and  G  k-)  components) ,  and 


we  divide  the  matrix  C(t,p)  into  compatible  blocks 

C(t,p) 


cn(t»4)  c12(t,^> 


Cg^(t,|A)  C00(tf\l) 


22' 


where  the  C.  .  are  continuous  and  bounded  for  0<t<T,0<p,<u  • 

From  (31) ,  (32)  we  obtain 
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=  Ax(t)u  +  uC;LI(t>p.)u  +  uC12(t,u)v  +  G2(u,v,t,M.)  , 

It  =  +  C22(t,u)v  +  ~G^(u,v,t,u)  , 

u(0,u)  =  v( Ojp.)  =  0  . 


(33) 

(3M 


The  functions  G2(u,v,t,u)  and  G^(u,v,t,u)  have  the  same  two 
properties  as  G(?,t,u)  . 

Let  us  denote  the  fundamental  matrices  of  the  homogeneous  linear 

systems 

^  d!  =  Al(t)u  ’  S  =  C22(t^)v  * 

by  U(t,s,u)  and  V(t,s,u)  ,  respectively  (U(s,s,u)  =  L’  .  ,  the 

m-  k 

(m-k)  X  (m-k) -identity  matrix,  and  V(s,s,u)  =  E^)  .  Clearly  the 
matrix  V(t,s,u)  is  bounded,  and  since  the  eigenvalues  of  A^( t)  have 
negative  real  parts,  it  follows  that  U(t,s,u)  satisfies  the  inequality 

||U(t,s,p.)||  <  c  exp[-*^~^-]  for  0<s<t<T,0<u<u0  . 

The  proof  of  this  estimate  is  given  in  Lemma  3.2  of  [13]  .  We  now 
replace  the  problem  (33) >  (3*0  with  the  equivalent  system  of  integral 
equations 

t 

u(t,p)  =  J*  U(t,s,u)  [C11(s,|i)u(s,u)  +  C.?(s,u)v(s,p.) 

0  ■L<; 

+  ^  G2(u»v»s,h)  Ids  , 
t 

v(t,u)  =  J  v(t,s,u)  [c01(s,n)u(.s,u)  +  -  Cr  (u,v,s,u)  Ids  , 

0  *■  U  0 


(35) 
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and  then  by  applying  the  method  of  successive  approximations  as  in 
[13;  Sec.  10],  it  is  possible  to  prove  that  the  solution  u(t,u) ,v(t,u) 
of  (35)  exists,  is  unique,  and  satisfies  the  estimate  u(t,u)  =  &(p.ml)  , 
v(t,u)  =  ©(uml)  . 

Thus  it  follows  directly  that 

|(t,u)  =  x(t,u)  -  Xwl(t,u)  =  &(uml).  (36) 

Since  Xml(t,p)  -  Xn(t,^)  =  |irWl(xrHl  +  TT^x)  =  ^(u^1)  ,  from  (36) 
we  obtain 

x(t,u)  -  xn(t,u)  =  ^(u1*1) 

and  Theorem  1.1  is  proved. 

Remarks .  1.  From  (28)  and  the  exponential  bound  (27)  for 

^*(*0  (i  =  0,1, ...,n)  it  follows  that  lim  x(t,u)  =  x  (t)  for  0  <  t  <  T 

U'-’O 

that  is,  the  solution  x(t,(j,)  of  the  problem  (1),  (2)  converges 
as  n  -*  0  to  one  of  the  solutions  of  the  degenerate  system.  For 
0<tg<t<T  (t^  fixed  as  p.  —  O)x^(t)  is  the  leading  term  of  the 
asymptotic  expansion.  In  certain  problems  (see,  for  example.  Subsection  3 
of  Section  3  below)  one  is  interested  in  a  precise  representation  of 
the  leading  term.  Then  to  determine  the  initial  values  <^(0)  of  the 
functions  Of^(t)  appearing  in  x^(t)  (see  (13))  it  is  convenient  to 
have  a  system  of  equations  involving  the  constants  c^  of  the  leading 
boundary  function  terms.  Such  a  system  can  be  obtained  by  taking  the 
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scalar  product  of  (l8)  vith  g.(0)  (j  =  after  which  the  tenr.s 

«J 

containing  vanish,  that  is  , 

k  a 

t,  <g,(0) ,e.(0))a  (0)  =  (g .( 0) ,x  >  (j  =  l,...,k)  .  (37) 

i=l  0  1  J 


2.  The  construction  of  the  asymptotic  expansion  and  the  bound 
on  the  remainder  term  have  been  obtained  under  the  condition  that  there 
exist  k  linearly  independent  eigenvectors  corresponding  to  the  zero 
eigenvalues  of  A(t)  .  In  the  case  when  the  number  of  linearly  independent 
vectors  corresponding  to  the  zero  eigenvalue  is  less  than  its  multiplicity 
the  asymptotic  expansion  will  contain  fractional  powers  of  u  .  We  will 
not  consider  the  general  problem  but  we  will  illustrate  the  occurrence 
of  fractional  powers  of  ^  in  a  model  system  of  two  equations, where 

(a11(t)  ai2(t)\ 

)  • 

*2l(t)  *22(t)S 


Suppose  that  det  A(t)  =  a^a^  -  a^a^  =  0  >  a]j_  +  a22  =  0  40(1 

a11  ^  0  for  0  <  t  <  T  .  Thus  X(t)  =  0  is  an  eigenvalue  of  A(t)  of 
multiplicity  two,  to  which  there  corresponds  only  one  eigenvector  (since 
a11  /  0)  .  In  this  case  (1)  has  the  form 


d*i 

^"dT  *  Vl  +  a12X2  +  Hfi(xp*2»t,n) 

<iXp 

“■JET  ■  a2ixi 4  a22*a  4 


(38) 
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Multiplying  the  first  equation  by  -agl  and  the  second  oy  and 

adding  we  obtain  the  following  system  of  equations  for  =  = 
y  =  -a2lx1  +  anx2  : 


dz 
^  dt 


& 

dt 


+  nf  (z,^^,t,u)  , 

_  J- 


u 


iu 


=  -a 


'21 


^(zj 


y+  a?,  z 

-  ,t,u)  +  a^_]_^2  a. 


^21* 


"11 


>t,u) 


11 


all  z  -  a. 


y^a21z 


21  11  a. 


'll 


that  is,  a  system  of  the  form 

u4r  =  a(t)y  +  yF(z,y,t,n)  >  f£  =  C(z,y,t,u) 


(39) 


dt 


The  behavior  of  solutions  of  system  (39)  depends  critically  on  the  sign 
of  a(t)  .  Let  us  transform  (39)  fey  setting  z  =  ,  Y  -  v/d 

obtain  the  system 


dz 


Vt-nr  =  a(t)z2  +  JZ  F(zlf  ’ 


dt 

dz 


(40) 


The 
the  form 
l-X 


^~dt  =  G(zl>^  z2,ty^  ' 

corresponding  characteristic  equation  (cf.  (3-21)  in  [131)1*8 


a(t) 


G  (<p("t)  ,0,^*°) 

Z1 


=  X2  -  a(t)G  («o(t),0,t,0)  =  0  , 

21 
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where  z^  =  cp(t)  is  a  root  of  the  equation  G(z^,0,t,0)  =  0  . 

If  aG  <0  then  X  will  be  purely  imaginary  and  in  the  system 
Z1 

one  usually  sees  oscillations.  Solutions  of  the  initial  value  problem 

do  not  have  a  limit  as  p  -•  0  ,  but  are  bounded  and  oscillate  with  a 

frequency  of  order  —  (cf.  [19]).  If  however  aG  >0  the  roots  X 

Zl 

will  have  different  signs  (the  so-called  conditionally  stable  case)  and 
consequently  solutions  will  not  generally  be  bounded  as  p  -  0  .  None¬ 
theless,  if  we  impose  appropriate  boundary  conditions  on  (38),  it  is 
possible  to  carry  out  the  construction  of  the  asymptotic  expansion  of 
the  solution  as  in  [I3,  Sec.  14].  Such  an  expansion,  as  can  be  seen 
from  the  way  (40)  is  written,  will  contain  powers  of  .  This  case  is 
considered  in  more  detail  in  [24]  . 

3.  Let  us  now  consider  the  equation 

U2  -^  =  A(t)x  +  pf(x,t,p)  . 

2 

It  differs  from  (1)  in  that  the  term  multiplying  the  derivative  is  p 
and  not  p  .  Thus  the  coefficient  multiplying  the  derivative  is  of  a 
greater  degree  of  smallness  than  the  coefficient  multiplying  the  non¬ 
linear  term  f(x,t,p)  .  As  a  result,  the  asymptotic  expansion  of  the 
solution  of  the  initial  value  problem  contains  [along  with  the  regular 
part  x(t,p)  =  xQ(t)  +  px^(t)  +  ...  and  the  boundary  function 
TTx(r,p)  =  TTqx(t)  +  pTT^x(t)  +  ...  ,  depending  on  T  =  t/p]  the  boundary 
function  Px(s,p)  =  PQx(s)  +  pP^x(s)  +  ...,  depending  on  s  «=  t/p  . 

The  conditions  ,.1  the  matrix  A(t)  are  the  same  as  those  in  Subsection  1. 


Consequently, 
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xQ(t)  =  e(t)a(t)  , 

where  the  function  a(t)  is  defined  as  in  Subsection  2  by  a  solvability 
condition  in  the  equation  for  x^(t)  \  this  time  however  we  do  not  obtain 
a  differential  equation  for  a(t)  but  rather  an  algebraic  equation 

g(t)f (e(t)a(t) ,t,C)  =  0  , 

where  g(t)  is  the  same  matrix  as  in  Subsection  2. 

TTqx(t)  is  obtained  not  from  a  differential  equation  but  from  the 
algebraic  equation  A(0)TTqX(t)  =  0  .  It  follows  that 

TT0x(t)  =  e(0)h(r)  , 

where  h(x)  satisfies  the  differential  equation 

(gCOeCOjJ^IfUQCO)  +  e(0)h(T),0,0)  -  f (xQ(0) ,0,0) ] 

obtained  from  a  solvability  condition  for  h^x(T)  . 

For  PqX(s)  we  obtain  the  equation 
dPQ(s) 

-i-.A(0)P0x 
with  the  initial  condition 

PQx(0)  *  x°  -  xQ(0)  -  e(0)h(0) 

(h(0)  is  as  yet  unknown).  Since  this  equation  is  like  (15)  its  general 
solution  can  be  written  in  the  form  of  (16)  : 


ft 

i 


2  6 


k  m 

Prtx(s)  =  T  c.e.(O)  +  t  c  .w .  ( s  >  exn(X  .  (0)  s) 

0  i=l  1  1  i=k+l  11 

From  the  condition  that  P^x(s)  -  0  as  s  —  ®  we  have  =  0( i  =  l,...,k' 
and  so  the  initial  conditions  imply 

m  o  _ 

e(0)h(0)  +  £  c.w.(O)  =  x  -  x_(G) 

.  ,  n  1  1  o 

i=k+l 

for  the  determination  of  h(0)  and  c  (i  =  k+l,...,m)  .  Likewise  the 
function  P^x(s)  will  be  determined  completely,  and  by  finding  initial 
conditions  for  h(f)  we  can  finally  determine  this  function  h(f)  from 
its  differential  equation. 

The  essential  role  in  the  construction  is  played  by  the  (k  x  k) -matrix 

(g(t)e(t))-1g(t)fx(x0(t) ,t,0)e(t)  . 

We  require  that  its  eigenvalues  \n(t)  satisfy 

Re  v^(t)  < 0  (i  =  l,...,k  v  0  <  t  <  T  )  . 

If  this  condition  is  fulfilled  together  with  certain  others  it  is 
possible  to  prove  the  validity  of  this  asymptotic  expansion  with  the 
boundary  functions  having  the  exponential  bounds 

I|TTjx(t)||  <  cexp(-HT)  (r  >  0)  •, 
ll^xCs)!!  <  c  exp(-Ks)  (s  >  0)  . 

A  more  detailed  consideration  of  this  problem  has  been  given  in  [5  ]  . 


§2  Difference  Equations  with  Small  Stepsize 

In  this  section  we  will  consider  the  difference  equation 

x(tw-M')  =  B(t)x(t)  +  uf(x(t),t,u)  ,  ^1) 

in  which  x  is  an  m-vector  and  the  argument  t  varies  discretely  with 
small  stepsize  u  ,  that  is,  t  =  0  ,  u,  2u,  ...  (t  <  T)  .  Such  variatirr.s 
of  the  argument  occur,  for  example,  in  difference  schemes  for  many  interro- 
differential  equations  (cf.  For  simplicity  of  notation  we  will 

write  x(t)  in  place  of  x(t,u)  .  We  prescribe  the  initial  condition 

x(0)  =  x°  .  (^2) 

I.  Suppose  that  the  (m  x  m) -matrix  B(t)  has  for  0  <  t  <  T  the 
eigenvalue  X(t)  =  0  of  multiplicity  k  to  which  there  correspond  for 
each  t  k  linearly  independent  eigenvectors  e^t)  (i  =  l,...,k)  , 
while  the  other  eigenvalues  satisfy  the  condition  |X^(t)|  <  1  . 

In  [I3]  it  was  shown  that  the  asymptotic  properties  of  the  difference 
equation  (bl)  in  the  noncritical  case  (k  =  0,  that  is,  all  IX^t)!  <1) 
are  similar  to  the  asymptotic  properties  of  the  differential  equation 
(1)  in  the  noncritical  case  (all  ReX^(t)  <  0)  .  We  will  consider  the 
critical  case  (k  ^  0)  ;  the  asymptotic  expansion  turns  out  to  be 
similar  to  the  one  for  the  differential  equation. 

If  we  define  the  matrix  A(t)  =  B(t)  -  K  ,  then  it  clearly 
satisfies  Conditions  II  and  III  of  §1  . 


II.  Suppose  that  B(t)  and  f(x,t,p)  satisfy  the  same  smoothness 


conditions  as  A(t)  and  f(x,t,p)  in  §1  (cf.  Condition  I)  . 

For  the  problem  (Ul),  (U2)  it  is  possible  to  construct  an  asymptotic 
expansion  in  the  parameter  p  and  to  give  an  estimate  for  the  remainder 
term  as  in  §1.  We  write  the  solution  in  the  form  (cf.  (6)  -  (8)) 


x(t)  =  xQ(t)  +  px^t)  +  ...  +  ^(t)  +  uTT^t)  +  ...  .  (1*3) 

Substituting  (1*3)  into  (Ul)  and  (U2)  we  obtain 


x0(t*u)  +  px-^t+p)  +  •••  +  t1'qX(t+1)  +  uTr1x(T+l)  +  ...  = 
Bx  +  TT(Bx)  +  p(f  +  TTf)  , 

xQ(0)  +  px.^0)  +  ...  +  TTqx(0)  +  pTTjx(O)  +  ...  =  x°  . 


(UU) 

0*5) 


The  right-hand  side  of  (UU)  reduces  to  the  same  form  as  in  (10)  .  Equating 
coefficients  of  like  powers  of  p  in  (UU)  and  (U5)  as  in  (10)  and  (11)  we 
obtain  for  xQ(t)  a  system  of  equations 

xQ(t)  =  B(t)xQ(t)  or  A(t)xQ(t)  =  0 

It  follows  that  the  representation  (13)  (or  (14))  is  correct  for  xQ( t)  , 
that  is, 

__  k 

x_(t)  »  I  a,(t)e  (t)  =  e(t)a(t)  . 
u  iel  1  1 


Remark.  The  expression  for  x^( t)  as  well  as  the  following  equations 
for  x^(t)  (i  *=  1,2,...)  which  appear  in  the  construction  of  the 
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asymptotic  expansion  of  the  solution  will  be  considered  not  only  for 
t  *=  0,  u>  2 p,...  but  also  for  al 1  t  in  the  interval  [0,T]  .  This 
is  necessary  as  seen  already  in  the  equation  for  x^(t)  and  in  the 
determination  of  the  subsequent  functions  x^(t)  (i  =  1,2,...)  .  The 
asymptotic  representation  (43)  naturally  involves  only  values  of  x^(t) 
for  t  =  0,  p,  2u,  ...  corresponding  to  the  discrete  variation  of  t 
in  (41)  . 

For  TTqX(t)  we  obtain  the  system 

ttox(t+i)  =  B(o)rrox(T)  ,  nQx(o)  =  x°  -  xQ(o)  .  (46) 

In  contrast  to  (15)  the  system  (46)  is  a  linear  constant  coefficient 
difference  system  in  which  the  argument  t  varies  discretely  with  unit 
steps.  Its  solution  can  be  constructed  in  a  manner  completely  analogous 
to  the  construction  of  the  solution  of  the  constant  coefficient  differen¬ 
tial  system  (cf.  [7,21]),  namely 

k  m 

TT  x(t)  =  t  c.e  (0)  +  I  c  w  (T)u  (t)  . 

u  i=l  1  1  i=k+ 1  1 

The  w^(t)  appearing  here  are  somewhat  different  than  those  in  the 
expression  (17),  but  as  before,  its  components  are  polynomials  in  T 
and  the  u^(t)  are  not  exponentials  as  in  (16)  but  solutions  of  the  scalar 
difference  equations 

u^(t+1)  =  Xi(0)ui(T),  u..(0)  =  1  (i  =  k+l,...,m)  . 
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[If  h^,...,hn  is  a  Jordan  chain  of  vectors  corresponding  to  the  eigen¬ 
values  of  the  matrix  B  ,  then  there  are  n  linearly  independent  solutions 
x  (t)  =  wr(T)XT  ,  (r  =  l,...,n)  ,  of  the  system  x(t+1)  =  Rx(t)  ,  where 
(cf.  (17)) 


V  (t)  = 
r 


T(T-l) 
. r-1 


(r-1): 


t(t-I) . . . (r-r+; 
X  <~(r~2) ! 


h„  4  . . .  4  h 

2  r 


Thus, 

U.(T)  =  [X.(0)]T  (t  =  0,1,2,...)  , 

from  which  it  follows  by  Condition  I  that  ju^(r)|  decays  exponentially 

with  increasing  T  .  Taking  this  into  account  and  proceeding  as  in  §1 

we  set  c^  =  0  for  i  =  l,...,k  and  further,  with  the  aid  of  equation  (18) 

(as  in  §1)  we  determine  0^(0)  (i  =  l,...,k)  and  c^i  =  k+l,...,m)  . 

For  x^(t)  we  obtain  the  system: 

_  _  to  (t) 

X1(t)  *  Bftjx^t)  +  f(xQ(t),t,0) - -pj— -  , 

coinciding  with  (19).  The  solvability  condition  for  this  system  leads  to 
a  system  of  differential  equations  (20)  for  a(t),  as  described  in  §1. 

III.  Suppose  that  Conditions  IV  and  V  of  §1  are  satisfied. 

The  construction  of  further  terms  is  almost  exactly  as  in  §1.  An 
unimportant  difference  occurs  only  in  the  fact  that  in  place  of  (23) 
we  have  a  nonhomogene ous  difference  equation 
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TMc(t4  1)  =  B(0)TTiX(t)  +  TB' {0)TTqx(t) 

+  f(x0(o)  +  rrox(T),o,o)  -  f(xQ(o),o,o)  , 

and  in  the  expression  for  TT^x(t)  we  have  u^(t)  in  place  of 
exp(^.(0)r)  •  Such  differences  occur  in  subsequent  ^-functions,  but 
for  them  the  exponential  estimates 

|)Trix(r)||  <  cexp(-aT)  (t  =  0,1,2,...) 

hold.  Let  us  denote  by  Xn(t,p)  the  n-th  partial  sum  of  the 
series  (L3)  . 

Theorem  1.2.  Under  Conditions  I- III  there  exist  positive  constants  u., 
and  c  such  that  for  0  <  p  <  pQ  a  solution  x(t)  of  the  problem  f U 1) , 
(42)  exists  in  the  interval  [0,T],  is  unique  and  satisfies  the  inequality 

*\i,  1 

||x(t)  -  Xn(t,n)||  <  cp  for  t  =  0,  p,  2p,  ...  (t  <  T) 

Proof.  Let  us  set  ?(t)  =  x(t)  -  X^-^tjp)  »  then  for  F(t)  we  obtain 
the  system 

?(t+p)  =  B(t)?(t)  +  pfx(t,p)F(t)  4  G(F,t,p)  , 

?(0)  =  0  , 


where 
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fx(t,n)  =  fx(x0(t)  +  ^(t)  ,t,n)  and 


G(e.t,u)»  B(t)xml(t,n)  +  nf(xml(t,u)  +  5,t,n) 


ufx(t,n)?  -  xwl(t,u)  . 


Note  that  G(E,t,y.)  has  the  same  two  properties  as  in  §1  . 


Let  us  now  set 


w(t)  =  T'1(t-u)^(t) 


) 


where  u  has  (m-k) -components,  v k-components,  and  the  matrix  T(t) 
puts  B(t)  into  the  block- diagonal  form 


T'1(t)B(t)T(t)  = 


Here  is  the  (k  x  k) -identity  matrix  and  B^( t)  has  eigenvalues 

X^(t)  satisfying  the  inequality  J X t) j  <1  . 

For  u  and  v  we  obtain  the  system  of  equations 


u(t+u)  =  B^tjuft)  +  C11(t,n)u(t)  +  C12(t,fi)v(t) 

+  G2(u,v,t,^)  ,  (47) 

v(t+p.)  =  v(t)  +  C21(t,p.)u{t)  +  C22(t,n)v(t) 

+  G3(u,v,t,n)  , 

u( 0)  *  0,  v(0)  =  0  ,  (48) 


* 


4 

i 
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where 


C(t,n)  = 


cn(t,n)  c12(t,u) 


C21(t,n)  C22(t,u) 


=  T-1(t)B(t)[T(t^)  -  T(t)  ]  4  nT'1(t)f  (t,4)  *  T(t-u) 


satisfies  the  estimate 


C(tfn)  =  ©(h)  , 


G^UjVjt.u)  = 


G2(u,7,t,H)\  /uv 

1  =  T  (t)G(T(t-n)  (  J  ,  t, 

G-,(u,v,t,u)/ 


has  the  same  two  properties  as  the  function  G(?,t,p.)  • 

Let  us  denote  by  U(t,s,u)  and  V(t,s,n)  the  corresponding  matrix 
solutions  of  the  homogeneous  difference  problems 

U(tf h,SjH)  ~  B^(t)U(t,s,u)  (t  =  5yS+H»  •••  )  » 

U(s,s,p.)  —  > 

and 

V(Um.,s,h)  =  (ER  +  C22(t,u))V(t,r.,p.)(t  »  s,&4h,...)  , 


V(s,s,h)  =  E.  . 


Since 


||^k+  C22(t,H,))V(t,s,M,)||  <  (l4cn)!|V{  t,s,n))| 
by  (49) ,  it  follows  that 


3^ 


II V(t, s,m.)II  <  k(l+cn)(t's)^  <  k(l+cn)T/u  <  c  , 

while  U(t,s,|i)  ,  by  virtue  of  the  fact  that  the  eigenvalues  of  the  matrix 
B-^(t)  have  modulus  less  than  one,  satisfies  the  inequality 

|IU(t,s,n)l!  <  C  exp(  )  for  t  =  s,s+m-  ,  ...  (50) 

M* 

(cf.  [13;  Lemma  6.2))  . 

Using  the  matrices  U(t,s,u)  and  V(t,s,q.)  we  convert  the 
system  (L7),  (L8)  into  the  equivalent  system  of  equations  (t  =  in, t = 0,1,2 , 

C 

u(t+u)  =  E  U(t+u,(Ul-i)u,u)  [Ci;L(t-in,4)u(t-iu)  +  C12(t-iu,u)v(t-iu) 
i=0 

+  G-(u(t-iii),v(t-iu),t-iu,4)  )  , 

t 

v(un)  =  E  V(Un,((rtl-i)n,u)  [c  (t-iM.,p.)u(t-i^) 
i=0 

+  G3(u(t-ip.)  ,v(t-i|A)  ,t-iu,^)  ]  . 

[The  solution  of  the  difference  problem 

z(tx-h)  ■  A(t)z(t)  +  b(t)  (t  =  0,h,2h, . . .) ,  z(0)  =  z°  , 


can  be  written  in  the  form  (t  =  th,  L  =  0,1,2,...) 


z(t+h)  «=  $(t+h,0)z°  +  E  *(bth,(t»-l-i)h)b(t-ih)  , 

i=0 


where 


t(Uh,c)  =  A(t)$(t,s)  (t  =  s,£Hh,...)  ,  $(s,s)  «=  b:  .  ) 


3(? 


By  applying  the  method  of  successive  approximation;  to  this  system  and 
using  the  two  properties  of  and  G^  as  well  as  the  estimates  (<’•‘9), 

(50),  it  is  not  difficult  to  prove  (analogously  to  [13,  §1'])  that 
a  solution  exists,  is  unique  and  satisfies  the  estimates  u(t)  =  fcfu  ) , 

TVL  1 

v(t)  =  &(u  )  .  Thus  Theorem  1.2  is  established. 

2 

Remark.  Suppose  that  the  stepsize  in  (11)  is  equal  to  u  and  not 
Mi  ,  that  is,  the  order  of  smallness  of  the  stepsize  is  greater  than  the 
coefficient  of  the  nonlinear  term  f(x(t),t,p.)  : 

x(t>fM2)  =  B(t)x(t)  +  uf  (x(  t)  ,t,u)  •  (hi) 

Then  besides  the  functions  x(t,u)  and  TxtTjU)  (t  ~  t/u)  the  boundary 

p 

function  Px(s,u)  (s  »  t/u  )  will  occur  in  the  conditions  determining 
the  asymptotic  solution  of  the  initial  value  problem,  that  is,  the 
indicated  difference  problem  is  similar  in  the  sense  of  asymptotic 
behavior  to  the  differential  problem  which  we  discussed  in  Kemark  3  of 
Subsection  U  in  §1.  Detailed  considerations  are  given  in  f‘;  '  . 


53  Applications 

1.  Difference  Formulas  for  the  Numerical  Integration  of  Differential 
Equations.  Difference  formulas  for  the  numerical  integration  of  initial 
value  problems  for  the  scalar  differential  equation 
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y’  =  f(y,t),  y(0)  =  y°  (0  <  t  <  T) 
can  be  written  in  the  form  (cf.  [2]) 


(52) 


+  at-iyk4t-i 4 


Vk  =  h(etfk+t+ 


(53) 


where  y  =  y( ih) ,  h  is  the  stepsise,  f  =  f(y^,ih)  , 

stants  defined  by  certain  equations,  one  of  which  is 
it  follows  that  the  characteristic  equation 


and  a.,  (3. 
1  1 

l 

la.  *  0  . 
i=0  1 


are  con- 

Hence, 


4 


,4-1 

at-lX  4 


+  on 


0 


corresponding  to  (53)  always  has  the  root  X  =  1  .  If  we  write  (53)  in 
the  form  of  a  system  of  difference  equations  of  the  type  {>(1)  (the  role 
of  p  is  played  by  h  )  then  the  matrix  U( t)  of  the  resulting  system 
has  the  eigenvalue  X  =  1  . 

Let  us  consider  in  more  detail  one  difference  scheme  of  the  type  (53) 


yk+2  '  2yk+l+  2yk  H^fk+l"3V 

(this  is  called  a  formula  of  extrapolation  type,  of  second  order  (t  =  2) 
with  two  steps).  Since  4  «  2  ,  to  use  this  formula  we  must  prescribe 
the  initial  values  y(0)  and  y(h)  ,  where  y(h)  is  nearly  equal  to 

y(0)  by  virtue  of  the  smallness  of  h  . 

0  0  12  2 

Let  us  set  y(0)  =  y  ,  y(h)  =  y  4  hy  4  h  y  4  . . .  ,  where 

1  2 

y  ,y  , are  certain  coefficients.  Moreover  let  us  set  t  =  kh, 
y  =  y(kh)  =  x^t),  yk+1  =  x^Uh)  =  x2(t)  .  Then  wc  obtain  a  system 
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of  two  difference  equations  of  the  form  (Ul)  for  x(t) 
(t  =  0,h,2h,...),  namely 


x1(Uh)  =  x2(t)  , 


x2(t+h)  =  •'§x1(t)  +  fX2^ 

+  Jh  [5f(x2(t),tfh)  -  3f(Xl(t),t) 


(5«0 


with  initial  conditions 


x.,  ( 0)  =  y°,  x  ( 0)  =  y°  +  hy 1  +  h2y2  +  . . .  . 


(55) 


/ 0 

The  matrix  3(t)  =  I  ^  ^  J  of  system  (51*)  has  the  eigenvalues  =  1  , 

V  2  2/ 

X2  =  |  ,  and  so  Condition  I  of  §2  holds.  A  difference  between  the 

problems  (5*0,  (55)  and  (41) ,  (42)  is  that  the  initial  value  x2(0) 

depends  on  h  ;  however,  this  dependence  introduces  only  an  insignificant 

change  in  the  initial  conditions  for  the  coefficients  of  the  asymptotic 

expansion  (43),  and  the  algorithm  for  the  construction  of  the  asymptotic 

solution  as  well  as  all  estimates  remain  valid.  _ 

/  xio(t)\ 

The  system  of  equations  corresponding  to  x  (t)  =[  1  will 

V  *20  <*y 

have  the  form  xQ(t)  =  B(t)xQ(t)  ,  that  is, 

X10(t)  =  X20(t)>  X20(t)  =  -l*10{t)  +  lX20(t)  ’ 

whence,  xQ(t)  =  a(t)e1  =  a(t)(^)  where  e^  =  (^)  is  the  eigenvector 
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of  the  matrix  B(t)  corresponding  to  =  1  and  a(t)  is  an  arbitrary 
scalar  function. 

f  Vi(T)\ 

The  system  of  equations  for  TT  x(t)  =  f  ]  has  the  form 

Vv2<tv 

(cf.  (U6))  n,Qx( T+ 1)  =  B( 0)tTqX( t)  ,  that  is, 


TT0Xl(T+1)  =  V2(T)  » 

TT0x2(T+l)  =  •|tT()x1(t)  +  |tT0x2(T)  ’  (56) 

with  initial  conditions 

Vi(°)  =  y°  -  a(°)»  ttqx2(o)  =  y°  -  a(o)  .  (57) 

The  general  solution  of  (56)  can  be  written  as 

V(T)  “  Clei  +  c2^  e2  *  cl^l'1  +  C2^  ^1/2) 

(t  =  0,1,2, • • •  )  , 

where  e2  =  (2/2)  the  eigenvec't°r  B(t)  corresponding  to  \ ^  =  l/2 
We  set  c^  =  0  in  order  that  TTqx(t)  —  0  as  T  -»  00  ,  and  substituting 
TTqX(t)  into  the  initial  conditions  (57)  we  obtain 

cg  ■  y°  -  a(0),  ~c2  =  y°  -  a(o)  •, 


whence,  c_  =  0  and  a(0)  =  yU  .  Thus  TT  x(t)  =  0 


Xi-(t) 


The  system  of  equations  corresponding  to  x.(t)  =[  I  is 

(  xa(t)  J 
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x-^t)  =  x21(t)  -  a'(t)  , 

x21(t)  =  -■|x11(t)  +  2^21^  '  a'  (t)  +  -|f(a(t)  ,t)  .  (58) 

Its  solvability  condition  implies  the  equation 

a'  =  f(a,t)  .  (59) 

The  initial  condition  for  a(t)  was  already  determined:  a( 0)  =  y°  . 

Given  that  we  can  find  a(t)  ,  the  construction  of  the  zeroth  approximation 
XQ(t,h)  =  xQ(t)  +  TTqx(t)  =  xQ(t)  is  complete.  If  the  construction  of 
the  asymptotic  solution  is  continued  then  boundary  functions  appear  in 
the  following  terms,  even  in  the  term  containing  the  arbitrary  y1  . 

We  note  that  the  zeroth  approximation  X1Q(t,h)  for  x^(t)  is 
equal  to  a(t),  as  we  would  expect  from  considering  the  exact  solution  of 
the  initial  value  problem  (52) .  The  solution  obtained  from  the  difference 
scheme  therefore  generally  differs  from  the  exact  solution  by  a  term  of 
order  h  . 

If  y^  were  not  prescribed  arbitrarily,  but  rather  so  that 
0  12  2 

y(h)  =  y  +hy  +hy  +  ...  differs  from  the  exact  solution  of  (52) 
by  a  term  of  order  h  (for  this  it  is  necessary  to  set  y\0)  =  f(y  ,0)) 

then  TT^x(t)  and  x^1(t)  are  equal  to  zero,  while  x^(t)  ,  found  from 
the  system  of  difference  equations,  differs  from  the  exact  solution  of 
the  problem  (52)  by  a  term  of  order  h  .  In  fact,  by  virtue  of  (59) > 
the  system  (58)  assumes  the  form 


uo 


xu(t)  =  x21(t)  -  a'(t)  , 

x21(t)  =  -|x1;L(t)  +  I *21^)  ■■|a,  » 


whence 

The  equation  for  the  scalar  function  @(t)  is  obtained  from  the  solvability 

/  *12* 

condition  for  the  system  corresponding  to  x  (t)  =  /  _  ],  namely 

V  x22(t) 

*l2(t)  =  X00(t)  -  ^x"n(t)  -I*,(t)  , 


22 


2  10 


11' 


*22(t)  =  *  |*22(t)  -|x20(t)  -  x^(t) 

+  i[5f  (t)x2i(t)  +  5ft(t)  -  3f  xn(t)] 


It  can  be  writ'  en  in  the  form 


=  f  (a(t) ,t)p  . 

Since  TTqx(t)  =  0  ,  the  equation  for  TT^x(t)  coincides  with  the 
equation  for  TTqx(t)  >  and  consequently. 


because  d1  =  0  in  order  that  T^xCt}  -  0  as  t  -®  .  Substituting 
for  TT^x(t)  and  x^(t)  into  the  initial  conditions 

TriX1^  0)  +  0)  =  °  »  ^i3^*0)  +  *21*  =  yl  » 


we  obtain 
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d2  +  p(0)  =  0  , 

|d2+3(0)  =  y1  -  a' (0)  =  y1  -  f(y°,0)  . 

Hence  if  y^"  =  f(y°,0),  d  =3(0)  ®  0  ,  so  TT^x(r)  =  0  ,  0(t)  =  0  and 
x11 (t)  =  0  ,  as  stated  above. 

Remark.  For  the  numerical  integration  of  the  singularly  perturbed 
initial  value  problem 

=  f(y,t),  y(0,p)  =  y°(0  <  t  <  T) 

2 

we  naturally  choose  a  stepsize  smaller  than  p  ,  for  example  h  =  p 

Then  if  we  write  the  difference  formula  in  the  form  (53)  with  step 
2 

h  =  p  ,  we  obtain  an  equation  which  reduces  to  a  system  of  the  form 

(51). 

2.  Markov  Chains.  A.  A  continuous- time,  homogeneous  Markov 
pro cess  with  finitely  many  states  m  can  be  described  by  the  system 

~  =  AP,  P(0)  =  P°  ,  (60) 

CLT 

where  P(t)  is  a  vector  with  components  p^(t) ,. . . ,pm(T)  (p^t)  is  the 

probability  of  being  in  the  i-th  state  at  time  T  ) ,  and  A  is  a  matrix 

m 

with  constant  elements  a  ,  satisfying  E  a.  .  =  0  .  It  is  known 

ij  i=l 

(cf. ,  for  example,  [20])  that  such  a  matrix  A  has  4  =  0  as  an  eigen¬ 
value  with  as  many  linearly  independent  eigenvectors  as  the  multiplicity 
of  the  root  4=0.  We  will  assume  that  the  other  eigenvalues  satisfy 
the  condition  Re4^<  0  (a  so-called  proper  chain)  . 
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In  the  study  of  Markov  chains  we  are  interested,  in  particular,  in 
the  limiting  behavior  of  the  probabilities  p^(t)  as  T  —  00  .  For  the 
system  (60)  the  limiting  values  for  p  -•  0  and  fixed  t  =  Tp  /  0  are 
found  as  the  components  of  the  solution  P(t,p)  of 

liff  =  AP,  P(0,p)  =  P°  . 

Suppose  that  X  =  0  has  the  eigenvectors  .  Then  it  follows 

-L  v  K 

from  $1  that 

k 

lim  P(t)  =  lim  P(t,p)  =  P  =  Z  a.(t)e.  =  ea(t)  , 

T-®  p-0  i=l  1  1 

where  a(t)  is  determined  from  equation  (20),  which  for  f  =  0  has 

the  form 

(ge)f^  =  0  ; 

whence,  ^  =  0  ,  that  is,  a(t)  s  a(0)  .  The  initial  value  a(0)  can 
be  determined  from  system  (37),  which  hero  assumes  the  (matrix)  form 

(ge)a(0)  =  P°  . 


In  particular,  if  the  chain  is  regular  (X  =  0  is  a  simple  root)  then 


1,  e,  =  e  *= 


,  while  it  is  possible  to  take  g,  =  g  =  (!,...,!) 


by  virtue  of  the  condition  Z  a 

i=0 


'ik 


0  .  Thus 


m  -  m  m  m  n 

a(o)  =  tpJ/te  .1/te,  since  Z  pT  =  1 
i*l  1  i=l  1  i=l  1  i=l  1 


I 


and 


p0  =  e/te  . 

U  i=l  1 

Therefore  the  limiting  probability  is  here  independent  of  the  initial 
state  P^  . 

B.  A  discrete  Markov  chain  can  be  described  by  the  difference 


system 


P(sH-l)  *=  BP(s),  P( 0)  =  P°  , 


where  s  =  0,1,2,...  is  the  number  of  trials.  Now  £  b  =  1  ,  and  so 

i=l  1K 

the  matrix  B  has  X  =  1  as  an  eigenvalue.  As  before  we  will 
assume  that  this  is  a  proper  chain,  that  is,  each  eigenvalue  different 
from  one  satisfies  |X_J  <1  . 

The  limiting  value  as  s-*®  of  the  probability  P(s)  coincides 
with  the  limiting  value  as  4  ~  0  for  fixed  t  ^  0  °f  the  solution 
P(t,p.)  of  the  system  (t  =  sn) 


P(t+4,4)  =  BP(t,n)  . 


It  follows  from  §2  that  the  limiting  value  of  P(t,u)  as  u  0  can  be 
described  as  in  case  A,  so  for  a  regular  chain 


where  e 


lim  P(s)  *  lim  P(t,ji)  =  P  =  e/  I  e . 
s—  4-0  U  i=l  1 

is  the  eigenvector  of  B  corresponding  to  X  =  1  . 


C.  Passage  from  a  discrete  process  to  a  continuous  one.  We  consider 
the  Markov  chain  (6l)  where  we  take  as  independent  variables  the  instants 
of  time  t  at  which  the  trials  occur.  We  suppose  that  the  trials  are 
separated  from  each  other  by  a  small  time  interval  p.  .  Introducing  t 
we  obtain  a  system  which  agrees  formally  with  (62),  but  for  which  we 
naturally  assume  that  the  transition  probabilities  bik  are  small  (of  order  u) 
for  i  4=  k  ,  but  nearly  one  for  i  =  k  ,  that  is, 

bik  =  ^aik(^ (i  +  k) ’  bii  s  1  +  *a..(u) >  aik(n)  ^  aik  » 

as  p.  -*  0  .  Then  (62)  assumes  the  form 

p±( -fc+M.)  =  Pi(t)  +  p,(ailp1(t)+  ...  +  aimPm( t) )  (i  =  l,...,m).  (63) 

This  system  is  of  the  type  (41)  where  B  =  Em  has  X  =  1  as  an  eigen¬ 
value  of  multiplicity  m  to  which  there  correspond  m  eigenvectors 


These  will  be  the  vectors  g.(j  =  l,...,m)  . 

It  follows  from  §2  that 

lim  P(t,p.)  =  ea(t)  =  E  a(t)  =  a(t)  , 

where  a(t)  is  determined  from  the  differential  equation  (20),  which 
in  this  case  assumes  the  form 


(A  is  the  matrix  with  elements  a.,  )  with  initial  condition 

iK 

0(0)  =  g(0)p°  =  P°  . 

Thus  the  limiting  value  P(t,u)  satisfies  a  system  of  differential 
equations  of  the  type  (60)  . 


Chapter  2 


Nonlinear  Singularly  Ferturbed  Equations 
in  the  Critical  Case:  Initial  Value 
Problems 

$1  Statement  of  the  Problem  and  Auxiliary  Results 
1.  Statement  of  the  Problem.  In  this  chapter  we  consider  problems 
analogous  to  those  of  Chapter  1;  however,  now  the  nonlinearity  in  the 
right-hand  side  of  the  equations  is  not  necessarily  small: 

^  9  *(*>*&)>  0  <  t  <T  ,  (1) 

x(0,p)  =  x°  ,  (2) 

where  x  is  an  m-dimensional  vector.  As  in  Chapter  1  we  will  consider 
cases  when  the  reduced  equation  does  not  have  an  isolated  solution,  but 
a  family  of  solutions.  Consequently  the  same  questions  arise  as  there. 
Under  what  conditions  will  the  solution  of  the  problem  (1),  (2)  converge 
as  pi  *•  0  to  one  of  the  solutions  of  this  family,  and  in  particular,  to 
which  one?  How  can  one  construct  the  asymptotic  expansion  of  the  solu¬ 
tion  to  an  arbitrary  order  in  p,  ,  uniformly  valid  for  t  in  [0,T|? 

The  conditions,  under  which  the  asymptotic  expansion  will  be  constructed, 
are  numbered  I,  II,  ... 

I.  Suppose  that  the  function  F(x,t,y)  is  sufficiently  smooth  in 
the  domain  D(x,t,p)  =  D(x,t)  X  [0,4q]  ,  where  D(x,t)  is  a  domain  in 
(x,t) -space  and  Uq 


is  a  constant. 


II.  Suppose  that  the  reduced  equation 


bl 


0  =  F(x,t,0) 

has  for  each  t  in  [0,T]  a  family  of  solutions 
x  =  cp(t,a1}...  ,0^)  =  <c(t,a)  , 

where  cp(t,a)  is  a  well-defined  function  of  t  and  the  arbitrary  para¬ 
meters  a1J***Jak  which  satisfies  the  following  conditions  in  the  domain 
D(t,a)  «  [0,T]  X  D(a)  : 

1)  the  function  cp(t,Gt)  is  stiff iciently  smooth; 

2)  the  rank  of  the  matrix  cpQ(t,a)  =  is  equal  to  k  , 

the  number  of  parameters. 

From  Condition  II  it  follows  that  F(<p(t,a)  ,t,0)  s  o  for  (t,a) 
in  D(t,a)  .  Differentiating  this  identity  with  respect  to  a  we  obtain 

Fx(cp(t,a)  ,t,0)cpa(t,a)  =  0  for  (t,a)  in  D(t,a)  . 

This  implies  that  the  matrix  Fx(cp(t,a)  ,t,0)  has  the  eigenvalue  X(t,a)  $  0 
and  that  the  columns  of  the  matrix  <PQ(t,ci)  are  eigenvectors  corresponding 
to  X  =  0  .  By  virtue  of  condition  2)  of  II  these  columns  are  linearly 
independent  since  the  multiplicity  of  X  =  0  is  not  less  than  k  . 

III.  Suppose  that  the  multiplicity  of  the  eigenvalue  X  =  0  is 
exactly  equal  to  k  and  that  the  remaining  eigenvalues  Xi(t,a)  of  the 
matrix  Fx(cp(t,a)  ,t,0)  satisfy 
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Re  X  i( t,a)  <  0  (3) 

in  D(t,a)  . 

Under  Conditions  I- III  and  certain  others  which  will  be  stated  below 
we  will  construct  the  formal  asymptotic  expansion  of  the  solution.  In 
§3  we  estimate  the  remainder  term,  while  in  the  next  subsection  we  obtain 
a  number  of  auxiliary  results  which  follow  from  analogous  results  in 
[13,  Subsection  3]  . 

2.  Auxiliary  Results.  A.  Stability  Manifolds.  An  equation  which 
will  play  an  important  role  in  the  construction  of  boundary  functions  is 

=  F(cp(0,a)  +  X,0,0)  ,  (U) 

where  a  is  a  parameter.  For  any  a  in  D(a)  this  equation  has  the 
rest  point  x  =  0  .  By  Condition  III  the  characteristic  equation  corre¬ 
sponding  to  this  rest  point,  det(Fx(co(0,a)  ,0, 0)  -  XE^)  =  0  ,  has  the  root 
X  =  0  of  multiplicity  k  and  (m-k)  roots  satisfying  condition  (3)  . 
Therefore,  the  rest  point  x  =  0  is  not  asymptotically  stable  in  the 
sense  of  Lyapunov,  that  is,  a  solution  with  initial  value  arbitrarily 
close  to  the  rest  point  will  not  necessarily  converge  to  it  as  T  -*  00  . 
However,  if  we  prescribe  special  initial  conditions,  then  the  solution 
will  converge  exponentially  to  the  rest  point  x  =  0  as  T  -  ®  .  Put 
precisely,  we  have 

Igmraa^.  In  a  sufficiently  small  neighborhood  of  the  point  x  =  0 
there  exists  an  (m-k) -dimensional  manifold  o>(a)  such  that  if  the 


initial  values  x(0)  belong  to  aid)  ,  then  one  can  find  positive 
constants  y  and  a  such  that  for  r  >  0  the  solution  x(t)  satisfies 
the  inequality 

I|x(t)||  <  v  exp(-  a t)  .  (5) 

Proof.  We  linearize  the  right-hand  side  of  (U)  with  respect  to  x 
and  write  the  system  (U)  as 

Jjf  =  A(a)x  +  G(x,a)  ,  (6) 

where 

A(a)  =  Fx(cp(0,a)  ,0,0)  and  G(x,a)  =  F(®(o,a)  -4  X,0,0)  -  A(a)x  . 

The  function  G(x,a)  has  the  following  two  important  properties: 

1.  G(0,a)  =  F(co(o,a)  ,0,0)  =  o  . 

2.  For  any  e  >  0  there  exists  a  5  >  0  (depending  on  e  and 
possibly  on  a)  such  that  if  |lx^||  <  6  and  || x^l I  <6  then 

||G(x1>a)  -  G(Xgfa)||  <  e|lx1  -  x^l  . 

This  inequality  is  established  by  elementary  means  using  Taylor's  formula, 
and  it  shows  that  for  sufficiently  small  ||xl|  G(x,a)  is  a  contraction 

operator. 

As  noted  above  the  matrix  A(a)  has  \  =  0  as  an  eigenvalue  of 
multiplicity  k  and  m-k  eigenvalues  ^(0,0)  which  satisfy  condition  (3) 
Thus  there  exists  a  matrix  B(a)  ,  as  smooth  as  A(a)  (cf.  [29])  ,  whicn 
reduces  A(a)  to  the  block-diagonal  form 


B_1(a)A(a)B(a)  = 


c(a)  o 


(7) 


^  *)• 

where  the  ((m-k)  x  (m-k) ) -matrix  C(a)  has  the  stable  eigenvalues 
Xi(0,a)  satisfying  condition  (3)  . 

Let  us  make  the  change  of  variables 

x  =  B(a)  Q  , 

where  u  has  (m-k)-  and  v  k-components  .  For  u  and  v  we  obtain 

=  C(o)u  +  G1(u,v,a) ,  =  G2(u,v,o)  ,  (8) 

in  which  G1  and  G 2  are  blocks  of  vectors  B  ^(ar)G(B(a)  (^)  ,a)  satisfying 
the  same  two  properties  as  G(x,a) 

We  introduce  the  system  of  integral  equations 
0  T  , 

u(t)  =  u(T,a)u  +  J*  u(T,a)u"  (s,a)G1(u(s),v(s),a)ds  , 

0 

and  (9) 

T 

v(t)  =  J  G2(u(s),v(s),a)ds  , 

QD 

where  U(T,a)  is  a  fundamental  matrix  of  the  system  ~  =  C(a)u(U(0,a)  =E  .  ) 

aT  m-k 

and  u°  is  an  arbitrary  constant  vector.  The  matrix  U(r,a)  satisfies 
the  inequality  ||U(T,a)U  ^(s,a)||  <M  exp(-H(T-s))  (where  the  positive  con¬ 
stants  M  and  h  can  depend  on  a  )  .  Every  solution  of  system  (9)  also 
satisfies  system  (8). 

We  apply  the  method  of  successive  approximations  to  (9) ,  replacing  u 
and  v  in  the  right-hand  side  of  (9)  by  ur  and  vr  ,  and  in  the  left-hand 


i 


f.l 

side  by  u  .  and  v  .  .  Taking  u_  =  v.  =  0  we  obtain  u't'  =  U'T,a  : 
'  tw- 1  W-i  0  U  -L 

and  v^(t)  =  0  ;  whence  , 

||u1(t)||  <  Mexp(-nT)  I'u0!!  <  M !l u°ll  exp(-cT)  ,  (10) 

where  a  is  any  number  in  the  interval  (0,k)  . 

Let  us  set  p  =  max(  m/(h-c) ,  l/c)  and  choose  e  >  0  so  small  that 

20e  =  q  <  1  .  Corresponding  to  this  £  is  a  certain  number  &  define  1 

by  the  second  property  of  the  functions  G.  and  G„  .  Let  us  pick  c  >  C 

-L  c 

so  that  the  inequality  (M/2) (l/(l-q)  +  l)p  <  5  holds.  Nov  consider 
those  u^  with  |!u°||  <  o  •  Using  (10)  and  Property  1  of  G-^  and  G0 
we  obtain 

T 

II  u  (t)  -  u  (t)||  <  r  Mexp(-n(T-s) )eM|  u)|exp( -os)  ds 
"  0 

<  (M/2)  q||  u^l  exp(-oT)  , 

CO 

||v2(t)  -  v1(t)||  <  [  £Kllu^!exp(-crs)  ds 

1  T 

<  (JV’2)qJlu^|!exp(-aT) 

It  is  easy  to  show  that  for  n  >  1 

llu^i^)  "  Un^T^  —  (H^2)<in|l^lexp(-oT)  , 
lluml(T)||  <  (M/2)(qn+  qn  1  +  ...  +  2)|!u°||exp(-rrT) 

<  6  exp(-cT) 


5? 


and  likewise  for  v  (t) .  The  uniform  convergence  in  r  of  the 
n 

successive  approximations  follows,  and  this  proves  the  existence  of  a 
solution  satisfying 

|Iu(t)||  <5exp(-ct)  ,  |iv(f)|!  <Sexp(-oT) 


Hence,  (5)  follows  directly. 

The  desired  manifold  a-(a)  has  the  form 

to(a)  =  [x:  x  =  B(a)(UQ)  ,  u  =  u°  ,  v°  =  J  G  (u(s)  ,v(s)  ,a)ds  , 

V  CD 

l!u°||  <  p}  . 

1 

This  concludes  the  proof  of  Lemma  1. 

If  we  consider  the  linear  approximation  for  system  (U) ,  that  is, 
if  we  set  G(x,a)  =  0  in  (6),  then  system  (8)  assumes  the  form 


du 

dT 


=  C(a)u 


and  consequently,  to  obtain  a  solution  which  converges  exponentially  to 
zero  as  T  -*  <*>  ,  it  is  necessary  to  take  v(0)  =  0  (whence  v(t)  =  0) 
and  u(0)  =  u°  ,  where  u°  is  arbitrary.  Let  us  denote  by  z  and  y 
the  upper  and  lower  blocks  of  the  vector  x  corresponding  to  the  dimen¬ 
sions  m-k  and  k  and  by  the  corresponding  blocks  of  the 

matrix  B(a)  *,  from  the  equation  (  )  =  B(a)(Q)  we  obtain  z  =  B^(a)u 
and  y  =  B21(r*)u  . 
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IV.  Suppose  that  det  B^Cci)  /  0  for  a  in  D(a)  . 

Then  for  the  linear  approximation  the  manifold  <c{a)  can  be 
written  as 

y  =  B21(a)B^(a)z  .  ( ID 

B.  Extension  of  the  Stability  Manifold. 

The  statement  of  Lemma  1  has  a  local  character.  If  we  continue  the  trajec¬ 
tories  originating  in  u:( a)  in  the  negative  T-direction  we  obtain  an 
extended  manifold  (denoted  by  0(a))  which  has  the  same  property  as  j(d), 
that  is,  trajectories  starting  in  0(a)  at  T  =  0  remain  in  0(a) 
for  T  >  0  and  converge  exponentially  to  the  rest  point  x  =  0  as 
T  —  <*>  .  In  some  cases  we  can  construct  0(a)  in  an  explicit  form 
(cf.§§4  and  5  below).  We  will  assume  in  the  present  chapter  that  0(a) 
admits  of  the  following  analytical  representation. 

V.  Suppose  that  in  some  domain  D(z,a)  the  manifold  0(a)  can  be 
represented  as 

y  =  P(z,a)  ,  (12) 

where  P(z,a)  is  a  sufficiently  smooth  function. 

Indeed,  the  definition  (12)  is  an  identity  along  trajectories 
x(t)  which  converge  exponentially  to  the  rest  point  x  =  0  as  T  —  00  , 
that  is,  the  manifold  0(a)  consists  of  such  trajectories.  Therefore, 
along  any  such  trajectory  ^  ^  ,  and  setting  ^z,0^  =  H(z,d)  , 


Therefore 
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It  is  clear  that  the  last  sum  taken  along  a  trajectory  x(t)  in  C.ia; 
is  equal  to  dH^/dT  .  Likewise,  along  the  indicated  trajectory  the 
equation  (-|^  F^)  =  cLH/dt  holds,  and  consequently,  from  (15)  we  have 

f  -  (F21*  F2£H'  -  H<FU*  F12K)  <16> 

along  any  trajectory  x(t)  in  Cl(ct)  . 


C.  The  Variational  System  on  the  Stability  Manifold.  We  consider 
now  a  nonhomogeneous  system  of  equations  whose  homogeneous  part  is  the 
variational  system  for  (4),  namely 

^  -  Fx(t)4  +  *<T)  ,  (IT) 

where  F^T)  =  F^Mo.a)  +  x(t),0,0),  x(t)  6  0(a)  for  some  constant 
vector  a  in  D(a) ,  and  i|i(t)  is  a  certain  function.  The  upper  and 
lower  blocks  of  A  having  dimensions  m-k  and  k  are  denoted  by 
A1  and  A2  ,  while  the  corresponding  blocks  of  ( t)  are  denoted  by 

^(t)  and  *2(t)  . 

Lemma  2 .  The  change  of  variables 


Ll  =  81»  L2  =  H(T)8l  +  8 2  > 


(18) 


where  H(t)  =  H(z(t),qO  and  z(t)  is  the  upper  block  of  x(t) 
forms  the  system  (17)  into  the  form 


dBi 

"d^  =  an(r)51  +  a12(T)52  +  ^(t)  , 


trans- 


(19) 
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a22(T)62+  (*2(t)  -  H(t)*1(t)  )  , 


where 


an(T)  =  Fn(T)  4  F12(t)H(t),  a12(T)  =  F12(t)  , 
a22(T)  -  F22(t)  -  H(t)F12(t)  ; 


(20) 


here  the  F. .  are  the  blocks  of  the  matrix  F  (t)  ,  as  in  (15)  . 
ij  x 

The  essence  of  the  lemma  is  that  by  this  change  of  variables  the 
equation  for  &2  can  be  separated  from  the  equation  for  5^  .  To  prove 
Lemma  2  it  is  necessary  to  write  system  (17)  in  block  form,  make  the 
change  of  variables,  and  use  the  identity  (16) . 

Suppose  now  that  the  nonhomogeneous  term  \|i(t)  in  system  (17)  is 
Fx(T)cpa(0,a)  .  Then  the  second  equation  of  (19)  assumes  the  form 
d& 

~dT  -  a22(T)&2  +  (Fx(t)coa(o,a)  ]2  -  H(T)  [Fx(T)cDa(0,a)  )1  (21) 

(here  and  below  the  indices  1  and  2  denote  the  upper  m-k  and  the  lower 
k  rows  of  the  indicated  matrix) . 

A  particular  solution  of  system  (17)  for  the  given  nonhomogeneous 
term  is  clearly  A  =  -cpa(0,a)  •  The  corresponding  particular  solution  of 
(21)  is 

62  =  H(T)  CcPa(°,a)  ]1  *  tcPa(°’a)  ]2  5  R(T’a)  • 


Thus  we  have  proved 
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Lemma  3-  The  matrix  6^  =  R(T,a)  is  a  solution  of  equation  (21) 
satisfying  the  initial  condition 

Bg(0)  =  H(0)[tpa(0,a)  ]x  -  [«Pa(0,a)]2  =  R(o,a)  . 


D.  We  now  obtain  a  number  of  important  results  for  the  matrices 
an(  t)  and  a22^T^  (20))*  T^t  us  denote  by  H(°»)  the  limiting  value 

of  H(t)  =  H(z(t) ,a)  as  t  -  ®  .  Since  (11)  is  the  linear  approximation 
for  (12)  we  have  that 


§f(o,a)  =  h(-)  =  B21(a)B^(a) 


(22) 


Let  us  denote  by  F.  .(“)  the  limiting  value  of  F.  .(t)  .  Then  F.  .(«■) 

ij 

are  the  blocks  of  the  matrix  F(«)  =  Fx(to(0,or)  ,0,0)  =  A(a)  appearing  in 

equation  (7) .  We  write  (7)  as 

/ c(a)  0’ 

F  (-)B(a)  =  B(a) ( 

\  0 

and  by  equating  blocks  with  index  11  we  obtain 


Fu(®)Bu(a)  +  F12(®)B21(a)  =  Bu(a)C(a)  . 

On  account  of  (20)  and  (22),  it  follows  that 

Fu(“)  +  F12(»)h(»)  =  a11(»)  =  Bu(a)c(a)B^(a)  . 

Clearly  the  eigenvalues  of  a^*)  coincide  with  those  of  C(a)  •, 
they  are  the  ^(0,0)  which  satisfy  condition  (3).  Thus,  the  fundamental 
matrix  $(t)  of  d&j/dT  =  a-^-r)^  ($(0)  =  Fm  k)  satisfies  (cf.  [13, 
(3.78)]) 
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||$(t)$  ^(s)||  <  c  exp( -k(t-s) )  for  0  <  s  <  T 


(23) 


Analogously  we  can  write  (7)  as 


B 


/  C(a)  0\ 

(a)F  (-)  =[  J  B  (or) 

\  °  °) 


and  by  equating  the  blocks  with  index  22  we  obtain 


+  (B_1(Q())22F22(»)  =  0 


Since 

(B'1(a))21=  -(B_1(a))22B2;L(a)B^(a)  (2h) 

and 

det(B~1(a)  )22  /  0  (  25) 

(cf.  [201  or  [13,  (U.55^  1) , 

f22(-i  -  E21(a)B^(«)r12(«)  =  j-22(«»>  -  h(»u12(») 

*  Sgg(® '  =  0  . 

Now  F  (t)  =  F..(<p(0,a)  +  x(  t) ,0,0)  converges  exponentially  to  V..(<*>) 
ij 

as  T  by  virtue  of  the  exponential  convergence  of  x(t)  to  zero, 

and  so  it  follows  that 

<  cexp(-KT)  for  T>0  .  (26' 

Lemma  U.  The  fundamental  matrix  Y(t)  for  db^dT  =  a22(T)52('F('0>  =  F  ) 

satisfies: 


,  i 
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(1)  Y(»)  =  limY(T)  exists^ 

T-® 

(2)  detY(®)/0  •, 


(3)  H*(t)  -  V (®)ll  <  cexp(-HT) 


Proof.  Consider  the  matrix  integral  equation 


¥(T)=Elc+J  a  (s)f(s)ds  . 


Applying  successive  approximations  to  it,  we  consider  the  sequence 


Vi(T)  "W  Ws>Vs>ds  -  VT>  ■-  E„ 


By  virtue  of  (26)  and  (28), 


-  Yq(t)||  <  J  jla22(s)|!ds  <  ( c/h) exp( -rt)  , 

T 

M  »  •  ^ 

l!Tg(t)  -  Y1(*r)|l  <  /  lla22(s)||||Y1(s)  -  Y0(s)|lds 


<  ~^(  c/k  )  2exp(  -  2kT  )  , 


II^Ct)  '  <  ■j^-(c/K)nexp(-nHT) 


<  — j-  for  a  =  c/h  ,  n  =  1,2, . 
™"  n  • 


Hence,  it  follows  that 


V’)  =  V 
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converges  as  n  **  ®  (uniformly  with  respect  to  t  >  0)  to  a  matrix 
'Y(t)  which  satisfies  equation  (27)  and  consequently,  dY/dr  =  a^li^Y 
and  Y (®)  =  E  •  It  follows  also  that  detY(O)  /  0  ,  for  otherwise 
detY(r)  =  0  ,  contradicting  the  relation  Y(®)  =  •  If  we  now  set 

Y(t)  =  Y(t)Y_1(0)  , 

we  obtain  a  fundamental  matrix  for  which  properties  (1),  (2)  and  (3) 
of  Lemma  4  are  satisfied. 


E.  Let  us  turn  now  to  the  matrix  R(t,o')  of  Lemma  3. 


Lemma  *3 .  detR(“,a)/0 


Proof.  Denote  the  upper  m-k  and  the  lower  k  rows  of  the  matrix 
1  -  ”2 

(29) 


B*1(a)cp  (0,a)  by  h.  and  h  respectively,  that  is, 


cr 


(4)  =  B’1(0,)'pa(0>a>  • 

CD* 


Substituting  a>a(0,ot)  =  B(a) 

-1 


into  Fx(®)ccfJ(0,a)  =  0  ,  multiplying 


on  the  left  by  B  (a)  and  taking  account  of  (7) ,  we  obtain 


(?  :)  ©  ■  -  • 

whence,  it  follows  that  C(a)h^  =  0  ,  so  that  h^  =  0  .  Since  the  rank  of 
(p^(°,a)  is  equal  to  k  ,  deth^/O  .  From  (29)  we  then  obtain 

hg  =  (B_1(a))21Ccpa(o,a)  ]1  +  (B"1(a))22(tpQ(o,a)  )2  . 
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By  virtue  of  (24)  and  (22), 

h2  =  -(B_1(a))22{H(®)  [cpa(0,a)  11  -  [coa(0,rr)  ]2) 

=  -(B*1(a))22R(®,a)  . 

Since  det(B_1(a) )  22  /  0  (cf.  (25)),  det  B(®  ,a)  /  0  . 

§2  An  Algorithm  for  the  Construction  of  the  Asymptotic  Expansion 
of  the  Solution  of  the  Initial  Value  Problem 

The  asymptotic  expansion  of  the  solution  of  problem  (1),  (2)  will 
be  constructed  in  the  form 

x(t,n)  =  x(t,u)  +  tTx(t,u)  (t  =  t/u)  ,  (30) 

where 

x(t,p.)  =  xQ(t)  +  ux-Jt)  4  ...  4  ^nxn(t)+  ...  , 

TTx(t,u)  =  TT  x(«T)  +  (iff  x(t)  +  ...  +  Hn,TnX(T)  4  ... 

By  substituting  (30)  into  (1)  and  representing  the  function  F  in  the 
form  F  =  F  4  ttf  just  as  in  Chapter  1,  we  obtain  a  sequence  of  equations 
for  the  determination  of  x^(t)  and  h\x(T)  (i  =  0,1,2,...) 

For  x^(t)  we  have 

F(xQ(t),0,0)  =  0  . 

By  virtue  of  Condition  I  solutions  of  this  equation  can  be  written  in 
the  form 


xQ(t)  =  ®(t,a(t))  , 


(3D 
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where  o(t)  is  an  as  yet  arbitrary  k-dimensional  vector  function. 
For  TTqx(t)  we  obtain 

dJT  x 

=  F(co(0,a(0) )  +  ttox,o,o)  , 


which  after  the  substitution  x  =  TT^x  is  made  coincides  exactly  with 
equation  (4)  for  a  =  a(0)  .  The  initial  condition  for  ^TqX(t)  is 
obtained  after  substituting  (30)  into  (2)  as 


rrQx(0)  =  x°  -  co(0,a(0)) 


z°  -cc1(o,a(  0)) 
y°  -  cc2(o,a(o) ) 


where 


is  the  upper  (m-k) -dimensional  block  of  x 


is  the 


lower  k-dimensional  block,  and  cp^,  co2  are  the  analogous  blocks  of  ®  . 
Both  the  equation  and  the  initial  condition  for  Pqx(t)  then  involve 
the  as  yet  arbitrary  vector  a( 0)  .  Let  us  use  this  arbitrariness  to 
guarantee  the  exponential  convergence  of  ^q^t)  to  zero  as  t  -  ®  . 

For  this  it  is  sufficient  to  require  ^hat  TT^O)  belong  to  Q(a(0))  , 
that  is,  that  tTqx( 0)  satisfies  (12).  This  gives 


y°  -  tp2(o,a(o))  =  p(z°  -  rai(o,a(o))  ,a(o))  . 


(32) 


Equation  (32)  is  a  k-dimensional  vector  equation  for  the  k  components 
of  the  vector  a(0)  . 


VI.  Suppose  that  equation  (32)  has  a  solution  a(0)  =  a°  . 

Taking  a(0)  *  a°  ,  TTqx(t)  belongs  to  n(a°)  for  T  >  0  ,  that  is, 
the  blocks  nQz(r)  and  tT^t)  of  ^^(t)  satisfy  ^^(t)  =  P(ttqz(t)  ,a°) 
and  consequently,  for  T  >  0 
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HtT0x(t)||  <  c  exp(-KT) 

Remark.  The  functional  determinant  (corresponding  to  (32))  is 

dcpo(0,a)  n  dcp.CC.O') 

&(a)  =  det(  - ^ - +  H(z  -  ^(0,00,0:) - ^ - 

PaU°  -  cp1(o,a),a))  • 

If  we  denote  H(tT0z(t)  ,a)  by  H(t)  ,H(t)  [c£>a(  0,a)  ]]_  -  [®a(0,a)]2  by 
R(T,a)  (as  in  §1)  ,  and  P(TT  z(«r)  ,0:)  by  P(T,a)  then 

A(a)  =  det[R(o,a)  -  Pa(o,a)  ]  . 

For  many  singular  perturbation  situations  analogous  to  this 
(cf.,  for  example,  [13,  §13,  Condition  III])  it  is  assumed  that  the  cor¬ 
responding  functional  determinant  is  nonzero.  In  the  present  problem  this 
requirement  is  unnecessary  since  it  is  not  difficult  to  prove  that 
PQ(T,a)  =  R(t, Qt)  -  Y(t)Y"1(«)R(®,0')  (Y(t)  being  the  fundamental  matrix 
from  Lemma  4).  Hence,  it  follows  that  A(a)  =  detf  (®)R(OD,a)/o  by 
virtue  of  Lemmas  4  and  5* 

Thus  the  function  is  completely  determined,  although  a(t) 

occurs  in  the  expression  (31)  for  x^(t)  •  We  only  know  a( 0)  =  a°  . 

The  function  a(t)  is  determined  completely  from  a  solvability  condition 
in  the  equation  for  x^(t)  . 

The  equation  for  x^(t)  has  the  form 


6h 


d x  _  _  _ 

-dT-  V*0(t),t,0)*1*  y*0<t),t,0) 

or 

Fx(®(t,of(t) )  ,t,0)x1  =  cca(t,a(t)  )j^  +  crt(t,a(t))  -  (33) 

Fu(o(t,a(t)),t,o)  . 

The  determinant  of  this  linear  algebraic  system  of  equations  is  equal 
to  zero.  For  the  solvability  of  this  system  it  is  necessary  and  suffi¬ 
cient  that  the  right-hand  side  be  orthogonal  to  the  eigenvectors 

* 

g  (t,a(t))(j  =  l,...,k)  of  the  adjoint  matrix  F  (co(  t,a(t) )  ,t,0)  cor- 
J  x 

responding  to  the  eigenvalue  X  =  0  .  Let  us  denote  by  g(t,a(t))  the 

(k * m) -matrix  whose  rows  are  the  g.(t,a(t))  .  Then  the  orthogonality 

J 

condition  can  be  written  as 


(g(t,a(t))®a(t,a(t)))-~  + 


(3M 


(g(t,a(t))  [®t(t,a(t))  -  F^(cp(t,a(t))  ,t,o)  ])  =  o  , 

where  (gco^)  denotes  the  ( k  X  k) -product  of  g  and  •,  analogous  meaning 
is  given  to  the  other  terms  in  (34).  As  noted  in  Chapter  1  detCgco^)  /  0  , 
and  so  (34)  can  be  solved  for  do/dt  : 


(35) 


VII.  Suppose  that  equation  (35)  together  with  the  initial  condition 
a(0)  =  a°  has  a  solution  a  =  a(t)  for  t  in  [0,T)  that  belongs  to 
D(a)  there,  where  D(a)  is  the  domain  in  Condition  II  . 


Through  a(t)  we  completely  determine  the  zero-th  term,  of  the 
approximation.  Let  us  introduce  the  curve  L  consisting  of  the  two 
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where  ^r(r)  is  known  and  such  that  IUr(T)ll  <  cexp(-*T)  . 

The  initial  condition  is 

fT^O)  =  -x1(0)  =  -cr>a(  0)0(0)  -  3^(0) 

Thus,  an  as  yet  arbitrary  vector  3(0)  appears  in  the  equation  and 
in  the  initial  condition  for  ^(t)  .  We  use  this  arbitrariness  to  guaran¬ 
tee  that  ^x(t)  decreases  exponentially  as  T  ■*  ®  .  Let  us  denote 
the  upper  and  lower  blocks  of  ^  x  by  and  f^y  ,  and  let 

V  =  &1  ,  TT^y  =  H(t)61  +  &2  , 

where  H(t)  =  H(TT0z(T),a)  .  Lemma  2  implies  the  equations 

d&, 

-jjr  =  a11(T)61+  a12(r)62  +  [Fx(T)cpa(0M13(O)  +  tj/*)  , 

d6-  _ 

a22(T)52  +  UFx(T)coa(0)]2  -  H(T)[Fx(T)coa(0)  1^3(0) 

+  tVg(T)  -  H(t)^(t))  (39) 

with  initial  conditions 

6^0)  =  -(®a(  0)^(0)  -  z1(0)  , 

62(0)  =  {H(0)[5o(0)  Jx  -  f5a(0)]2J3(0)  +  {H( 0)2^0)  -  yx(0)}  •  (40) 

Using  Lemma  3  and  introducing  &2  =  {H(0)z^  (0)  -  yx(0)}  we  have 
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62(t)  =  R(t,u3)0(O)  +  y(T)5° 


+  f  y(T)y-1(s)[t0(s)  -  H(s)+.(s)]ds  . 

0 

By  requiring  that  b^T)  -*  0  as  T  ®  we  obtain 

00 

R(“,a°)P(0)  =  -y(«){6°+  J  y  ( s)  ( s)  -  H(s)ijr  (s)]dsj  .  (^< 

^  c  0  ^ 

By  virtue  of  Lemma  5  this  equation  is  uniquely  solvable  for  3  =  3(0)  . 
Substituting  into  (Ll)  and  using  the  exponential  decay  of  y(r)  and  H(t) 
we  have 


|62(t)||  <  cexp(-KT)  for  T  >  0 


Since  F  (T)tp  (0)  =  [F  (t)  -  F  (0)  ]cp  (0)  satisfies  the  same  exponen- 

X  U  X  X  Q. 

tial  estimate  , 


db.. 

IF  =  aH(T)6l  +  *1(T)  » 

where  ||^(t)||  <  cexp(-HT)  .  Thus 


5x(t)  =  $(t)&1(0)  +  J  $(▼)$  1(s)f1(s)ds  , 

where  the  fundamental  matrix  $(t)  satisfies  (23),  and  so 


||61(t)||  <  cexp(-HT)+  J  c exp(-*(T-s)  )c exp(-ns)ds 
<  c  exp(-nT)  . 


(45) 


?x(t)  x2  +  2(x1#  Pxx(t)Xl)  +  F  (t)  +  2  Fpu(t). 


Here  (x.  ,  F  (t)x.)  is  the  vector  whose  components  are  the  scalar  products 

J-  XX  x 

_  _ 0  _  ®  2 _ P  4  4  _ 4 _ 4  _ 

<x. ,  F  x.>  «  £  OF  /3x  3xJVt)x  xJ  (£=l,...,m).  Substituting  x.  from 

i,j=l 

formula  (36)  and  writing  the  solvability  condition  for  equation  (45)  in  a 
form  analogous  to  the  condition  (34),  we  obtain  the  equation 


d6 

dt 


f1(6,t). 


(46) 


At  first  glance  it  may  seem  that  f^(8,t)  depends  quadratically  on  8(t) 

as  a  result  of  the  term  (x, ,F  (t)x,)  in  (45),  However  this  is  not  the  case, 

1  xx  1 

since 


<  gj(t,a(t)) ,  (<J,a(t)  8(t) ,  Fxx(t)  <t>a(t)  8(t))  >  =  0  (47) 


for  j  *  In  order  to  verify  (47)  we  differentiate  the  identity 

F(<J>(t,a)  ,t,0)  =  0  twice  with  respect  to  the  components  a  and  a  of  the 

p  q 

vector  a(p,q  =  l,.,.,k)  and  obtain 


(2±-  v  1  +  F  0  v  ,  o 

l3a.  *  xx  3a  '  x  3a  3a 


Then,  forming  the  scalar  product  with  g,(t,a)  and  noting  that  <  g. ,  F 
2  J  J  x 

3 

g---g—  >  *  0  since  F  *g,  ■  0,  we  see  that 

p  q  3 


<  *j ’  P«^  >  "  0 
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for  j  ,  P.  <?  *  1 . k. 

Thus  equation  (46)  is  linear,  that  is,  f^(6,t)  =  A(t)  6 ( t )  +  B^(t),  where 
expressions  for  A(t)  and  B^(t)  are  obtained  easily  from  (45).  By  virtue  of 
this  linearity,  there  is  a  unique  solution  of  (46)  which  exists  in  [ 0 , T] 
and  satisfies  8(0)  =  g^. 

Thus  the  term  of  order  y  in  the  asymptotic  expansion  is  completely 
determined.  The  determination  of  the  successive  terms  of  the  expansion  proceeds 
in  a  manner  analogous  to  that  of  x^(t)  and  H^x(t).  At  the  i-th  step  an 
arbitrary  function  (say  y(t))  appears  in  the  expression  for  x^(t).  First 
we  determine  y(0)  from  the  condition  that  IT  x(t)  -*•  0  as  T  -*■  °°:  the  equation 
for  y(0)  is  of  the  same  type  as  (42),  with  det  R(°°,a^)  4  0,  Then  from  the 

solvability  condition  in  the  equation  for  x^+^(t)  we  obtain  an  equation  for 
y(t)  like  (46),  namely 

-  A(t)  y  +  B^t), 

which  determines  y(t)  uniquely  in  [0,T]. 

Thus  it  is  possible,  under  Conditions  I -VIII,  to  construct  arbitrarily 
many  terms  of  the  series  (30). 


§3  An  Estimate  of  the  Remainder  Term 

As  in  Chapter  1  let  us  first  make  a  more  precise  formulation  of 
Condition  I  concerning  the  smoothness  of  the  function  F(x,t,u)  •  We  note 
that  it  is  possible  to  take  an  arbitrary  6-tube  of  the  curve  L 
(cf.  VIII)  to  be  the  domain  D(x,t)  . 

I.  Suppose  that  the  function  F(x,t,p.)  has  continuous  partial 
derivatives  with  respect  to  each  argument  up  to  order  (m-2)  inclusive 
in  the  domain  D(x,t,|j,)  =  D(x,t)  X  [0,uQ] 

Let  us  set 

k  .  _ 

x.(t,u.)  =  l/(x  (t)  +  n  x(t))  . 

i=0  1 


Theorem  2.  Under  Conditions  i-vill  there  exist  Dositive  constants 
and  c  such  that  for  0  <  p,  <  the  solution  x(t,u)  of  the 
problem  (1),  (2)  exists  in  the  interval  [0,7],  is  unique  and  satisfies 
the  inequality 

||x(t,u)  -  Xn(t,n)||  <  ciiml  (0  <  t  <  T)  . 

Proof.  Substituting  x  =  X  .  +  E  into  equation  (1)  we  obtain 

± 

an  equation 

^dt  *  +  >  (U8) 


where 
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and 


GU,t,n)  =  F(Xml(t,n)  + 


dX^  ,(t,u) 


The  function  G(£,t,u)  has  the  following  two  important  properties 
which  can  be  established  just  as  easily  as  in  [13,  §10,  Subsection  1] 

namely 


1.  G(0,t,n)  =  Mn**2)  . 

2.  If  II  ?j,  (t,ii)||  <  c^p,2  and  ||?2(t,y.)||  <  c^2  for  0  <  t  <  T 

and  0  <  p,  <  p,1  (for  some  constants  c^  and  p  ) ,  then  there  exist 
constants  cQ  and  pQ  <  such  that  for  0  <  t  <  T  and  0  <  u  <  pQ 

Il0(51*t,p.)  -  G(?2,t,n)||  <  c^2znax  tl^  -  ^l  (t9' 


(for  pQ  note  the  remark  in  Subsection  3,  §1,  Chapter  1)  .  In  conformity 
with  Chapter  1  G(^,t,p)  is  a  contraction  operator  with  contraction 
coefficient  of  order  Q(p2)  for  e  =  ®(p2) 

We  now  introduce  the  change  of  variables 

?  -  T(t)  (")  , 


where 


T_1(t)Fx(t)T(t)  = 


an(t} 


> 


0 


0 
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for  Fx(t)  =  F  (Xg(t) ,t,0)  .  Here  the  eigenvalues  X^(t,a(t))  of  the 
((m-k)  x  (m-k) ) -matrix  a11 (t)  satisfy  condition  (3)  •  Then 


T_1(t)F  (t,n)T(t)  =1 


ai;L(t,g.) 


a21(t,p.) 


^(t.y.) 

a22(t,u) 


) 


so 

lla^t,^)  -  a11( t)||  <  c(exp(-Kt/p.)+  4)  , 


while  the  other  blocks  a^(t,p.)  satisfy 
II aik(t>M-)H  <  c(exp(-Ht/p.)  +  4)  . 

The  system  for  u  and  v  has  the  form 

V> =  a11(t,4)u+  al2(t,4)v  -  nb^tju  -  Ml>12(t)v+  (T"1C)1  , 


»» =  a21(t,4)u+  a22(t,n)v  -  4b21(t)u  -  4*>22(t)v  +  (T_1G)2  , 


21' 


22 


21' 


where  the  b^(t)  are  the  blocks  of  T*^(t) 


dT(t) 


dt 


Note  that 


Aik(t,4)  =  a^(t,4)  -  4bi^(t)  satisfy  the  same  inequalities  as  the 

aik(t,4)  • 


Suppose  now  that  U(t,s,4)  and  V(t,s,4)  are  the  fundamental 
matrices  of  the  homogeneous  systems 


'■‘dt  =  Au^*^0  (U(S»S>M.)  =  Em_k)  , 
^dt  =  A22^»^)v  (v(s»s>^)  ■  Er) 


(50) 
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By  the  properties  of  a^(t),  A.^(t,p)  and  A02(t,p)  ,  these  fundamental 
matrices  satisfy  for  0  <  s  <  t  <  T  and  0  <  p  <  p^ 

||U(t,s,p)||  <  c  exp(-K(t-s)/p) ,  I! V( t , s ,n) IS  <  c  . 

Using  the  fundamental  matrix  V(t,s,p)  and  the  trivial  initial 
values  of  u  and  v  (as  well  as  of  5),  we  can  express  the  second  equation 
in  (50)  as  the  integral  equation 

t 

v(t,p)  =  f  K  (t,s,p)u(s,p)ds  +  Q^(u,v,t,p)  .  (5I1 

0  * 

Here  the  kernel 

K2(t,s,p)  =  p  1V(t,s,p)A21(s,p) 
clearly  satisfies  the  inequality 

||K2(t,s,p)||  <  c [p  ^xpC-Ks/p)  +  11  ,  (52) 

while  the  integral  operator 

-1  1  -1 

Qg(u,v,t,p)  =  P  J  V(t,s,p)(T"  G)  ds 
c  0  c- 

by  the  two  properties  of  satisfies  the  estimate  Q^(0909t9\i)  = 

1 

)  and  is  a  contraction  operator  with  an  order  ft(p)  contraction 

2 

coefficient  for  u  and  v  of  order  (t(p  )  . 

By  substituting  (51)  into  the  first  equation  of  (50)  and  using 
U(t,s,p),  we  obtain  the  integral  equation 
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u(t,u)  =  P  K  (t,s,n)u(s,n)ds  +  Q  (u,v,t,u)  , 

0 

where  the  kernel 

,  t 

K,(t,s,n)  =  H  f  U(t,pJu)A12(p,^)K2(p,s,u)dp 

X  s 

satisfies  the  same  inequality  as  K2(tjS,p)  (of.  (52))  >  while  the 
integral  operator 

,  t  -1 

Q  (u,v,t,u)  *|i  J  U(t,s,p) [Al2(s,p)Q2(u,v,s,p)  +  (T  G)1]ds 

1  °  0 

has  the  same  two  properties  as  Q2(u,v,t,|i)  . 

Let  us  denote  by  R(t,s,p.)  the  resolvent  kernel  of  K^tjS,^)  .  It 
satisfies  the  same  estimates  as  the  kernel  itself.  We  can  express  (53) 
as  the  equivalent  equation 

t 

u(t,n)  =  Q1(u,v,tJ^)  +  J*  R(t,s,vi)Q1(u,v,s,p.)ds 

0 

h  S1(u,v,t,n)  , 

where  the  integral  operator  S1(u,v,t,Mi)  has  the  same  two  properties 
as  Q^UjVjt.p.)  . 

Substituting  (54)  into  (51)  we  obtain 
t 

v(t,n)  =  J  K  (t,s,p.)S1(u,v,s,ti)ds  +  Q  (u,v,t,p) 

0 

E  s2(u,v,t,p.)  , 


7i; 

where  S2  has  the  properties  of  •  Therefore  we  can  apply  the  method 

of  successive  approximations  to  the  system  (5t),  (55)  (with  uq  =  vc  = 

and  easily  show  as  in  {13,  §10]  that  for  sufficiently  small  u  a  solution 

u(t,|i),  v(t,p.)  exists  in  the  interval  [0,T]  ,  is  unique  and  satisfies 

the  estimates  u(t,u)  =  ©(m,™”^)  and  v(t,M.)  =  ©(ii'"*''")  .  Hence,  it  follows 

directly  that  E;(t,|i)  =  x(t,p.)  -  X^^tjU.)  =  ©(u™^)  >  so  x(t,u)  -  Xn(t,u'  - 

rv^  1 

&(p,  )  and  this  proves  the  theorem. 

§4  Special  Cases 

1.  Consider  the  system  of  equations 
-  A(y,t)  2  +  B(y,t) 

(0<t<T)  (56) 

^  =  C(y,t)z  +  D(y,t) 

with  the  infinitely  large  (as  4  -*  0)  initial  condition 

z(0,p,)  =  z°/\i  ,  y(0,n)  =  y°  .  (57) 

In  the  special  case  that  z  and  y  are  scalar  functions  with  C(y,t)  =  1 
and  D(y,t)  »  0  this  problem  was  considered  in  detail  in  [13,  §l6]  • 

Suppose  now  that  in  the  system  (56)  z  is  an  (m-k) -dimensional  vcct  r 

and  y  a  k-dimensional  vector.  Let  us  introduce  in  place  of  z  the 
function  nz  (which  we  will  again  denote  by  z  )  ,  then  (56) >  (57)  takes 
the  form 
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U-fl  =  A(y,t)  z  +  uB(y,t)  , 

|i ^  =  C(y,t)z  +  nD(y,t)  , 
z(0,u)  =  z°  ,  y(0,u)  =  y°  . 

This  is  a  problem  of  the  form  (1),  (2)  for  the  m-dimensional  vector 
function  x  =  (  )  .  The  reduced  system 

A(y,t)z  =  0  ,  C(y,t)z  =  0 

has  the  family  of  solutions 

z  =  0  ,  y  =  a  s  cp(t,a) 


The  matrix  Fx(co(t,a)  ,t,0)  can  be  written  in  block  form  as 


and,  consequently,  Condition  III  is  satisfied  provided  the  eigenvalues 
X^tjCt)  (i  =  1,  ...k)  of  the  matrix  A(a,t)  satisfy  the  inequality  (3), 
that  is, 


Re Xi(t,a)  < 0 

The  matrix  q?a(t,a)  ,  consisting  of  the  eigenvectors  corresponding  to 
X  s  0  ,  now  has  the  form 


(58) 


77 


and  the  system  (U)  can  be  written  as 

■||  =  A(oy,0)  2  ,  =  C(Ofy,0)z  . 


(59' 


Suppose  that  z  and  y  are  scalars,  that  is,  m-k  =  k  =  1  .  Then 
condition  (58)  reduces  to  A(y, t)  <  0  .  If  we  assume  that  C(y,t)  is  of 
constant  sign  then  from  (59)  we  obtain  an  explicit  representation  of  the 
manifold  fl(a)  ,  namely 


z  =  ry  Moy.jO) 

Jq  C(oy,0) 


dy 


5 


which  agrees  exactly  with  the  formula  in  (12)  . 

If  C(y,t)  s  1  and  D(y,t)  =  0  (this  case  was  discussed  in 
[13,  §16]),  then  =  (°)  ,  g  =  (1,-A),  ,  and  equation  (3U) 

assumes  the  form 


(60) 


A(a,t)^+  B(a,t)  =  0  . 

This  agrees  precisely  with  equation  (^*385)  of  [I3,§l6]  ,  with  a(t) 
playing  the  role  of  yQ(t) .  The  initial  condition  for  a(t)  is 
determined  from  equation  (32)  which,  in  the  present  case,  through  (60) 
can  be  written  as 

0  7°-a(°) 

z  =  J  A(a(0)  +  y,0)dy  * 

0 

hence,  we  obtain 

0  X°  0 

ZU  =  J  A(T1 , 0) dLTl  or  zU  +  J  ACn,0)dTl  =  0  . 

a(0)  y° 


The  succeeding  equations  agree  exactly  with  the  formulation  in  equation 
(4.395)  of  [13,  §16]  ,  and  from  them  we  can  determine  y  ( 0)  .  Thus  a 
coincides  with  y  (0)  ,  that  is,  the  formal  construction  of  §2  reduces 
to  the  results  obtained  in  [13,  §16]  • 


where  co(t,a)  is  a  z-root  of  the  equation  K(z,a,t)  =  0  .  It  is  possible 
to  develop  the  construction  of  the  asymptotic  solution  of  the  problem  (63) , 
(64)  by  the  scheme  of  §2  ,  which  after  a  finite  number  of  calculations 
gives  the  same  result  as  in  [13,  Chapter  3]  . 
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Thus  we  are  led  to  the  system  of  differential  equations 


t  Y..w  (t  =  . 

i=l  lC  1 


(65' 


Under  actual  conditions  the  rate  constants  differ  widely  from  each 
other.  This  property  can  be  expressed  by  means  of  a  small  parameter  u  • 


+ 

r 

+ 

Suppose  then 

that  k .  =  a  k.  (i=l,...,n<n)  . 

i  1  ’  ’ 

Then 

we  have 

w. 

1 

—  Cx.n  rr.  —  3... 

, +  ll  1m  -  ll 

=  uw.  =  k.  x,  ...  x  -  k.  x.  ... 

1  1  1  m  1  1 

&im 

X 

m 

(i  =  1, 

So  the  system  (65)  can  be  written  as 
dx  n  _  n 

»■?’  'Vi"!  Vu“i 

“1  i.»l 

Setting  p,  =  0  we  obtain  the  reduced  system 


.  ,m) 


(66) 


0  =  t  Y ..  W  (t  =  1,  • .  •  ,m) 
i-1  Xl  1 


(67) 


In  practice  it  often  happens  that  system  (67)  has  a  family  of  solutions 
which  depend  on  one  or  more  arbitrary  parameters,  and  thus,  the  problem 
reduces  to  a  singularly  perturbed  equation  (66)  in  the  critical  case. 

One  method  for  determining  approximate  solutions  of  the  equations  of 
chemical  kinetics  containing  a  small  parameter  is  known  in  physical 
chemistry  as  the  method  of  quasi- stationary  concentrations  of  Semenov - 
Bodenstein.  A  number  of  works  are  devoted  to  questions  involving  the 
mathematical  justification  of  this  method  (that  is,  to  a  justification 
of  the  passage  to  the  limit  as  p  -  0)  ;  cf.,  for  example,  [17,18]  . 
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We  now  discuss  an  example  of  an  actual,  chemical  reaction  for  which  the 
calculations  can  be  carried  out  using  the  asymptotic  methods  presented  in 
this  chapter.  This  system  is 


dx^  dx 

-jr  -  -klxi  *  kix2  >  -sr  *  klxi  - k:x 


1-2  -  k3x2x3  -  • 


dx  +  dx^ 

“dt  =  -k3X2X3  ’  "dt~  =  -kUX2XU  * 


(68) 


[Such  a  system  occurs  in  investigations  of  the  reaction  kinetics  of 
organometallic  compounds  and  was  proposed  by  A.N.  Kashinym,  a  colleague 
of  ours  in  the  chemistry  department  of  Moscow  State  University.  ]  The  rate 
constants  have  the  orders  of  magnitude 


k*~10  ,  k^~  109  ,  k*  ~  1010  ,  k*  ~  108  ,  k”  =  kj^  =  0  . 

Dividing  each  of  these  equations  by  k^  and  making  the  substitutions 
H  =  l/k^  ,  a  =  k*  ,  b  =  k~/k^  and  =  c  ,  we  obtain 


dx. 


dxr 


~=  -nax1+bx2  ,  ^  axx  -  bx2 -cx2x3  -  x^  , 


dx, 


1 - _cx2x3  -Vu 


(69) 


The  reduced  system 


0  =  bx2  ,  0  =  -bx2  -  c  x^  -  x2x^  , 
0  *  -c^  ,  0  =  -x£x4 


has  a  family  of  solutions  depending  on  three  arbitrary  parameters,  namely 
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x,  = 


=  °1  »  X2  =  °  »  X 3  =  a2  »  xl+  =  a3  ‘ 


'70) 


Here  the  matrix  F  (co(t,a)  ,t,0)  is 


-b  -  ca2  - 


-  ca 


0 

0 

0 

0 


o\ 

0 

0 

°l 


and  its  eigenvalues  are  X^  =  =  0  and  X^  =  -b  -  ca^  -  .  Since 

b  >  0  ,  c  >  0  and  a 2  ,  a are  nonnegative  (which  makes  sense  physically) , 
it  follows  that  X^  <  0  .  Thus  Conditions  I-III  of  §1  are  satisfied. 

The  system  (U)  now  has  the  form 


-cx2(a2+x3) 


-  cx2(a2+x3) 

’  Sr  =  “  x2 


*  X2(a3+xU) 

<a3+V  . 


y 


This  system  is  sufficiently  simple  that  it  can  be  integrated,  and  we  obtain 
for  the  manifold  Q(a)  the  exact  representation 


C  Ai 

x2  *  "  X1  +  (exP(_f  Xx)  ■  +  (exp(--~ )- 1  )a3  , 

c  X1 

x3  =  (exp(-^x1)  -l)a2  ,  xh  =  (exp( -  — )  -  1)<*3  . 


(71) 


Thujas  in  (12),  the  lower  block  (consisting of  three  components)  of  the 
vector  x  is  expressed  in  terms  of  the  upper  block  (consisting  of  one 
component) .  By  the  same  token  Condition  V  is  satisfied.  An  elementary 
argument  verifies  Condition  IV  . 


_ ISFk. 
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Suppose  now  that  the  system  (69)  is  furnished  with  the  initial  condi¬ 
tion  x(0,|i)  =  .  Then  the  vector  equation  (32)  assumes  the  form 

x°  =  -(x^-aJO))  +  (exp[--^(x°-  0^(0))  ]  -  i}a„(o) 

+  {exp[--g  (x°  -  0^(0) )  ]  -  l}a^(0)  ,  (72' 

x°-a2(o)  =  {exp[--g (x° - 0^(0) )  ]  -  l}a2(o)  , 

x°-a3(0)  =  {exp[--g(x°-a1(0))  ]  -  l}a^(0)  . 

By  setting  t  =  x°  -  0^(0)  we  can  determine  a2(0)  and  a^(0)  in  terms 
of  t  ,  namely 

a2(0)  =  x°exp(*^t)  ,  0^(0)  =  x°exp(*»^t)  . 

Substituting  into  the  first  equation  of  (72)  we  obtain  an  equation  for  t 

0  /C^v  0  ,  1^,  0  0  0 

t  +  x^  exp(#»—  t)  +  x^  exp(»—  t)  =  x^  +  x^  -  x2  . 

Elementary  considerations  show  that  this  equation  has  a  unique  solution 
for  x°  >  0  and  x°  >  0  .  Thus  0^(0),  a2(0)  and  0^(0)  are  uniquely 
determined  from  (72),  that  is,  Condition  VI  holds. 

One  can  also  write  equations  for  a^(t),  a2(t)  and  o^(t)  .  In 
the  present  case  it  is  a  matter  of  integrating  by  quadratures.  Thus  one 
can  determine  x1Q(t)  =  a^t),  x2Q(t)  =  °>  =  a2^  311(1  xi,0(t)  • 

a^t)  . 

The  determination  of  ttqx(t)  reduces  to  the  integration  of  the 


scalar  equation 


s* 


dVi 

dT 


=  b(_TT0xi  +  texp(-|TTAx1)  -  l]aA(0) 


b  or 


fexp(  -  --TT^x^)  -  l)a^(O)} 


(73) 


with  tTqx1(  0)  =  x°  -  0^(0)  by  quadratures.  After  determining  ^x^t)  the 
remaining  functions  h^x^x)  ( i  =  2,3,M  are  found  by  means  of  the  equation 
for  fi(a)  (cf.  (71))  . 

Using  the  scheme  of  §2  we  can  also  construct  the  successive  terms 
of  the  asymptotic  expansion. 


2.  Equations  of  a  Nonequilibrium  Gas.  The  following  equations  are 

valid  for  a  spatially  homogeneous  gas  with  a  distribution  of  velocities  at 
equilibrium: 

dn 

U-^r=  £iQ(n,T)  T  M.£i;L(n,T)  ( i  =  1, . . .  ,N)  ,  (7M 

f=-fVn»T>  • 

Here  (i  =  1, ...,N)  denotes  the  density  of  those  particles  with 
internal  energy  ei  and  T  is  the  translational  temperature  (25]  . 

To  these  are  added  certain  conditions  n^(0,n)  =  n^  and  T(0,u)  =  T°  . 
Ei0(n,T)  characterizes  the  change  in  n^  as  a  result  of  exchanges  of 
energy  in  collisions,  while  E^(n,T)  characterizes  the  change  in  n^ 
as  a  result  of  the  transfer  of  internal  energy  to  the  energy  of  trans¬ 
lational  motion.  The  small  parameter  \i  signifies  that  the  transfer 
of  internal  energy  to  translational  energy  is  considerably  less  likely 
than  the  exchange  of  internal  energy  as  a  result  of  collisions,  l-'or 
particles  of  equal  mass  we  have 
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Ci0(n,T> 

tu'n,T) 


*  «kt(T)  (Vi  •  niV  > 

k,C,m 

Ae=0 


=  Z  p“(T)n.n  -  P?v(T)n  n  ) 
,  ,  k£  k  L  1m i  m 

KjVtin 

Ae/O 


,kfc. 


and 


Vn-T»  -  1  ei|pS(T|Vi  -  pki>v»> 

k,i,m,i 


Here  Ae  =  e.  +  e  -  e  -  e  ,  while  Q^T)  and  P^f(T)  are  the 

probability  of  the  exchange  of  internal  energy  in  collisions  and  the 

probability  of  the  transfer  of  internal  energy  to  translational  energy. 

The  reduced  equation  £  (n,  T)  =  0  has  solutions  obtained  from 

the  condition  that  n.  n,  =  n.n  •,  whence,  taking  note  of  the  fact  that 

k  l  i  m 

Ae  =  e^^  +  effi  -  efe  -  e^  =  0  we  have  that  n^  =  aexp(pe^)  for  arbitrary 
parameters  a  and  p  .  The  Boltzmann  distribution,  in  which  -l/p 
denotes  the  internal  temperature  and  which  depends  on  t  ,  is  found 
by  means  of  the  following  approximations  which  agree  with  the  general 
rules  stated  above  in  §2  . 

The  system  of  equations  for  the  h^-functions  has  the  form 

dTToni 

-  Cio(n(0)  +  V’  T(0))  (i  =  • 


The  law  of  conservation  of  particles  and  the  law  of  conservation  of 
energy  can  themselves  be  represented  by  two  first  integrals  of  this 
system,  namely 
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Z  (ru(0)  +  Tt^n^)  =  const.  =  Z  n^O) 

and 

Z  e^(n^(0)  +  TT^n^)  =  const.  =  ^  eji^(O)  , 
and  consequently,  the  equations 


Z  TT  n  =  0 
i  01 


Z  e.TT  n. 
i  1 01 


0 


furnish  a  2-dimensional  stability  manifold,  while  the  equations 


Z  [n?  -  a(0)exp(p(0)e  ) 1  =  0  ,  Z  e  [n?  -  a(0)exp(p(0)e. ) ]  =  0 
i  1  i  1  1 

lead  to  a  determination  of  a(0)  and  3(0)  .  We  note  that  Z  n?  =  1 

i  1 

implies  immediately  that  a(0)  can  be  determined  in  terms  of  3(0)  , 
that  is,  a( 0)  =  l/  texp(3(0)ei)  ,  after  which  3(0)  can  be  determined 
from  the  second  equation. 

Equations  for  a(t)  and  p(t)  can  be  obtained  from  the  general 
rule  of  §2  involving  orthogonality  conditions.  In  the  present  case 
it  is  clear  that  system  (7^)  has  a  first  integral  of  the  form 

2  0  2  0 
T  +  — te.n.  =  const.  =  T  +  —  te.n;  ,  Z  n.  =  1  . 

3^11  3  i  1  i  i  i 

Whence,  by  virtue  of  the  fact  that  the  ^-function  converges  to  zero 
as  T  *•  •  ,  we  have  that 

T(t)  +  Ite^t)  =  T°4  |Cein°  ,  Z  n^t)  =  1  . 


(75) 


From  the  second  equation  in  (75)  cr(t)  =  l/  t  exp(P't)e  '  ,  after  vhi 

i 

the  first  equation  gives  the  connection  between  T(t)  and  Si t)  •  S 
stituting  T(t)  and  n(t)  ,  expressed  in  terms  of  0(t)  ,  into  the 
second  equation  of  (7^)  we  obtain  a  differential  equation  for  3! t  [ 
We  note  finally  that  higher  approximations  are  also  constructed 


in  [25]  . 
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Chapter  3 

Boundary  Value  Problems  for  Singularly 
Perturbed  Equations  of  Conditionally 
Stable  Type  in  the  Critical  Case 

In  the  previous  chapters  we  assumed  that  the  matrix  F  (cp( t,a)  ,t,0) 
evaluated  along  a  family  of  solutions  of  the  reduced  equation  had  the 
eigenvalue  X  =  0  of  multiplicity  k  arid  that  its  other  eigenvalues 
satisfied  the  inequality  Re  X  < 0  (Condition  III)  .  However,  it  fre¬ 
quently  happens  in  applied  problems  (cf.  §3)  that  this  matrix  also  has 
eigenvalues  satisfying  Re  X  >  0  in  addition  to  those  with  X  =  0  and 
ReX<0  .  Such  cases  are  naturally  called  cases  of  critical  conditional 
stability,  and  we  shall  investigate  below  the  associated  boundary  value 
problems  (as  opposed  to  the  initial  value  problems  of  Chapters  1  and  2). 

We  examine  such  problems  in  this  chapter  as  well  e 3  applications  of 
our  asymptotic  analysis  to  same  conciete  systems.  In  order  to  do  this 
we  will  make  extensive  use  of  the  ideas,  methods  and  results  of  [13,  §11+]  , 
where  we  investigated  boundary  value  problems  in  the  "ordinary"  conditionally 
stable  cases  (that  is,  X  s  0  is  absent) . 

The  systems  of  equations  considered  in  this  chapter  do  not  have 
the  same  general  form  as  those  in  Chapter  2  •,  instead,  we  study  several 
important  special  cases. 

«1  Boundary  Value  Problems  for  Quasiline ar  Systems 

1.  Statement  of  the  Problem.  In  this  subsection  we  consider  the 


system  of  equations 
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|j,-^  =  A(u,t)y  +  u  B( u,t)  , 

u*.z>  (0<t<l)  (1) 

UIt  =  C^U’t)y  +  ^  D(U5t)  » 

where  z  and  y  are  scalar  functions  and  u  is  a  k-dimensional  vector 
function.  In  this  case  the  system  is  quasilinear  because  it  is  linear  with 
respect  to  z  and  y  .  The  choice  of  such  a  system  is  motivated  in  part 
by  its  occurrence  in  the  study  of  applied  problems  from  semiconductor  theory. 
We  prescribe  for  (1)  the  following  boundary  conditions: 

z(0,p.)  =  z°  ,  z(1,m.)  =  z1  ,  u(0,n)  =  u°  .  (2) 

I.  Suppose  that  the  functions  A(u,t)  ,  B(u,t)  ,  C(u,t)  and  D(u,t') 
are  sufficiently  smooth  in  some  domain  G(u,t)  . 

II.  Suppose  that  A(u,t)  >  0  in  G(u,t)  . 

It  is  clear  that  the  reduced  system 

A(u,t)y  =  0  ,  z  =  0  ,  C(u,t)y  =  0 

has  the  family  of  solutions 

z*0,y=0,u  =  a  , 

where  a  is  an  arbitrary  k-dimensional  vector.  The  matrix  ^(evaluated 
at  z*y=u  =  |j,=  0)  ,  fo” 


i 

f 

i 
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is  equal  to  the  block  matrix 


/o  A(a,t) 

(  1  ° 

\o  c(a,t) 


It  is  easy  to  see  that  has  X  =  0  as  an  eigenvalue  of  multi¬ 

plicity  k  as  well  as  two  eigenvalues  of  opposite  signs  in  the  domain 
G(a,t)  ,  namely  X  (a,t)  =  —yA{ci,t)  .  Thus  we  have  indeed  a  critical 
conditionally  stable  case.  This  leads  to  boundary  layers  at  both  ends 
of  the  interval  [0,1]  . 


2.  Construction  of  the  Asymptotic  Expansion  of  the  Solution. 
The  asymptotic  expansion  of  the  solution  of  problem  (1),  (2)  will  be 
constructed  in  the  form 


x(t,n)  =  x(t,(i)  +  tTx(t0,h)  +  Qx(t1}(a) 

(3) 

(tq  =  t/\l  ,  T1  =  (t-l)/u)  , 

where 

x(t,n)  =  xQ(t)  +  ux1(t)+  ...  +  4nxn(t)  +  ...  , 

1*(V0  =  V(T0}  +  UV(T0)  +  ***  +  ^q)  ■*■•••>  (M 

Qx(t,U.)  =  QqXC^)  +  hQlx(t1)  +  ...  +  itn^1x(T1)  +  ... 


TTx(t0,u)  and  Qx(t^,h)  represent  boundary  series  at  the  left  and  the 
right  ends  of  the  interval  [0,1]  respectively. 


Substituting  (3)  into  (1)  and  replacing  the  right-hand  side  F  t 
the  sum  F  +  hF  +  Qf  (as  in  [13,  §1U,  Subsection  5]! ,  we  obtain  a 
sequence  of  equations  for  determining  x^t),  TTix(T0)  and 
(i  =  0,1,...)  . 

For  Xg(t)  we  have  the  reduced  equation 

A(u0,t)yQ  =  0  ,  zQ  =  0  ,  cCaQ,t)yQ  =  0  , 

from  which  we  obtain 


zo  =  0  »  yo  =  0  »  uo  =  ’ 

where  a(t)  is  an  as  yet  arbitrary  k-dimensional  vector  function. 
For  there  is  the  system  of  equations 

z  dTT  v 

-  a(o(o>  ♦  V,o)V  ,  . 

cUT  u 

-5^-  -  +  v»°> V  • 

The  initial  condition  for  TJqx(Tq)  is  obtained  after  substituting  (3) 
into  (2)  and  has  the  form 


ttqz(o)  =  z°  ,  TrQu(o)  =  u°  -  a(o)  . 


As  usual,  we  also  require  that  tTqX(t0)  -  0  as  tq  -*  ®  ,  that  is, 
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The  as  yet  arbitrary  vector  a{0)  appears  both  in  the  equation  ic' 
and  in  the  initial  condition  (7) •  Moreover,  the  initial  value  of 
TT^(t0)  is  as  yet  arbitrary.  We  will  take  advantage  of  this  arbitrari¬ 
ness  in  order  to  guarantee  that  condition  (S)  is  satisfied. 

To  this  end  let  us  first  describe  the  stability  manifold  0^  for 
system  (6)  •,  it  is  analogous  to  the  one  which  figured  in  our  discussions 
in  Chapter  2.  From  (6)  we  have  that 

drTQu  c(a(o)  +  rrQu,o) 


dTToz  A(a(0)  +  n0u,0) 


Let  us  denote  by 


rru  =  u  (a(o) ,  rr  Z) 


the  solution  of  this  system  such  that  ttqU  =  0  for  ^ Qz  =  0  ,  that  is, 
UQ(a(0),0)  =  0  .  By  virtue  of  Conditions  I  and  II  this  solution  exists 
and  is  unique  in  a  certain  neighborhood  of  the  point  tt^z  =  o  .  Substituting 
it  into  the  first  two  equations  of  (6)  we  obtain  the  system  of  equations 

dff  z 

-gp-  =  Ma(0)  +  uQ(a(o)  ,ttoz)  ,o)TToy  , 

dJT  y  (10) 

The  rest  point  TT^y  =  o,  tt^z  =  0  of  this  system  is  a  saddle  (that  is, 
conditionally  stable) ,  since  the  roots  of  the  corresponding  characteristic 
equation  are  clearly  equal  to  — ,/A(a(0)  ,0)  ,  and  by  virtue  of  Condition  II, 

are  real  and  have  opposite  signs.  System  (10)  can  be  integrated  in  an 
elementary  fashion  by  quadratures,  and  for  stability  as  tq  "  we  :^a*-n 
the  equation  of  the  separatrix  of  the  saddle  as 


1.0 


l.l 


JfflM 

HIM 

*  m 

III2-2 

t  m 

11111== 

i_  |4Q 

u  *** 

II  2.0 

m 

1.4 

III 16 

■■■  ■  ■ 
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tt  y  = 

<r 


■  ( 2  r 


V 


e  df 

A(a(o)  +  u0(a(o),?) ,0) 


1/2 


sgfl(nQz) 


(11) 


The  proper  choice  of  the  sign  in  front  of  the  square  root  (in  the  present 
case,  minus)  was  easy  to  make  from  an  analysis  of  the  phase  plane  of  the 
variables  tTqzj  •  After  linearizing  the  right-hand  side  of  (11)  with 

respect  to  TqZ  we  obtain 


V  =  -n0z/v/A(a(°,0) 


The  formulas  (9)  and  (11)  give  an  analytic  representation  of  the  one¬ 
dimensional  manifold  0Q  having  the  property  that  if  the  initial  value 
TTqX(O)  belongs  to  ,  then  TTqx(‘tq'1  belongs  to  nQ  for  tq  >  0  .  For 
such  a  tT0x(t0)  the  inequality 


HTr0x(-r0)tl  <  c  exp(-KT0)  (tq>0) 


(12) 


is  satisfied,  which  implies  the  validity  of  condition  (8)  . 

Thus,  in  order  that  a  solution  of  system  (6)  satisfy  condition  (8) 
it  is  necessary  to  require  that  the  initial  value  TfQx(0)  belong  to  . 

IIIq.  Suppose  that  the  values  of  ttqz  =  z°  belong  to  the  domain  of 
definition  of  the  solution  (9)  . 

Substituting  (7)  into  (9)  we  obtain  the  equation 

u°  -  a(o)  =  u0(a(0),  z°)  ,  (13) 


which  represents  a  system  of  k  scalar  equations  in  the  k  unknown 
components  of  the  vector  a(0)  . 
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IV.  Suppose  that  equation  ( 13)  has  a  solution  a(0)  =  aC  . 

If  we  take  a(0)  =  a®  and  define  the  initial  value  TTQy( 0)  by- 
means  of  equation  (11)  (for  this  we  must  put  TT^z  =  z^  and  a(0)  =  ct° 
in  the  right-hand  side  of  (11)),  then  hQx(0)  belongs  to  Qq  ,  and  con¬ 
sequently,  TTqX(Tq)  satisfies  the  inequality  (12)  and  condition  (8^  . 

We  note  that  for  the  actual  determination  of  TTqx(tq)  it  is 
necessary  to  substitute  (11)  into  the  first  equation  of  (10)  and  to 
solve  the  resulting  scalar  equation  for  tt^Tq)  with  the  initial  con¬ 
dition  h^z(O)  =  z°  .  The  functions  TT^u(Tq)  and  n^(rQ)  are  determined 
by  formulas  (9)  and  (11)  once  is  found. 

Thus  ^^(Tq)  is  completely  determined,  while  for  the  as  yet 
unknown  function  uQ  *  a(t)  we  have  the  initial  value  a°  .  The 
function  a(t)  is  determined  completely  by  the  following  steps. 

The  equation  for  x^(t)  (that  is,  for  z^(t),  y^t)  and  u^(t)) 
has  the  form 


~dF  =  A(u0’t)yl  +  V  Vt)y0Ul  +  B(u0,t}  ’ 


“■“Q  _  _  _ _  _ 

—  *  C(u0,t)yi  +  Cu(u0,t)yoUl  ♦  D(u0,t)  . 

Hence,  by  virtue  of  (5)  we  have  that 


^1  *  0  ’  yl  *  -B(a(t),t)/A(a(t),t) 


(1U) 


da 

dt 


-c(a,t)B(a,t)/A(a,t)  +  D(a,t)  . 


(15) 


Equation  (15)  is  a  differential  equation  for  the  unknown  function  a(t)  . 


V.  Suppose  that  equation  (15)  together  with  the  initial  condition 
a(0)  =  a°  (see  IV)  has  the  solution  a  =  a(t)  for  0  <  t  <  1  . 

Thus  x  (t)  is  completely  determined.  Concerning  x^(t),  the  formula 
(14)  defines  z^t)  and  y^t)  ,  while  u^(t)  is  as  yet  undetermined, 
that  is,  it  is  possible  to  set  u^t)  =  3(t)  ,  where  0(t)  is  an  as  yet 
arbitrary  k-dimensional  vector  function. 

For  tt^x(Tq)  we  have  the  system  of  equations 


=  A(T0)rriy  +  Au(T0)Troy(TTiu+  0(0))  +  cp^)  , 

-jfc-v  > 

dTT  u 

=  c(T0)rriy  +  0(0))  +  cp2(t0)  , 

where  A(tq)  ■  A(a°+  TfQu(T0) ,0)  and  analogous  meanings  are  attached  to 
the  terns  Au(tq)  ,  C(tq)  and  Cu(tq)  ,  while  cp^Tq)  and  co2(tq)  can 
be  expressed  in  terms  of  known  functions  and  satisfy  the  exponential 


estimates  IIcp^Tq)!!  <  cexp(-KTQ)  .  The  supplementary  conditions  for 
^ix( tq)  have  the  form 


^*(0)  =  0  ,  nlU(0)  =  -0(0)  , 

1TjX(»)  =  0 


(17) 


As  in  the  case  of  TTqX(Tq)  we  can  take  advantage  of  the  arbitrariness 
of  0(0)  and  choose  it  so  that  condition  (17)  is  satisfied.  Let  us  make 
the  change  of  variables  in  the  system  (16) 


(18) 


C(t0) 

V  -  51 »  V  -  82  *  V  =  51 +  63 


It  is  easy  to  verify  that  we  obtain  the  system 


dTo  ’  A(,o> 


V<VS1  +  A(T0)S2 


*  *U(VV(V(V9(0|)  +  W 


c 

d8  C(T  )A  (T  ) 

=  lCu(’o>  A(T0)  -)V(tO»63^(0» 

C(r  ) 

+  [c02(TO)  -  aRJ  1  • 

We  note  that  C(t0)/a(tq)  plays  the  role  of  H(t)  in  Chapter  2,  and 
that  the  application  of  this  change  of  variables,  as  in  Chapter  2,  leads 
to  a  system  (19)  in  which  the  equation  for  8^  can  be  separated  from 
those  for  6^  and  6^  • 

The  supplementary  conditions  for  6^  ,  8^  and  6^  are 

B1(0)  «  0  ,  83(0)  =  -3(0)  ,  8i(-)  =  0  (i  =  1,2,3)  • 

The  solution  of  the  third  equation  in  (19)  and  the  initial  condition 
83(0)  =  -3(0)  can  be  written  as 


8  (tq)  *  -3(0)  +  f  Y(T0)Y'1(s)cp  (s)ds  , 


where  Y( rQ)  is  a  fundamental  matrix  of  the  corresponding  homogeneous 
system  (Y(0)  =  E^)  ,  having  the  same  properties  as  the  function  Y(t) 
in  Lemma  U  of  Chapter  2.  The  function 


®3iv  =  VV  ■  A(0  wivno 


W-i 


satisfies  an  exponential  estimate. 

The  condition  5  (®)  =  0  uniquely  determines  0(0)  : 

•  -i 

0(0)  =  Y(-)J*  *' (s)cd  (s)ds  • 

0  * 

By  virtue  of  the  exponential  convergence  of  Y(Tq)  to  Y(®)  as  Tq  ® 
(Lemma  U,  Chapter  2)  we  obtain  for  &3(t0)  the  estimate 

IIB  (T0)||  <  cexp(-HTQ)  (tq>0)  . 


Having  defined  6^(tq)  we  now  write  the  first  two  equations  in 


(19)  as 


d5  A  (t  )C(t  ) 

-  a(t0)  ■  V'To>6i  *  A<V5a  *  *<V 


dB, 

_ _ £ 

dT, 


“  ^  » 


where  ♦(▼q)  is  311  exponentially  decreasing  function.  The  homogeneous 
system  corresponding  to  (20)  is  the  variational  system  of  (10)  .  Hence, 
by  virtue  of  Lemma  U-5  in  [13],  it  follows  that  there  exists  a  unique  solu 
tion  of  (20),  which  satisfies  the  conditions  6^0)  =  0  ,  6^®)  =  0  and 
which  is  exponentially  decreasing,  that  is, 


ll»1(T0)H  -  ce*p(-*T0)  (Tq  >  0  ;  i  =  1,2)  . 

Thus  TLjx(Tq)  is  completely  determined  and  satisfies  an  exponential 
estimate,  while  for  the  as  yet  unknown  function  0(t)  we  have  found  the 
initial,  value  0(0)  .  The  complete  determination  of  0(t)  follows  by 
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steps  analogous  to  those  used  for  the  determination  of  a(t)  ,  except  that 
for  P(t)  we  obtain  a  linear  differential  equation.  Thus  it  is  possible 
to  construct  the  terms  x^(t)  and  TT^x(Tq)  up  to  an  arbitrary  order  n  . 

The  determination  of  the  right  boundary  functions  Q^x(t^)  is 
analogous  to  the  determination  of  the  left  boundary  functions  tt.x(tq)  . 
For  QqX(t^)  we  have 

dQoz  -  dQr/ 

dTx  =  A^U0^  +  ^0U,1^0y  *  dT^  =  SdZ  * 

dQou 

—  =  c(u0(i)  +  Q0u,i)Q0y 

and  the  supplementary  conditions 

QqZ(O)  =  z1  ,  QgX(-®)  =  0  . 

A  fundamental  difference  between  system  (21)  and  the  analogous  system  (6) 
for  TTqX(Tq)  is  that  Uq(1)  is  a  known  quantity,  while  at  the  same  stage 
Uq(0)  =  a(0)  in  system  (6)  is  as  yet  arbitrary.  By  using  this  arbitrari¬ 
ness  to  choose  a(0)  in  a  special  way  (equation  (13)),  we  were  able  to 
satisfy  conditions  (7)  and  (8).  In  system  (21)  there  is  no  such 
k-dimensional  parameter,  but  the  number  of  supplementary  conditions  for 
QqX(t^)  is  clearly  less  than  the  k  identities  obtained  by  comparing 
(7) ,  (8)  since  QqU(O)  is  not  specified. 

Fran  the  first  two  equations  in  (21)  we  have 

dQgU  C(u0(l) +  QqU,1) 
dV  A(u0(l)  +  QgU.l) 

Let  us  denote  by 


(22) 
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Q0u=U1(Q0z)  (23) 

the  solution  of  equation  (22)  satisfying  the  condition  Qqu  =  0  for 
QqZ  =  0  .  Prom  this  we  obtain  (cf.  (H)) 

,V  -  de  V2 

=  (2  J  — — - ^ - )  sgn(Q  z)  .  (24) 

0  A(u0(l)+  1^(5)  ,1) 

The  formulas  (23)  and  (24)  provide  an  analytic  representation  of  the  one¬ 
dimensional  manifold  fi  which  is  analogous  to  that  for  the  manifold  C.Q  • 
It  is  natural  then  to  require  that 

III^  .  The  valueSof  Q^z  =  z^"  belong#  to  the  domain  of  definition 
of  the  solution  (23)  . 

The  initial  values  QqU(O)  and  Q^r(O)  are  determined  by  the 
formulas  (23),  (24)  for  Q^z  =  z1  ,  while  the  solution  qqx(t1)  of  system 
(21)  with  these  initial  conditions  belongs  to  oi  for  <  0  and 
satisfies  the  inequality 
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and  the  supplementary  conditions 

Q1z{0)  =  -  IjU)  =  0  ,  Q1x(-«)  =  0  ,  (27) 

where  A^)  =  A(uq(1)  +  QquCt^jI)  and  analogous  meanings  are  ascribed 

to  A  (t, ),  C(f,  )  and  C  (t.  )  .  The  functions  )  and  f  (t  )  are 

u  1  1  u  1  11  2  d 

known  and  satisfy  an  exponential  estimate  of  the  type  (25) • 

By  means  of  a  change  of  variables  like  (18)  it  is  a  simple  matter  to 
prove  that  Q^x(t^)  exists  and  satisfies  an  exponential  estimate  like  (25)- 
Succeeding  terms  Q^x(t^)  follow  in  an  analogous  manner. 

3.  An  Estimate  of  the  Remainder  Term.  Let  us  introduce  in  the  space 
of  the  variables  (x,t)  a  curve  L  composed  of  the  three  pieces: 

Lx  *  {(x,t):  x  -  xQ(0)  +  nox( tq) ( Tq >  0) ;  t  »  0}  , 

L2  =  {(x,t):  x  =  xQ(t);  0  <  t  <  l)  , 

=  {(x,t):  x  =  x0(l)  +  QQx(T;L)(T1<0)i  t  =  1}  . 

We  denote  by  t  the  projection  of  this  curve  onto  the  space  of  the 
variables  (u,t)  .  It  is  possible  to  take  for  the  domain  G(u,t) 
occurring  in  Condition  I  an  arbitrary  6- tube  of  the  curve  t  .  More 
precisely, 

I.  Suppose  that  the  functions  A(u,t)  and  C(u,t) (B(u,t)  and 
D(u,t))  have  continuous  partial  derivatives  with  respect  to  each  argu¬ 
ment  up  to  order  (w2)((wl))  inclusive  in  some  6-tube  of  the  curve  t  . 
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Having  defined  the  terns  of  the  series  (1)  up  to  order  (rn-1)  inclusive, 
let  us  now  denote  by  X^(t,u)  the  k-th  partial  sum  of  the  series  (3),  that. 


=  t  u  (x.(t)  +  TTix( tq)  +  Qix(r1)) 


Theorem  3 •  Under  Conditions  I-V  there  exist  positive  constants  u 
and  c  such  that  for  0  <  u  <  there  exists  a  unique  solution  x(t,tJ 
of  the  boundary  value  problem  (1),  (2)  lying  in  a  cu-tube  of  the  curve  i, 
and  satisfying  the  inequality 

m.  1 

||x(t,u)  -  Xn(t,u)|l  <  cp,  (0  <  t  <  1)  .  (29) 

Proof.  Let  us  set  C  -  z  -  Z...  ,  ,  T|  =  y  -  Y  ,  and  w  =  u  -  U  .  , 

n+ x  ro-x  n+ x 

where  z,  y  and  u  is  the  unknown  solution  of  the  problem  (1),  (2),  and 
Y^^  and  are  the  partial  sums  determined  by  (28).  By  sub¬ 

stituting  these  into  (l),  (2)  we  obtain  for  £  ,Tj  and  w  the  boundary 
value  problem 


^dt  =  A^U0»t^+  Au(U0,t)Y0(t>ti)w+  G1(TJ,w,t,g.)  , 


n-ft  =  C  +  » 


^  =  C(U0,t)71  +  Cu(U0,t)Y0(t,ji)w  +  G  (Tl,w,t,u)  , 


C(0,u)  ,  and  w(0,n)  are  known  and  of  order  (Mt/**2)  • 
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In  the  equations  of  (30)  we  have  isolated  the  linear  terms  whose  coefficient 
are  evaluated  at  the  zeroth  approximation.  The  functions 


G1(T)Jw,t,n)  =  A(Uml  +  w,t)(Ynfl  +  Tj) 

dU  l 

.  >a(Vl  ♦  ».t)  - 
-  *(00,t)11  -  Au(U0,t)Y0(t,u)v 

and  GgCTIjWjtjp,)  ,  which  is  defined  analogously,  satisfy  the  following 
two  important  properties: 

1.  Gi(0,0,t,u)  =  <9(nW2)  ; 

2.  G^(‘n,w,t,u)  is  a  contraction  operator  with  contraction  coefficient 
of  order  G(p.)  for  T)  and  w  of  order  0(p.)  . 

It  is  necessary  to  transform  the  Gi(11,w,t,M.)  into  a  different  form 
for  the  subsequent  analysis.  We  begin  with  the  identity 


A(U 


n+l 


+  =  A^rwl*^  +  VUi»l>t)w 

+  fA(Vi+  w’t}  • 

s  A<Vi’t}  +  Au(um.i»t)w  +  . 


Here  the  function  q^(w,t,ii)  is  clearly  a  contraction  operator  with  con¬ 
traction  coefficient  of  order  for  w  of  order  G(u)  ;  moreover  , 

<^(0,^)  *  0  .  By  expressing  C(Uml+ w,t)  in  an  analogous  form  (cor¬ 
responding  to  q^(w,t,n)  there  is  a  contraction  operator  which  we  denote 


1C3 


ty  qgCwjtjp.))  and  doing  the  same  for  15(1!^^+  w,t)  and  D(Un+1+w,t1  , 
the  functions  G^(7],w,t,p.)  can  be  reduced  to  the  form 

Gi(Tl,w,t,n)  =  n  ai(t,p,)Tl+  nbi(t,|i)w  +  c±( t ,(j.) T^w 

+  qi(wJt,u)Y0(t,4)  +  Qi(T),w,t,u)  , 


where  a.  ,  b.  and  c.  are  certain  bounded  functions  or  matrices  (Here 
11  i 

and  below  for  ease  of  writing  we  will  denote  a  function  or  matrix  by  one 

and  the  same  symbol  us  since  only  the  boundedness  of  this  quantity  is 

important  to  us.),  and  Q^(T),w,t,y,}  is  a  contraction  operator  with  con- 

2 

traction  coefficient  of  order  )  for  71  and  w  of  order  te(q.)  • 

In  addition,  Q^OjC^tjM.)  =  ft(u  )  . 

Let  us  now  replace  w(t,u)  in  the  system  (30)  by  the  function 
P(t,|i)  ,  where  w=  E+  (C(U0>t)/A(U0,t) )  Q  .  An  elementary  calculation 
shows  that  this  system  of  equations  assumes  the  form 


A  (U  t)C(U  t) 

-  A(U0,t>  *  ‘“V”' 

c(u  ,t) 

+  [coY0(t,u)?4  G1(Tl,F+-^-_  C,t,u)  )  , 


w!t  =  c  +  , 


ia  ^  =  h(t,y,)  p  +  [g(  f  ,71  »t  ,M>)  +  tl(  ^»C,fe»lA)  Yn(  t,y,) 


(31) 


Q(  ?  ,7), C  ,fe»u)  1  » 


where 
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h(t,u)  =  $(u+  exp(-K  t/u)  +  exp(-H(l-t)/|i) )  , 


(32) 


G(e  =  4uJ|  +  ua<  +  P*T\  +  > 

c(u0,t)  c(u0,t) 

q(*»e»*»n)  =  -  ty  ^(^  +  A(uQ,t)  £,t,p,) 

c(uQ,t) 

+  q2^+  A(U0,t)  » 

C(U  ,t)  C(U  ,t) 

Q(?»’n,C»t»M.)  =  -  A(uo/t)  +  A(U0,t)  C  ,t,u) 


c(u0,t) 

A(U0,t) 


The  operator  q(£3C,tj4)  is  a  contraction  with  contraction  coefficient  of 
order  $(4)  for  E  and  £  of  order  <$(4)  satisfying  q(0,0,t,u)  =  0  , 


while  the  operator  Q  is  also  a  contraction  with  contraction  coefficient 
o 

of  order  &(u  )  for  E  ,  1)  and  £  of  order  Q(p.)  satisfying  Q(0,0,0,t,4)  = 

«(u"2)  • 


We  will  consider  the  terms  contained  in  the  square  brackets  of  the 
equations  in  (31)  as  nonhomogeneous  terms  by  passing  from  system  ( 3D  to 
an  equivalent  system  of  integral  equations.  Let  us  denote  by  r(t,s,4) 
the  Green's  matrix  for  the  boundary  value  problem  consisting  of  the  first 
two  equations  in  (31)  together  with  the  boundary  conditions  5(0,4)  = 
g(l,>i)  =  0  .  It  is  possible  to  prove  as  in  [8)  that  the  Green’s  matrix 
exists  and  satisfies  the  estimate 


r(t,s,4)  =  &(exp(-K|t-s|/n))  . 


The  solution  of  the  corresponding  homogeneous  system  and  the  boundary 
conditions  £(0, u)  =  ®(u  ),  £(l,u)  =  ©(u  )  has  the  same  order  of 

smallness  as  the  boundary  values.  In  place  of  the  first  two  equations  in 
(31)  we  have  therefore  the  integral  equation 


C(t,,°Y  .  j-VVu.s.u)  ( 

Tl(t,u)  J  0  \ 


u>Y0(s,y.)  e(s,u)  +  G^ 


SgU,T!,C,t,ii)  j  . 


Let  us  denote  by  H(t,s,p.)  the  fundamental  matrix  of  the  homogeneous 
system  u  df/dt  =  h(t,n)f  .  By  virtue  of  (32)  H(t,s,u)  is  bounded.  The 
initial  condition  for  P(t,p,)  is  clearly  of  the  same  type  as  that  for 
w(t,p),  that  is,  p(0,u)  =  <0(un+2)  •  Therefore  the  last  equation  in  (31"! 
can  be  written  as  the  integral  equation 

t 

?(t,ii)  =  ©(u  )  +  p  J  H(t,s,u)[G(F,Tl,C,s,n)  + 

°  (3^ 

q(f,C,s,M.)Y0(s,u)  +  Q(?,T1,C,s,4)  ]ds  . 


The  operator 


s  h"1  J  H(t,s,g,)Q  ds 
-L  n 


by  virtue  of  the  properties  of  Q  is  a  contraction  with  contraction 
coefficient  of  order  &(ji)  for  C,  Tj  and  Q  of  order  ©(u)  •,  moreover, 
R1(0,0,0,t,pi)  =  ©(uIHl)  •  Since 
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Y0(t,u)  =  S(exp(-Kt/u)  +  exp(-H(l-t)/p.) )  , 


and  therefore 

t  t 

J*  H(t,s,p,)Yn(s,(i)ds  =  J  fi(exp(-Ks/u)  +  exp( -k(  l-s)/p.) )  ds 

W  />»  A 


=  6(u)  j 


the  operator 


H(t,s,ii)q(e,g,s,u)  YQ(s,u)ds 


has  the  same  properties  as  •  Let  us  now  set  R(  >M-)  = 

R1  +  R^  +  ©(m.  )  and  substitute  into  the  expression  for  G  the  values  of 

£  and  7]  from  formula  (33)  .  Then  in  place  of  (3*0  we  obtain  the  equation 


t(t,n)  =  J*  H(t,s,u)(u£1  +  cuS2  +  ^  %S2  +  ^  S-^Jds 


(35) 


+  P(e,T),c,t,n) 


By  taking  account  of  the  estimate  for  the  Green’s  function,  namely 

r(t,s,u)  =  &(exp(-*|t-s|/u)  , 

the  estimate  for  YQ(t,u)  and  the  fact  that 
t  i 

f  f  u"  exp(-K|s-p|u)  [exp(-np/ii)  +  exp(-K(l-p)/kiM«tp  ds 
'0*0 


ft(u)  , 


it  is  easy  to  show  that  the  first  term  in  the  right-hand  side  of  (35)  is  a 
contraction  operator  of  the  same  type  as  the  second  term  R(F,Tl,C,t,u)  . 
Thus,  equation  (35)  can  be  written  as 


*(t,u)  =  T1(e,Ti,c,t,p,)  ,  (36) 

where  the  operator  T  is  a  contraction  with  contraction 

coefficient  of  order  for  €,T1  and  Q  of  order  fc(y)  ■,  moreover, 

T1(0,0,0,t,|j,)  =  P(uml)  . 

Substituting  (36)  into  (33)  we  obtain  the  equations 


C(t,u)  =  si(T1,T),cft,n)  =  T2(r,Ti,c,t,u) 
Tl(t,u)  =  S^T^Tl.C.t.u)  a  T  (f,Tl,C,t,u) 


(37! 


in  which  the  operators  T 2  and  are  similar  to  T^  . 

We  now  apply  to  the  system  (36),  (37)  the  method  of  successive 
app-oximations  as  in  [13 ] .  It  is  possible  to  prove  that  for  sufficiently 
small  u  a  unique  solution  exists  in  a  certain  cy-  tube  of  the  curve 
5  =  T)  =  C  =  0,  and  satisfies  the  estimates  £  =  ^(u,1^1)  ,  T|  =  fc(y,ml)  > 


C  =  f>(M.n+1)  •  Hence,  it  follows  also  that  w  =  f»(urHJ')  . 


n+ 1. 


Thus  z  -  2ml,  y  -  Yjh1  and  u  -  U  are  all  of  order  fl>(nm;L)  , 

•yi  1 

and  since  -  Xr  =  ft(y  )  the  inequality  (29)  is  established. 

This  completes  the  proof  of  the  theorem. 
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§2  Other  Boundary  Value  Problems 

1.  Boundary  Value  Problems  of  a  More  General  Type.  In  §1  we  considered 
a  problem  with  the  boundary  conditions  (2)  .  Using  the  results  for  this 
problem  it  is  possible  to  consider  more  general  boundary  conditions.  The 
corresponding  constructions  are  analogous  to  those  which  were  performed 
in  detail  in  [13,  §13]  for  the  case  Re\^<0  and  in  [23]  for  the  con¬ 
ditionally  stable  case.  Therefore  we  confine  ~uT'<'elves  to  a  brief 
description  of  the  constructive  scheme. 

Suppose  that  the  boundary  conditions  for  system  (1)  are  of  the  form 

R(x(0,u) ,  x(l,u^ ) 

h  R(z(0,u),  y(0,u),  u(0,|i)  ,  z(l,n),  y(l,u),  u(l,u))  =  o  ,  (38) 

in  which  the  dimension  of  the  vector  R  is  equal  to  k+2  ,  the  dimension 
of  x  .  We  consider  as  an  auxiliary  problem  the  boundary  value  problem 
(1),  (2)  with  as  yet  arbitrary  values  of  z°,  z1  and  u°  .  We  propose 
to  select  z^,  z^  and  u^  so  that  the  solution  of  the  problem  (1),  (2) 
satisfies  the  condition  in  (38)  .  This  device  was  used  in  [13,  §13]  . 

Let  us  seek  z  ,  z  and  u*^  in  the  form  of  power  series  in  |i  ; 
for  example, 

0  0  0  2  0 
z  =  zQ  +  u  z1  +  u  zg  +  . . . 

Under  Conditions  I  -  V  we  can  construct  an  asymptotic  expansion  of  the 
solution  and  substitute  it  into  equation  (38)  .  By  further  decomposing 
R(x(0,u) ,  x(l,u) )  into  a  power  series  in  y,  we  obtain  equations  for  the 
terms  in  the  series  for  2°,  z1  and  u°  .  Thus,  in  the  zeroth 


10? 

approximation,  we  have  the  equation  (for  simplicity  of  notation  we  omit 
the  lower  index  0  ,  that  is,  we  write  z°  in  place  of  z^  for  example) 

R(z°,TT0y(0) ,  u°,  z1  ,  QqyCO),  uQ(l)  +  Qqu(  0) )  =  0  .  (39) 

We  note  that  TTQy(0)  as  defined  by  formula  (11)  for  =  z°  is 

a  function  of  ?P  and  a(0)  .  In  turn  or(0)  is  defined  by  equation  (13) 
as  a  function  of  z^  and  vP  .  Thus  ^yfO)  is  a  known  function  of  7? 
and  .  Similarly,  uQ(l)  is  a  known  function  of  z°  and  u°  ,  while 

QqU(O)  and  0^(0)  are  defined  by  the  formulas  (.23)  and  (24)  for  QqZ  =  z~ . 

Hence  it  follows  that  QqU(O)  and  Q^yfO)  are  known  functions  of 

1  0  0  0  0 
z  ,  z  and  u  .  This  dependence  on  z  and  u  results  from  the  fact 

that  uQ(l)  enters  into  equations  (23)  and  (24).  Thus  the  equation  (39) 

is  a  (k+ 2) -dimensional  vector  equation  in  the  (k+2)  unknowns:  z°,  21 
and  k  components  of  the  vector  u°  .  If  the  equation  (39)  has  a  solu¬ 
tion  z^  =  Zq  ,  z1  =  Zq  and  u°  =  u°  and  if  the  corresponding  functional 
determinant  D(R)/D(z^,z^,u^)  is  not  zero  at  the  point  (z^  ,  z^  ,  u^)  , 

then  each  of  the  succeeding  equations  can  be  solved  for  z?  ,  z^  , 

(i  =  1,2,...)  .  Moreover,  for  sufficiently  small  values  of  p,  there 
exists  in  a  certain  6-tube  of  the  point  (z^,  Zq,u°)  a  unique  point 

(z^(n),  z1(u)  ,  u^(u))  such  that  the  solution  of  the  equation  (1)  and  the 
boundary  condition 

s(0,n)  *  z°(u)  ,  z(l,u)  =  zX(m.)  ,  u(0,n)  =  u°(u)  (40) 


no 


satisfies  the  boundary  condition  ( 38) .  The  formally  constructed  series 

0  0  11  0  0 
zQ  +  nz1  +  . . . ,  zQ+  nzx  +  •  • . ,  uQ  +  +  ... 

are  asymptotic  series  for  z^(p.)  ,  z^(y.)  and  u*^( p,)  ,  while  the  asymptotic 
expansion  of  the  solution  of  the  problem  (1),  (40)  serves  as  an  asymptotic 
expansion  for  the  basic  problem  (1),  (38).  The  proofs  of  these  assertions 
can  be  given  without  difficulty  by  using  the  methods  in  [13,  §13]  . 

2.  A  Class  of  Boundary  Value  Problems  Reducible  to  a  Type  Already 
Considered.  Suppose  that  a  singularly  perturbed  system  has  the  form 

^It  =  F(Ujt)  1  It  *  c(u»t)y  +  D(u,t)  ,  (41) 

where  y  is  a  scalar  and  u  a  vector,  and  suppose  that  a  certain  boundary 
value  problem  is  posed  for  the  system  (1*1)  .  For  definiteness  we  win  con¬ 
sider  the  fonowing  boundary  conditions 

u(0,vt)  =  u°  ,  y(l,u)  =  y1  •  (1*2) 

(It  is  of  course  possible  to  consider  other  types.)  The  peculiar  thing 
about  the  system  (1*1)  is  that  the  function  F  does  not  depend  on  y  , 
and  therefore  the  usual  algorithm  for  the  construction  of  the  asymptotic 
solution  of  a  singularly  perturbed  problem  is  inapplicable  here.  This 
fonows  because  for  4  =  0  the  equation  F(u,t)  =  0  cannot  be  solved  for 
y  .  One  way  of  circumventing  this  difficulty  is  the  following.  We 
differentiate  the  first  equation  in  (41)  and  use  the  second  equation  to 


obtain 


Ill 


U2  =  F  (u,t)  [C(u,t)y  +  D(u,t)  ]  +  F  (u,t) 
dt*  U  C 

=  A(u,t)y  4  B(u,t) 

If  we  now  introduce  the  new  variable  z  =  >  then  we  are  led  to  the 

Uw 

system 

^dt  =  A(u,t)y  4  B(u,t),  =  z  , 

=  C(u,t)y  4  D(u,t)  .  (13 

It  is  necessary  to  prescribe  for  the  system  (U3)  another  condition  besides 
the  boundary  conditions  (42) ,  since  as  a  result  of  differentiation  the 
order  of  the  system  has  increased  by  one.  This  condition  is  obtained  from 
the  first  equation  in  (41)  by  setting  t  =  0  : 

Thus  the  boundary  conditions  for  system  (43)  have  a  singularity  as 
li  0.  However,  it  is  possible  to  remove  this  singularity  (cf.  Chapter  2, 
§4,  Subsection  1)  by  introducing  the  new  variables  z  =  yz  ,  y  =  uy  •  For 
these  new  variables  we  have 

Wdt  =  A(u»t)5?+  4B(u,t),  y,-^  =  z  ,  ()<4) 

^  -  C(u,t)y  4  yD(u,t)  , 

which  coincides  with  (1)  except  for  notation.  The  boundary  conditions  for 
the  new  variables  are  now  regular  in  y  ,  and  so  we  can  use  the  method 
described  in  Subsection  1  for  the  construction  of  the  asymptotic  expansion 
of  the  solution. 


We  note  that  the  passage  from  system  (1*1)  to  system  (1*1*)  is  not  the 
only  one  which  allows  the  construction  of  the  asymptotic  expansion.  It 
is  possible  to  apply  to  the  original  problem  (1*1)  ,  (1+2)  a  certain  modi¬ 
fied  algorithm  for  the  construction  of  the  expansion  in  the  form  of  a 
regular  part  and  a  boundary  part.  For  simplicity  let  us  consider  the  case 
when  both  y  and  u  are  scalar  functions.  We  seek  a  solution  of  the 
problem  (1*1) ,  (1*2)  in  the  following  form 

y(t,u)  =  y0(t)  +  M-yx(t)  +  ...  +  mT^^Tq)  +  ^Qy(  rQ)  +  ...  + 

Qoy(Ti)  +  +  ...  , 

u(t,n)  =  uQ(t)  +  uu^t)  +  ...  +  ttqu(  tq)  +  uttiU(t0)  +  ...  + 

Q0u(Ti)  +  +  ...  . 

Then  for  yQ(t)  and  uQ(t)  we  ot'tain  system 
—  du^  _  _ 

0  =  F(uQ,t) ,  —  =  C(u0,t)y0  +  D(uQjt)  .  (1*5) 

Suppose  that  cp(t)  =  u^(t)  is  a  certain  root  of  the  first  equation  in 
(45).  Substituting  it  into  the  second  equation  gives  yQ  : 

*0 

Thus  the  leading  term  of  the  regular  part  of  the  asymptotic  expansion  is 
determined  in  a  rather  unusual  way.  For  TT  ^y(  tq)  and  t*^u(Tq)  we 
obtain  the  system 
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dTT.y 

-  F(»(0)  +  ttqu,0) 


d"  u 

=  c^(°) +  v,0)TT-iy 


together  with  the  supplementary  conditions 


Tf0u(0)  =  u°  -  cp( 0)  ,  rr_iy(-)  =  0  ,  ttou(»)  =  o  . 


( 16) 


(17) 


If  F  (cp(0),0)  C(cd(0),0)  >  0  (In  order  to  obtain  an  asymptotic  expansion 
of  the  desired  form  it  is  natural  to  require  that  F  (co(t)  ,t)C(tp( t)  ,t)  >  0  •, 
such  a  condition  appears  in  the  first  approach  considered  for  the  system 
(ll)),  then  the  rest  point  n  ^y  =  u  =  0  of  system  (16)  will  be  con¬ 
ditionally  stable.  The  stability  separatrix  for  ~  oc  Is  described  by 
the  equation 

TT  u  .  / 

ff-l SSorl  ll 0)  dF)  sgnC(<p(0),0)  sgnTT0u  . 

0 

By  substituting  rrQu(  0)  =  u°  -  cp(  0)  into  the  right-hand  side  we 

obtain  the  initial  value  ft  ^y(O)  .  The  solution  of  system  (16)  with 

these  initial  values  satisfies  each  of  the  conditions  in  (17)  and  an 
exponential  estimate. 

The  function  Q^u(t^)  is  found  to  be  identically  zero  (T’his  is 
quite  natural  since  the  function  u  is  not  given  at  the  point  t  =  1  .), 

while  for  and  Q^u(t^)  we  obtain  the  linear  system 

dQ  u 

along  with  the  supplementary  conditions 


11U 


Since  F^(cp(l),l)  C(®(1),1)  >  0  the  characteristic  equation  for  system  («8' 
has  roots  X,  „  =  — y?  C  of  opposite  signs,  and  the  solution  satisfying  the 
conditions  in  (U9)  has  the  form 


=  (y1  -  yQ(l))exp(XT1)  , 

(y1-y0(D^  _ 

v(Ti^ =  F^imr  exp(XTi^  for  x  =  A? 


One  can  now  construct  the  succeeding  terms  of  the  asymptotic  expansion 
in  a  similar  manner. 


3.  Boundary  Value  Problems  for  a  Weakly  Nonlinear  Equation.  In 
Chapter  1  we  considered  the  initial  value  problem  for  the  weakly  nonlinear 
equation 

=  A(t)x  +  pX(x,t,n)  (0  <  t  <  T) 

under  the  assumption  that  the  matrix  A(t)  had  the  eigenvalue  X(t)  =  0 
of  multiplicity  k  and  that  the  remaining  eigenvalues  X^(t)  satisfied 
ReX^(t)  <0  .  If  now  in  addition  to  the  zero  eigenvalue  of  multiplicity 
k  the  matrix  A(t)  has  m^  eigenvalues  X^(t)  such  that  ReX^(t)<0 
and  m^  eigenvalues  X^t)  such  that  ReX^(t)>0  (with  k+m^+m2  =  m)  , 
then  ve  would  like  to  obtain  the  same  qualitative  results  as  in  Chapter  1. 
In  other  words,  we  want  to  determine  which  solution  of  the  reduced 
equation  A(t)x  =  0  is  the  limit  as  ^  -  0  of  the  solution  x(t,p^  • 

To  achieve  this  it  is  necessary  to  consider  in  place  of  an  initial  value 
problem  a  boundary  value  problem  in  which  at  least  components  of 
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x  are  prescribed  at  t  =  0  and  at  least  m2  components  are  prescribed  at 
t  =  T  .  The  asymptotic  expansion  of  the  solution  will  have  the  form  (3), 
but  the  terms  of  the  expansion  are  determined  using  certain  modifications 
in  the  construction  procedure  which  are  analogous  to  those  employed  in 
Chapter  1.  The  details  are  given  in  [16],  where  a  similar  boundary  value 
problem  for  a  weakly  nonlinear  system  of  difference  equations  is  considered. 


53  Applications 

1.  A  Problem  from  the  Theory  of  Transistors.  We  first  make  some 
explanatory  remarks  of  a  physical  nature.  Consider  a  contact  (w  =  01 
between  two  semiconductors  of  different  types,  leading  to  a  one-dimens ional 
problem.  To  the  left  of  the  contact  (-t  <  w  <  0)  we  place  a  semiconductor 
of  p-type,  while  to  the  right  (0  <w  <  t)  a  semiconductor  of  n-type, 
that  is,  a  (p-n)  junction.  Such  a  semiconductor  scheme  can  be  described 
by  a  system  of  equations,  consisting  of  Poisson's  equation 


(Here  E  is  the  polar  electric  voltage,  p,  n,  N*  and  are  the 

respective  concentrations  of  holes,  electrons,  donors  and  acceptors,  q 
is  the  electron  charge,  and  e  is  the  dielectric  permeability. )  and 
the  equations  for  the  holes  (ip)  and  *-he  electron  current  (in) 

ip  -  WpEP  -  «»,  , 

ln  ■  4  • 


(50) 


(5D 
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(Here  u,  ,  u  are  the  mobilities  of  holes  and  electrons,  respectively,  and 
p  n 

D  ,  D  are  the  diffusion  coefficients  of  holes  and  electrons.) 

p»  r 

It  is  known  that  =  M- T/^>n  =  <l/^T  »  where  k  is  Boltzmann's 

constant  and  T  is  the  temperature.  We  will  assume  that  the  problem  is 
stationary  in  time  and  that  there  are  no  externally  generated  sources.  Then 


from  the  continuity  equation  it  follows  that  i  and  i^  are  constant. 

We  will  consider  the  following  special  case,  namely,  to  the  left  of  the 
contact  =  N,N*  =  0  ,  while  to  the  right  =  0,  N*  =  N  .  Let  us  now 


introduce  the  dimensionless  variables 


t  =  w/t,  y  =  Etq/kT,  v1  =  p/N,  v 2  =  n/N  , 

V'Vp  -  <=!•  yA,Bn  -  v  “2  *  ,m/’2“2  • 

Then  the  system  (50)  -  (51)  can  be  written  in  dimensionless  form  as 
=  V1  •  V2+  N(t)  ’ 


dt  =  Viy  "  cl*  dt  "  "V  T  w2  * 


=  -vny  +  c. 


-1,  -1  <  t  <  0 

,  and  u  is  a  small  quantity  of  order 

1,  0  <  t  <  1 

It  is  possible  to  consider  various  kinds  of  boundary  conditions  for 
the  system  (52) .  We  restrict  ourselves  here  to  one  of  the  simplest, 
known  as  the  symmetric  case,  in  which  c^  =  c^  *  c  is  a  given  constant, 


and  so  the  boundary  conditions  are  given  separately  for  the  intervals 
[-1,0]  and  [0,1]  •,  namely, 
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v1(-l)  =  1  ,  v2(-l)  =  0  ,  v1(0)  =  v2(0)  ,  (53) 

and 

Vi«»  =  v2(0)  ,  v^l)  =  0  ,  v2(l)  =  1  .  (5M 


In  the  present  case  the  problem  (52),  (53)  reduces  to  the  problem  (f-2)  ,  (5I 
under  the  changes  of  variables  t  -  -t  ,  v^  -*  v^  and  v2  "*  •  Therefore 
it  suffices  to  consider  only  one  of  these  problems.  We  will  consider  the 

system 


v2+  ! 


-v2y  + 


( o  <  t  <  1) 
C  > 


(55) 


along  with  the  conditions  in  ( 54)  . 

By  introducing  the  new  variables  up  =  vi  +  v2  »  u2  =  vi  ~  v2  we 
obtain  the  system 


U2+  1 


9 


du^ 

dt 


^2 

dt 


V 


-  2c 


(56) 


which  is  clearly  of  the  same  type  as  (41)  since  the  right-hand  side  of  the 
first  equation  does  not  contain  y  .  Proceeding  as  in  Subsection  2  of 
j2  ,  that  is,  differentiating  the  first  equation  and  introducing  the  new 
variables 


y s  » z  = 


9 


A 
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we  obtain  the  system 


dz  dy 

^dt  =  V  ■  »  “dt  =  2  » 

dux  du2 

V  >  V  •  ^  • 


(57) 


This  system  is  of  the  type  studied  in  §1  ,  where  the  dimension  of  the 
vector  u  is  now  two.  The  boundary  conditions  (fjl)  in  the  new  variables 
have  the  form 


Ug(0)  =  0  ,  u^l)  =  1  ,  u2(l)  =  -1  .  (58) 

In  addition,  it  is  necessary  to  supply  a  further  condition,  which  is 
obtained  from  the  first  equation  (5 6)  in  the  two  forms 


z(0)  =  1  , 

(59) 

z(l)  =  0  . 

(60) 

It  is  easy  to  see  that  of  the  five  conditions  (58) ,  (59)  and  (60)  we 
need  only  consider  the  four 

z(0)  =  1  ,  z(l)  =  0  ,  u^l)  =  1  ,  u2(l)  =  -1  ,  (61) 

since  the  condition  Ug(0)  =  0  is  automatically  satisfied.  Indeed,  from 
the  first  and  last  equations  in  (57)  it  follows  that  dz/dt  =  du^/dt  ; 
whence,  noting  the  boundary  conditions  z(l)  =  0  ,  Ug(l)  =  -1  ,  we  obtain 
that  z  =  Ug  +  1  .  Hence,  by  virtue  of  the  condition  z(0)  =  1  it 
follows  that  Ug(0)  =  0  . 
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Thus  we  have  a  boundary  value  problem  (57) ,  (6l)  which  is  of  the  same 
type  as  the  problem  (1),  (2)  in  §1  ,  the  only  exception  being  that  the 
boundary  condition  for  u  is  given  at  t  =  1  rather  than  at  t  =  0  . 

Let  us  now  construct  the  asymptotic  expansion  of  the  solution  as  in 
§1  •  We  have  first  that 


Z0  =  °  ’  y0  =  °  ’  U1  =  aL^t)  ’  U2  =  a2(t)  ’ 


where 


and  a2  are  as  yet  arbitrary  functions.  Since  u  is  prescribe: 


at  t  =  1  we  consider  first  the  system  of  equations  for  Q.x(t  ' 

0  1 

(t^  =  (t-l)/|i)  (Note  that  the  roles  of  Nx  and  Qx  in  the  present 
problem  are  interchanged  relative  to  §1  . ) 

dQ0Z  dV 

dTx  =  +  Q0Ul)Q0y  ’  "dT^  =  Q0Z  » 


dQ0U~l 

d*r. 


-1  dQ  u 

=  (a2(l)  +  Q0u2)Q0y  ,  dT^  =  (0^(1)  +  Qq^)^  * 

The  supplementary  conditions  for  Q^x  are 

Q0z(0)  =  0  ,  QqU^O)  =  1  -  OjU)  ,  QqU^O)  =  -1  ’  a2(1)  1 

QqX(-®)  =  0  . 


(62) 


(63) 


Let  us  construct  the  manifold  for  the  system  (62). 

from  (62)  the  equations 

dSpUj  ag(i)-.yg  av2 


We  obtain 


dV  '  °id)  +  Q0\  ’  dQ0z 


=  1 


120 


The  solution  of  this  system  and  the  initial  condition  Qqu  =  0  for  Q^z  =  0 

is 

2  2  V2 

=  -0^(1)  +  (sgna1(l))  (a£(l)+  2a2(l)QQz+  qqz)  f  (61) 

Q0u2  =  V  • 

Equation  (61)  is  analogous  to  equation  (9)  of  §1  ,  while  the  equation 

V  -  1/2 

(ly  =  (sgna  (1))  (2  J  — ~ ' — ‘ - p-y 7^)  sgn(Q  z)  (65) 

V  1  o  (^(D+2a2(iK+  k2)1^  0 

is  analogous  to  equation  (11)  .  The  equations  (61)  and  (65)  provide  an 

analytic  representation  of  the  manifold  0^  . 

By  substituting  the  boundary  values  (63)  into  (61)  we  obtain  equations 

for  0^(1)  and  a2(l)  ,  that  is,  1  -  0.(1)  =0,-1-  ag(l)  =  0  .  Hence 

OjU)  =  1  ,  a2(l)  =  -1  ,  (66) 

and  the  function  Qqx(t^)  is  easily  seen  to  be  identically  zero  in  the 
present  case. 

The  system  of  equations  for  a^(t)  and  a0(t)  is  obtained  in  the  usual 
way  and  has  the  form  (cf.  (15)  in  §1) 
dCL  dQL 

-ST  -  2c  V°i  ’  “5T  =  0  • 

The  solution  of  this  system  arid  the  supplementary  conditions  (66)  is 
OjU)  =  (l+lcd-t))1/2  ,  a2(t)  =  .1  . 
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We  have  for  TTqx(t0)  the  system 


dh  z  dTT  v 

-57-  + 


dTT  u.  dTT  u 

^  -  I'l^oVV  >  -  (0i(0)  *  Vi'V 

1/2 

(for  a,  (0)  =  (1+  l*c)  '  )  and  the  supplementary  conditions 


(67' 


TT02(0)  =  1  ,  ttox(«)  =  0 


The  manifold  QQ  for  (67)  is  represented  by 

l/2  2 

TT^  =  -(1+Uc)  '  +  (l+4c  -  2TT0z+  TT0z)  , 


V2-V  > 


TT„Z 


V  -  -(2  1° 


1/2 


(l+tc-2!|+?  ) 


2-372)  sgnC^z) 


(68) 


It  is  necessary  now  to  substitute  the  last  equation  of  (68)  into  the 
first  equation  of  (67)  and  to  solve  the  resulting  differential  equation 
for  TTqZ(t0)  together  with  the  initial  condition  TTQz(0)  =  1  .  After  this 
the  rest  of  the  function  TTqx(t0)  is  determined  from  (68)  . 

The  construction  of  the  succeeding  terms  of  the  asymptotic  expansion 
can  be  executed  as  in  §1. 

A  comparison  with  experiments  shows  that  that  the  application  of  the 
asymptotic  method  under  consideration  is  suitable  already  in  the  zeroth 
approximation  with  a  high  degree  of  accuracy  in  processes  involving  tran¬ 
sistors.  A  more  detailed  physical  analysis  of  the  mathematical  results 
is  given  in  [14]  . 
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2.  Some  Control  Problems.  Certain  problems  in  control  theory  appear 
as  singularly  perturbed  boundary  value  problems  of  the  critical  condition¬ 
ally  stable  type.  As  an  example,  consider  the  linear  problem  of  Mayer - 
Bolza  [22] 

^It  =  All^Z  +  Ai2^y  +  ]Vt)u  » 

=  A^tjz  +  A22^^y  +  Vt)u  » 
z(0,4)  =  z°  ,  y(0,n)  =  y°  , 

min{d  x(l,^)  +  -|j  [x  (t,p,)F(t)x{t,n)  +  u  (t,n)R(t) u(t,n) ]dt]  . 

Here  z  and  y  are  M-  and  m-dimensional  phase  vectors,  respectively, 

* 

u  is  the  control,  x  denotes  z  and  y  taken  together,  denotes  transpose, 
and  F(t)  ,  R ( t )  are  symmetric  matrices. 

Suppose  that  no  supplementary  conditions  involving  bounds  on  the 
control  are  imposed.  Then  the  problem  becomes  a  classical  problem  in 
the  calculus  of  variations.  By  applying  the  method  of  Lagrange  multipliers, 
we  can  reduce  the  problem  to  the  following  one  for  the  auxiliary  Lagrangean 
vector  functions  t,jj.)  and  Xg(t,M,)  of  dimension  M  and  m  ,  respec¬ 
tively,  namely 

dt  *  A11Z  +  B1R  ^®i*i  +  Ai2y  +  B1B  ^2*2  ’ 

*1  * 

**  dt  *  Fllz  "  AllXl  +  F12y  '  A21X2 


9 
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4?  =  A„,z  +  B0R  ^X.,  +  A22y  +  B2R  j 


dt  21 


11 


dX, 


dt  =  F212  "  A12X1  +  F22y  "  A22X2  ’ 

z(0,^)  =  z°  ,  X1(l,n)  =  -d^/u  ,  y(0,u)  =  y°  ,  Xg(l,n)  =  -dg 


Here  we  have  denoted  by  F. .  and  d.  the  appropriate  blocks  of  the  matrix 

lj  1 

F  and  the  vector  d  • 

Thus  we  obtain  a  problem  with  a  singular  boundary  condition  at  t  =  1. 
By  changing  variables  as  in  Subsection  2  of  §2  we  are  led  to  a  problem 
having  nonsingular  boundary  conditions  and  a  corresponding  matrix  with 
zero  eigenvalues.  In  order  to  apply  our  asymptotic  methods,  it  is  neces¬ 
sary  to  assume  that  the  matrix 


has  M  eigenvalues  with  negative  real  parts.  Then  we  obtain  a  con¬ 
ditionally  stable  system  in  the  critical  case. 

Similar  kinds  of  systems  occur  in  other,  more  complicated  problems 
of  optimal  control  (cf.,  for  example,  [1])  . 


§4  The  Case  of  an  Incomplete  Set  of  Eigenvectors 


In  Chapter  1  we  showed  for  the  system  of  two  linear  equations  (38) 
that  if  the  number  of  eigenvectors  corresponding  to  X  =  0  is  less  than 
the  multiplicity  of  this  eigenvalue,  then  in  order  to  obtain  a  solution 
bounded  as  p,  “•  0  we  must  pose,  in  general,  a  boundary  value  problem. 
Moreover,  the  asymptotic  expansion  of  the  solution  will  contain  fractional 
powers  of  p. 

We  now  consider  this  question  for  a  certain  nonlinear  system. 

1.  A  System  of  Two  Nonlinear  Equations.  Let  us  consider  the  system 

V-  -  ^(z.y)  +  vif^z.y.t), 

d  0  <  t  <  1  (69) 

^  dt  ”  F2(z,y)  +  ^2  <*»*•*>» 

where  the  functions  F^(z,y)  and  fi(z,y)  (i=l,2)  are  sufficiently  smooth 
for  z  in  (z^z^,  y  in  (y1>y2)  and  t  in  [0,1). 

Suppose  that  the  equation  F^(z,y)  =  0  has  a  root  y  =  cp(z)  in  (y^,y2) 
for  z  in  (z^,z2)  such  that  F^y(z,cp(z))  ^  0,  and  also  that  F2(z,tp(z))  =  0 
for  x  in  (z^,z2>.  Then  the  reduced  system  corresponding  to  (69)  has 
the  family  of  solutions 

G>C  q>( ay)  *  “  in  <Z1»Z2^* 

and 

det  F  ■ 
x 


f12(z.<P<*»  Fly(z,qi(z)) 
T2t(zM*))  F2y(z,«p(z)) 


(71) 


* 
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for  z  in  (z^.z^).  In  addition,  we  assume  that 


plz(z.9(0)  +  F2y(z,<p(z))  =  0  in  (zltz2).  (72) 

It  follows  from  (71)  and  (72)  that  the  matrix  has  a  zero  eigenvalue 

of  multiplicity  two  (X^  =  =  0)  to  which  corresponds  the  single  eigcn- 

vector  V(C)K 

Let  us  introduce  new  variables  z  and  w  =  F^(z,y)  in  (69),  and 
note  that  in  a  neighborhood  of  y  =  cp(z)  we  can  define  a  function 
y  *  Y(z,w)  for  which  Y(z,0)  ®  q>(z).  Then  the  system  (69)  becomes 


^  *  w  +  p,f(z>w>t), 

^  dt  m  A<2*w>w  +  M.g(z,w,t), 


(73) 


where  f(z,w,t)  *  f1(z,Y(z  ,w)  ,t) ,  A(z,w)  »  (Fu  +  ply  (z,Y(z,w)) ,  and 
g(z,w,t)  -  (Flzfl  +  Fly  f2)(z,Y(z,w),t).  We  have  that 


and  so 


lim 

w-*0 


F2(z,Y) 

F^z.Y) 


=  lim 
vO 


F2v<z»Y> 

Fly(z  *Y) 


p2v  (*,?(*>) 
Fiy(z.«P(z>)  * 


A(z,0)  -  Flz(z,q>(z))  +  F2y(z,tp(z))  -  0 

by  virtue  of  (72).  Thus,  the  order  of  A(z,w)  as  w-»0  is  at  least  w 
and  we  can  write  (73)  in  the  form 

||  ■  v  +  nf  (z,w,t), 

**  ft  "  B<*>w>w2  +  4g(*»w,t), 


v  -  0 

g(*,0,t)  ■  0. 

We  assume  Chat  Che  second  of  these  equations  has  a  root  z  ■>  z(t) 
in  (z^.Zj)  and  that  gz(z,0,t)  >  0  for  t  in  [0,1].  The  characteristic 
equation  corresponding  to 

-A  1 

-  0 

gr(z,0,t)  -A 


defines  a  pair  of  characteristic  values  with  opposite  signs,  namely 
\2  “  ±  <gz(*,0,t))1/2. 

Thus  the  system  (74)  is  of  conditionally  stable  type  with  small  parameter 
to  which  we  can  apply  the  theory  developed  in  [13,  Sec.  14]. 

We  consider  now  various  supplementary  conditions  for  the  system 


Let  us  first  prescribe  the  initial  conditions 


Then,  in  general,  the  solution  of  this  initial  value  problem  for  the 
conditionally  stable  system  (74)  is  unbounded  as  n  -•  0.  (This 
phenomenon  occurs  even  if  the  term  1  / /\i.  is  absent  from  the  expression 
for  v(0,*i).) 

2°.  Let  us  now  prescribe  the  boundary  conditions 

2(0,H)  -  2°,  y(l,n)  -  yl,  (76) 

which  we  write  as 

z(0, n)  -  2°,  &  v(l,n)  -  Fl(2(l,n),y1), 

in  the  variables  z  and  v.  This  boundary  value  problem,  in  principle, 
admits  a  solution  bounded  as  u  “•  0,  and  its  asymptotic  expansion,  which 
features  right  and  left  boundary  layers,  consists  of  powers  of  Sii.  The 
question  of  the  existence  of  such  a  solution  is  investigated  by  the 
method  in  113,  Sec.  13]  (cf.  also  [23]  and  §2  of  this  chapter). 

These  remarks  also  apply  to  the  more  general  boundary  conditions 
R(2(0,n),  z(l,l0,  y(0,|i),  y(l,n))  -  0. 

As  *o  Illustration  of  the  theory  consider  the  system 


^Zt"v  +  ‘^z’  at  “  *<z+y)2  ■  (z+y)  + 

which  can  be  written  as  (cf.  (74)) 

v^dt  =  V  +  v^2’  ^  ft  =  +  2  +  1  + 

The  solution  of  this  system  satisfying  the  boundary  conditions  (76)  has 
the  following  asymptotic  representation  valid  to  order  G(v^i): 

*  «  z  +  nQz  +  QqZ,  y  -  y  +  nQy  +  0Qy. 

Here  z  **  -y  ■  -(1+t)  and  II^z  *  -Il^y ,  QqZ  =  -Qgy  are  found  by 
quadratures  from 

n  z  ■  -sgn  n  z(|  exp{-2n0z]  +  nQz  -  |)1/2, 

0 

nQz(0)  *  z°  +  1,  T0  -  t//n; 

~~  QjjZ  -  sgn  QQz(i  expt-20^]  +  0Qz  -  \)lf2 , 

00z(0)  <*  2  -  yl,  rx  -  (1  -  t)Mi. 

2.  A  Problem  Arising  in  the  Theory  of  Singular  Optimal  Control.  [This 
problem  has  been  investigated  by  M.  G.  Dimitriev.)  Suppose  that  it  is 
required  to  minimize  the  functional 

J  -  (77) 


•long  trajectories  of  the  system  of  two  equations 


f  ^ (X1 »x2  1  ^  ^ * 


0  <  t  <  1 


f2^Xl,X2,t^  +  ^3^x2 ,t^u* 


Xj^O)  =  x°,  x2(0)  =  x°. 


There  is  no  bound  on  the  control  u  of  the  type  involving  a  closure 
inequality.  Now  this  problem  need  not  have  a  solution  in  the  class  of 
continuous  functions  u  and  so  there  is  the  question  of  the  construction 
of  a  generalized  solution  of  the  problem  (77),  (78). 

To  this  end  we  introduce  a  regularized  problem,  that  is,  in  place  of 
the  functional  (77)  we  consider  the  functional 


J  »  cp(x1(l,^))  +  £ 


J  U2dt. 


If  we  then  introduce  conjugate  variables  (Lagrange  multipliers),  we  can 
reduce  this  problem,  as  in  Section  3  of  the  present  chapter,  to  a 
system  of  differential  equations  with  boundary  conditions,  namely 


dt  * 


dx2  2 
V  ft-  *  ^2(xl,x2,t)  +  f3(x2,t)'*,2’ 


dt  "  "^1x^1  +  f2x1*2)’ 


4 


dt 

dt 


*1 


-  -  -H(f,  4  +  f,  f_)  -  f. 

1x2t1  2x2  2  3x2 

(0,4)  *  x°,  x2(0,4)  =  x°, 

(1,4)  =  -®  (x  (1,4)),  *,(1,4) 
1 


=  0. 


(81) 


The  right-hand  sides  of  the  equations  for  the  fast  variables 
and  ^2  have  the  same  properties  as  in  the  system  (69).  (The  presence 
of  the  slow  variables  x^  and  \Jr ^  offers  no  complication;  cf.  Section  5 
of  Chapter  2.)  In  fact,  if  we  assume  that  f^CXj.t)  4  0,  then  setting 
4  *  0  in  the  second  and  the  fourth  equations  of  (80)  we  obtain  *2  =  0, 
while  x2  remains  undetermined.  The  determinant  (71)  here  has  the 
form 


0 

0 


and  consequently,  ^  **  0;  thus  condition  (72)  holds. 

The  change  of  variable  described  in  Subsection  1,  having  the  form 
2 

v  ■  leads  to  a  problem  in  which  the  characteristic  values  are 

again  equal  to  zero.  However,  by  using  the  change  of  variable 
v  ■  we  obtain  in  place  of  (80),  (81)  the  following  problem 


XjCO.lO  *  x°,  x2(0,ti)  =  x°, 

♦j^l.lO  *  -9X  (XjCl.li)),  v(l,n)  *  0. 

This  problem  is  conditionally  stable  if  we  assume  that 

<^20^10  <  °* 

An  investigation  of  this  problem  by  the  methods  discussed  in 
Subsection  1  reveals  that  the  leading  term  of  the  asymptotic  expansion 
for  x2(t,J*)  has  a  boundary  layer  at  the  left  endpoint.  The  asymptotic 
expansion  of  the  control  u(t,lO  to  order  &(>40  is  found  from  the  formula 
u  *  v/^fi,  once  we  have  found  the  expansion  of  v  to  order  After 

the  corresponding  calculations  have  been  made,  we  find  that  to  order 
®(<40  the  control  u(t,l*)  has  the  form 

V<V 

u(t,n)  -  +  vL(t)  +  nlV(T0)  +  q  v(T  >,  (85) 

while  the  optimal  trajectory  is 


f‘V2 


(83) 


(84) 


xx(t,n)  =  x1Q(t),  x2(t,^)  -  *20(t)  +  n0X2(T0)- 


(86) 


From  the  relation  (85)  it  follows  that  the  leading  term  of  the  asymptot 
expansion  as  ^  -•  C  has  the  character  of  a  6-function. 

Remarks.  1.  Using  singular  perturbation  theory,  we  can  define  a 
class  of  singular  functions  in  which  the  problem  (77),  (78)  is 
solvable . 

2.  The  method  discussed  above  can  be  applied  in  the  vector  case. 
For  example,  if  u  and  x^  are  k-dimensional  vectors  and  f^  a 
(k  X  k) -matrix,  then  the  problem  of  zero  characteristic  values  can  be 
eliminated  by  means  of  the  change  of  variable  v  * 


Chapter  t 


Singularly  Perturbed 
Integro-differential  Equations 
in  the  Critical  Case 

U  Statement  of  the  Problem  and  Auxiliary  Results 
1.  Statement  of  the  Problem.  The  results  of  Chapter  1  can  be  extended  in 
a  natural  way  to  cases  where  A(t)  is  no  longer  a  matrix,  but  rather  a 
more  complex  linear  operator.  In  this  chapter  such  a  generalization  is 
made  in  the  following  direction:  x  is  assumed  to  be  a  scalar  function 
of  two  variables  t  and  s  as  well  as  of  the  parameter  y.  ,  while  A 
is  assumed  to  be  an  integral  operator,  integration  being  with  respect  to  s 
(Analogous  problems  in  a  more  abstract  form  are  considered  in  [28].) 

Thus,  we  consider  the  equation 

U  =  -[x(t,s,y)  -  J  K(s,o)x(t,o,y)dn  ] 

(1) 

+  yf(x,t,s,y,)  (0  <  t  <  T,  a  <  s  <  b)  , 
which  can  be  written  in  the  simple  form 

u||  =  Ax  +  yf  S  -  [X-Bx]  +  yf  ,  (1*) 

and  we  prescribe  the  initial  condition 
x(0,s,y)  =  x°(s)  . 


(2) 


Let  us  assume  that  X  =  1  is  an  eigenvalue  of  the  operator  E  . 
the  operator  A  has  the  eigenvalue  X  =  0  ,  and  so  the  reduced  equati  r. 
Ax  =  0  has  a  family  of  solutions  which  depends  on  one  or  more  arbitrary 
functions  of  t  .  The  same  questions  arise  here  as  in  the  previous 
chapters:  To  which  member  of  this  family  does  the  solution  of  the  prat  It 
(1)  ,  (2)  converge  as  y,  -  0  ,  that  is,  how  do  we  determine  the  function  r 
t  in  the  family  of  solutions  for  the  reduced  equation  which  provider  the 
limiting  solution  of  (1),  (2)?  What  does  the  asymptotic  expansion  of  the 
solution  of  this  problem  with  respect  to  y,  look  like? 

Suppose  that  the  following  conditions  are  satisfied: 

I.  K(s  ,a)  is  continuous  in  the  square  R  =  fa<s<b  ,  a<o<b} 
f(x,t,s,u)  is  continuous  with  respect  to  s  and  sufficiently  smooth 
with  respect  to  x,t  and  u  in  a  domain  D(x,t,s,y,)  =  D(x,t,s)  x  fO>uQ] 
where  D(x,t,s)  is  some  domain  in  the  space  of  the  variables  (x,t,s), 
and  x  (s)  is  continuous  for  a  <  s  <  b  . 

II.  The  kernel  K(s,o)  is  symmetric,  that  is,  K(s,o)  =  K(o,s) 
in  R  . 

III.  The  eigenvalues  of  the  operator  B  are  such  that 

=  . . .  =  X^  =  1  ,  X^  <  1  for  i  =  k  +  1,  . . . 

The  corresponding  set  of  eigenfunctions  {cp^(s)}  is  assumed  to  be  ortho¬ 
normal,  that  is, 


In  what  follows  we  will  denote  the  scalar  product  of  two  functions  u(sl 
and  v(s)  by 


(As  we  will  see  below,  this  is  an  equation  for  a  boundary  function.) 

We  make  the  following  assumptions: 

1°.  y°(s)  is  continuous  for  a  <  s  <  b  . 

2°.  The  operator  B  satisfies  Conditions  I- III. 

3°.  g(T,s)  is  continuous  for  t  >  0  and  a  <  s  <  b  ,  and  satisfies 
the  estimate  |g(T,s)|  <  cexp(-KT)  ,  where  c  and  k  are  positive  con¬ 


stants  . 


1 


0,  . 

We  now  pose  the  following  questions:  How  does  one  choose  y  (s)  in 
order  that  the  solution  of  the  problem  (3),  (U)  will  converge  to  zero  as 
T  -*  «  ?  Will  this  convergence  be  of  exponential  type,  that  is,  do  we 
have  an  estimate  of  the  same  type  as  that  for  g  in  3°  ?  The  answers 
to  these  questions  are  given  by 

Lemma  1.  Suppose  that  Conditions  1°  -  3°  are  satisfied,  and  that 
the  function  y°(s)  is  such  that 

CD 

<yn(s),  C0.(s)>  =  -  f  (g(f,s),  w.(s))dT  (i  =  l,...k).  (5) 

u  1  •  Q  1 

Then  the  solution  y(T,s)  of  the  problem  (3),  (1+)  exists  and  is  continuous 
for  T  >  0  and  a  <  s  <  b  ,  and  converges  to  zero  as  t  -»  ®  uniformly  for 
s  in  [a,b]  ;  moreover, 

|y(T,s)|  <  cexp(-nT)  .  (6) 

(We  note,  as  ir  previous  chapters,  that  the  constants  c  and  h  in 
various  estimates  of  the  type  (6)  sire  not  generally  the  same,  despite 
the  fact  that  they  are  denoted  by  the  same  letters.) 


Proof.  It  is  not  difficult  to  prove  the  existence  and  the  continuity 
of  the  solution  in  the  stated  domain  if  we  pass  from  (3),  (4)  to  the  inte¬ 
gral  equation 


b 

*  IJ  K(.  ,o)y(0,o)do  +  g(9,s)  ]d0 

'  a 


(7) 


7*nrx 


1 


13 


) 


and  therefore  it  follows  from  (11)  that 


00  ^  CD 

=  -  r  <g,CD.)iT+  f  <g,cp.)dT  =  -  r  <g,co  >dT  . 

“  q  -*■  q  1  1 

In  view  of  Condition  3^  we  have  that  52^(t)  »  consequently, 

satisfy  an  exponential  estimate,  namely 


62(t,s'. 


| s2i( T)  I  <  c  exp(-KT)  ,  |52(t,s)  |  <  c  exp(-nT) 
(t  >  0  ,  a  <  s  <  b)  . 


(12; 


From  (10)  we  obtain  now  for  5^  the  equation 


t&i  -  B81]  +  g1  , 


(13) 


where 


k  d8  k 

gl  *  g  -  1t'1~dT-<Pi  =  6  "  **  »  lej  <  cexp(-KT)  . 


Taking  the  scalar  product  of  (13)  with  2£^  ,  we  obtain 


^  <61,51>  =  -2[<&1,81)  -  <B1,B81>  ]  +  2(8^) 


By  the  Hilbert-Schmidt  Theorem  B61  =  E  X,  <  6_ ,  cp.  >  tp.  ,  where 

i=l  11  1 

the  summation  begins  with  i  =  k+  1  because  of  (9)«  Hence, 


=  E  X  <8.,cp.)2  .  (15) 

A  A  i=k+l  1  x  1 

Since  X^  <  1  for  i  *  k+  1,...  and  since  X  =  1  is  not  a  limit  point  of 
the  spectrum  of  B  ,  there  exists  a  positive  constant  X  <  1  such  that 
<  X  for  i  «  k+1,...  .  Further,  Bessel's  inequality 


£  >  <(B.,6  >  implies  that 

i=k+l  1  1  1  1 

<61,B81>  <  X  v 


<6^5^  -  <61,B81>  >  (1-X)  (b^bj  ,• 


where  X  <  1  . 


Let  us  rewrite  equation  (1U)  as 


dT^l*^?  =  -2D(t)<B1,B1>  +  2<&1>g1)  , 


where 


<81,81>  -  <81,B51> 


D<T) '  — ^ 


>  1  -  X  >  0  . 


Dividing  by  2(5^6^)' 


11  11  (b^)1' 


and  integrating,  we  obtain 


<51»&1>:1^2  =  <61(O»s),81(O,s))l/2exp(-fD(0)d0)  + 

JTexp(-  JTD(0)d0)<81,g1  X6i"6i>'VadS 


By  virtue  of  Cauchy's  inequality  and  Condition  3 


l<81*g1>l  5  <B1»S1>l/2<g1.61>l/2  <  <81,&1>l/2exp(-HT) 


<8i,Bi>^2  <  c  exp(-KT)  (t  >  0,  a  <  s  <  b) 


Thus 


!i 


14] 


By  substituting  (20)  into  (3),  (U)  and  equating  coefficients  of  like 
powers  of  p,  ,  separately  for  coefficients  depending  on  (t,s)  and  those 
depending  on  (t,s)  ,  we  obtain  equations  and  supplementary  conditions 
for  the  determination  of  x^(t,s)  and  TT\x(t,s)  .  By  substituting  (20) 
into  the  nonlinear  function  f(x,t,s,u)  we  can  express  f  as  f  +  TTf  , 
just  as  in  previous  chapters. 

We  obtain  first  the  equation 
b 

xQ(t,s)  =  J*  K(s,o)x0(t,o)do  • 

a 

Bu  virtue  of  Condition  III  its  general  solution  is 
__  k 

x  (t,s)  =  r  a.(t)q>  (s)  ,  (21) 

u  i=l  1  1 

where  the  a^(t)  are  as  yet  arbitrary  functions. 

The  equation  and  the  initial  condition  for  tTxq(t,s)  are 

Btr  x(t,s)  b 

- &r - =  ’  [V(T>s)  '  J  k(s,o)tr0x(T,o) do ]  , 

a 

q  _  n  k 

ff0x(°>s)  =  x  (s)  -  xo(0,s)  *  x  (s)  -  Z  a.(0)co  (s)  . 

i=0  1  1 

In  addition,  as  is  our  custom,  we  impose  the  restriction  that 

TT0x(t,s)  -  0  as  T  -  •  .  (22) 

Thus  we  have  for  tT^x(t,s)  a  problem  of  the  type  considered  in  Subsection  2 
of  5l  for  the  case  g(t,s)  =  0  .  By  virtue  of  Lemma  1,  TrQx(T,s)  satisfies 
the  condition  in  (22),  and  moreover,  the  inequality 


1 


| TTqX ( t ,  s )  |  <  c  exp(-KT)  (t  >  0,  a  <  s  <  b)  ,  (25 

provided  the  initial  value  TTqX(0,s)  satisfies  the  condition  in  (5),  that 

is, 

k 

<x  (s)  -  Z  a.(0)co.(s)  ,co.  (s))  =  0  (i  =  1, . . .  ,k)  .  ( 2U ) 

J=1  J  J  1 

Hence , 

0^(0)  =  <x°(s)  ^(s))  (i  *  1 . k)  .  (2-y 


As  was  the  case  in  previous  chapters,  we  now  determine  the  function: 
O^t)  completely  by  considering  the  equation  for  x^(t,s),  namely 


axQ(t,s) 

__ 


=  -  [x^tjS)  -  r  K(s,o)x1(t,o)do  1 


or 


+  f(xQ(t,s) ,t,s,0) 


x^( t , s)  -  J*  K(s,a)x1(t,o)dc  = 

a 

k  k  da.  (t) 

f(  Z  a.(t)©.(s)  ,t,s,0)  -  Z  — tt —  cp. ( s)  s  t ( t, s)  . 
i«=l  1  1  i=l  1 


(26) 


For  the  solvability  of  this  inhomogeneous  problem  it  is  necessary  and 
sufficient  that  the  right-hand  side  f(t,s)  be  orthogonal  to  each  cpi( s) 
(i  ■  1, ...,k)  .  This  orthogonality  condition  is  itself  represented  by  a 
system  of  differential  equations  for  a^t)  ,  namely 


k 

<f(  Z  a.(t)<j>.(s)  ,t,s,0) ,  <p.(s)>  (i 
J-l  J  J  1 


1, • . •  ,k) 


(27) 


IV.  Suppose  that  the  system  (27)  together  with  the  initial  conditions 
(25)  has  the  solutions  =  cr(t)  for  0  <  t  <  T  . 

Thus  the  function  XQ(t,s)  is  determined  completely  by  the  formula  (211 , 
and  the  construction  of  the  zeroth  term  in  the  expansion  is  finished. 

Let  us  now  introduce  in  the  space  of  the  variables  (x,t,s)  a  surface 
L  which  consists  of  the  two  parts: 

Lx  =  f(x,t,s):  x  =  x0(0,s)  +  TTqX ( t> s )  ( t  >  o)  ;  t  =  0;  a<s  <b}  , 

l>2  =  {(x»t,s):  x  =  x0(t,s)v  0  <  t  <  T;  a  <  s  <  b]  . 

It  is  natural  to  require  that  the  following  condition  holds: 

V.  Suppose  that  the  surface  L  belongs  to  the  domain  D(x,t,s) 
which  appears  in  Condition  I. 

The  general  solution  of  equation  (26)  can  be  written  as 
_  k 

0i(t)cpi(s)  +  x^t.s)  ,  (28) 

where  the  p^(t)  are  as  yet  arbitrary  functions  and  x^tjs)  is  a 
particular  solution  of  (26)  which,  for  example,  has  the  form 

-  X 

Mt,s)  *  ♦(tjs)  +  t  T-r~  =  <*(t,s)  ,cp.(s)^  . 

A  i=k+l  •L’Ai  1  1  1  1 

The  equation  and  the  supplementary  conditions  for  7T^x(t,s)  are 
dTT.X^s)  b 

- §5 -  =  -  [TT^t.s)  -  J  K(  S,o)rr1X(T,o)do  1  +  g(T,s)  , 


where 


g(T,s)  =  TT^f  =  f(xQ(0,s)  +  TTqX(t} s)  ,  0,  s ,  0)  -  f (xQ(0,s) ,0,s,0)  , 
^(0,5)  =  -  x^O.s)  , 

tT^x(t,s)  -*  0  as  T  -  •  .  (29) 

By  virtue  of  (23)  the  function  g(T,s)  satisfies  the  exponential  estimate 
I g( s)  |  <  cexp(-HT)  (t  >  0,  a  <  s  <  b)  .  Thus  the  function  tt^x(t,s)  is 
the  solution  of  a  problem  of  the  type  considered  in  Subsection  2  of  §1  . 

And  by  virtue  of  Lemma  1  it  satisfies  condition  (29)  and  an  exponential 
estimate,  provided  that 

_  "  )~ 

<  -X  (0,s)  ,©  (s)>  =  -f  <g(T,s)  ,ep.(s)  /  (IT  (i  =  l,...,k)  . 

j.  x  ■  0  x 

By  inserting  here  the  expression  (28)  for  x^OjS)  ,  we  obtain  the  values 
Pi(0) ,  namely 

O 

3,(0)  =  r  <g(T,s) >®.(s))dT  (  i  =  1, . . . ,k)  . 

’  0  1 

Thus,  ^x^s)  is  completely  determined,  and  we  have  found  the  initial 

values  3^(0)  .  The  functions  0^(t)  are  determined  completely  from  a 

solvability  condition  for  the  integral  equation  defining  x0(t,s),  in  a 

maimer  analogous  to  that  for  the  determination  of  the  functions  a.(t)  . 

i 

It  turns  out  that  the  p^(t)  satisfy  the  system  of  linear  differential 
equations 

)c 

=  t  bij(t)3j(t)  +  fx(t)  (i  =  l,...,k),  (30) 

vhere  b^(t)  «  (fx(xQ(t,s)  jtjSjOjcp^s)  and  f^(t)  is  a  known  function. 


The  determination  of  the  remaining  terms  in  the  expansion  proceeds 
according  to  considerations  analogous  to  those  used  for  the  determination 


of  x^(t,s)  and  tT^x(t,s)  .  At  the  i-th  step  the  expression  for  xi(t,s) 
contains  k  arbitrary  functions  (let  us  denote  them  by  y^(t),  i  =  l,...,k) 
Now  the  problem  for  h^Tjs)  is  analogous  to  the  one  for  tt^x(t,s)  , 
while  the  values  y^(0)  are  found  from  a  condition  like  (5)  •  Finally, 
we  obtain  a  system  of  linear  differential  equations  like  (30)  for  y^(t) 
from  a  solvability  condition  in  the  equation  for  x^+^(t,s)  .  It  follows 
also  that  each  ^-function  satisfies  an  exponential  estimate 

|tT^x(t,s)  I  <  c  exp(-«T)  (t  >  0,  a  <  s  <  b) 

2.  An  Estimate  of  the  Remainder  Term.  Let  us  first  make  more  pre¬ 
cise  the  requirement  involving  the  smoothness  of  f(x,t,s,p)  (cf.  I)  . 

It  is  possible  to  take  as  the  domain  D(x,t,s)  occurring  in  Condition  I 
an  arbitrary  8-tube  of  the  surface  L  (cf.  V)  .  We  then  require  that 
f(x,t,s,p)  have  continuous  partial  derivatives  up  to  order  n+  2  inclu¬ 
sive  with  respect  to  x,t  and  p  in  the  domain  D(x,t,s,p)  =  D(x,t,s,p) 
x[0,Pq]  .  We  have  determined  the  terms  of  the  series  (20)  to  order  n+  1 
inclusive,  and  let  us  denote  by  X^(t,s,p)  the  k-th  partial  sum  of 
(20),  that  is, 

k  1  - 

\(t»s,p)  =  E  p  (x.(t,s)  +  .x(t,s))  . 

K  i=0  1  1 

Let  us  now  introduce  the  norm  of  a  function  y(t,s,p)  by  ||y(t,s,p)||  * 
8«p{|y(t,s,p)  | :  0<t<T,  a<s<b1  . 


Theorem  U.  Under  Conditions  I-V  there  exist  positive  constants 


and  c  such  that  for  0  <  U  <  the  solution  x(t,s,u)  of  the  problem 
(1) ,  (2)  exists  in  the  domain  {0  <  t  <  T,  a  <  s  <  b]  ,  is  unique  and 
satisfies  the  inequality 

||x(t,s,n)  -  Xn(t,s,u)!l  <  c^1  . 


Proof.  (Suggested  by  A.  Kasanov.)  Set  P(t,s5u)  =  x(t,s,u)  -X  „(t,s 

rn*  x 

Then  by  substituting  x  =  X^^  +  into  (1),  (2)  we  obtain  for  f  the  pr 

lem 


u||  =  -[e-Be]  +  ufx(t,s,u)e  +  c(rft,s,u)  , 


(3D 


?(0,s,^)  =  0  , 


(3?) 


where 


fx(t,s,n) 


=  fx(xQ(t,s)  +  TTox(T,s),t,s,u)  and 


G(f,t,s,u)  =  -[X^  -  BXih1J  +  uf(Xrifl+  f,t,s,n) 


ufx(t,s,u)r  -  ji- 


dxnfl(t,s,li) 


dt 


As  in  the  previous  chapters,  it  is  not  difficult  to  show  that 
G(^,t,s,ii)  has  the  following  two  properties: 


2.  G(S,t,s,n)  is  a  contraction  operator  with  contraction  coefficient 

p 

of  order  G(u  )  for  ||5j|  =  &(u)  .  This  means  that  if  j!  £^(  t, s,u) 11  < c^u 
811(1  —  cili  »  then  there  are  constants  cQ  and  uQ  <  u  such 

that  for  0  <  u  < 

llG(e1,t,saM.)  -  G(e2,tfs,u)||  <  CqU2!!^  -  Fjl  .  (33) 

See  our  remark  in  Subsection  3,  §1  of  Chapter  1  regarding  the  constant  uQ  • 
Let  us  now  write  (31)  as 

U?  =  G(e;,t,s,n)  ,  (3U) 

where  U  is  the  operator  defined  by  US  =  |j,-~  +  [p-BP]  -  uf^F  . 

We  consider  first  the  auxiliary  linear  problem 

Uy  =  a(t,s,ii),  y(0,s,u)  =  0  ,  (35) 

where  a(t,s,u)  is  a  given  continuous  function. 

Lemma  2.  For  each  sufficiently  small  value  of  g.(0  <  u  <  uQ)  the 
solution  y(t,s,u)  of  problem  (35)  exists  in  the  domain  {0  <  t  <  T  , 
a  <  8  5  »  ls  unique,  continuous  and  satisfies  the  inequality 

lly(t,s,u)||  <  My,'1)|a(t,s,u)||  ,  (36) 

where  the  constant  M  >  0  is  independent  of  w  . 


Remark.  The  constant  in  this  estimate,  because  of  its  importance  for 


later  discussion,  is  denoted  by  M  rather  than  by  c  . 

The  existence,  uniqueness  and  continuity  of  y(t,s,u)  can  be 
established  without  difficulty  by  passing  from  (35)  bo  the  integral  equation 

y(t,s,u)  =  r  expor1r  (-1+  Uf  )d©)  de  (37^ 

0  e  x  u  * 

and  applying  the  method  of  successive  approximations. 

Let  us  prove  that  the  estimate  (36)  is  valid.  By  taking  the  scalar 
product  of  equation  (35)  with  y  we  obtain  the  equation 

•§u~[£<y,y>  =  -<y,y>  +  P<y,y>  +  uQ(y,y)  +  R(y,y>1//2  ,  (38) 

where 


We  know  already  from  an  analogous  situation  that 


<y,By>  =  <y,  t  \  <y,©.)co.>  =  t  X.<y,q>.>2 
i=l  1  11  i=l  1  1 

CD 

<  t  <y,cp.)2  <  <y,y>  , 

i=l  1 

so  P  <  1  .  Moreover,  it  is  clear  that 

||qj|  <  ||fxll  and  ||R||  <  lloll(b-a)1/2  . 
Dividing  (38)  by  (y,y)^2  we  obtain 


(39) 
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u-^  (yjy'*1/2  =  (-i+  p+  uQXyjy^1^2  +  R  , 
and  so 

(yjy^1/2  =  J  exp()i  1r  (-1+  P+uQ)d0)  •  —  d© 

0  ’  0  u 

By  virtue  of  the  properties  just  noted  P  and  Q  are  exponentially  bounded, 
and  consequently,  by  taking  account  of  the  inequality  in  (39)  for  R  ,  we 
have  the  estimate 

<y,y>1//2  <  cy,  ^id!  .  (40) 

From  the  obvious  inequality 

-1  ^  -1 
exp(y,  f  (-1+uf  )d0)  <  exp(-u  n(t-0))  , 

J0  x 

(37)  implies  now  that 

|(y(t,s,u)!l  <  c(|IBy)l  +  Hall)  .  (4l) 

Since 

| By |  =  |  f  K(s,o)y(t,a,u)  do|  <  ^K,K>1//2<y,y>1//2  <  ^lldl  , 

’  a  ^ 

it  follows  from  (40)  that  ||Byjj  <  ^11  all  ,  and  so, 

Ily1l<3la|l  • 

If  we  replace  the  constant  c  by  M  then  the  lemma  is  proved. 

We  turn  now  to  equation  (34)  and  the  initial  condition  (32).  We  will 
prove  the  existence  of  a  unique  solution  satisfying  the  estimate 


II  til  =  Mu*1) 


for  all  sufficiently  small  p  .  This  will  establish  the  theorem. 

'  In  order  to  apply  the  method  of  successive  approximations  we  set 
»  0  and  define  inductively  by  =  G(fk,t,s,u)  .  By  vir¬ 

tue  of  Property  1  for  G  ||G( t , s ,ut) II  <  c \i*2  ,  while  from  Lemma  2 

c 

ii.  H  -  M  n+2  ..  m-1  1  rnl  , 

II  Ijjl  <  ~  ^U  =  Mcp,  •  (^3 

I’M*  1 

We  continue  now  by  induction  on  indices  knowing  that  ||5j!  <  c^u  and 
||  l[^||  <  c^p  .  If  we  assume  that 

ll?iH  <  c^1  (i  =  0,1,... k)  ,  (kb 

then  we  must  prove  that  (U4)  is  also  valid  for  i  =  k  +  1  .  Clearly 

^Vl’V  "  G<?k,t,s,p)  '  » 

m.  1 

and  since  ||  f^H  <  c^  <  c^  and  Hf^ll  <  c^p  for  0  <  p  <  <  1  , 

we  can  apply  the  inequality  (33)  to  the  difference  in  the  right-hand  side. 
In  addition,  by  using  Lemma  2,  we  obtain 

-] I  C0U^  h. '  -k-l'1  =  Mc0^^k'  Vl11  • 

We  now  choose  p^  so  small  that  the  inequality  Mc^p^  <  -|  is  satisfied. 
Then,  taking  account  of  (U3) ,  it  follows  that 

i  -  hJ  S  •  •  • 

(^5) 

Sy"V*ol'  -7"¥5-^ri  'Z1  • 


Hence 


lyfvi-y  *  “5k-  w  ♦  - + "  fi" 


,11  1.  TM-l  . 

ici’ 


iw- 1 


Thus,  (44)  and  therefore,  (U5)  are  valid  for  any  number  k  .  By 
virtue  of  <1+5)  the  series  of  terms  E|k+1  -  converges  uniformly  with 
respect  to  (t,s),  that  is,  the  sequence  {^}  converges  uniformly. 

This,  in  turn,  implies  the  existence  of  a  solution  of  the  problem  (34), 
(32).  To  see  this,  simply  write  equation  (34)  as  an  integral  equation 
(from  (37)  it  suffices  to  write  y  for  E  and  a  for  G  ) .  The  unique¬ 
ness  of  the  solution  of  this  integral  equation  can  be  proved  in  the  usual 
way,  if  we  note  that  the  operator  BE  +  G(E,t,s,p.)  is  Lipschitzian. 
Finally,  the  estimate  (42)  follows  from  the  fact  that  each  Ek  satisfies 
the  inequality  (44) .  This  concludes  the  proof  of  Theorem  4. 


3.  Concluding  Remarks.  1.  Our  results  can  be  extended,  under  cer¬ 
tain  additional  assumptions,  to  the  more  general  equation 

ir^t^5^)  r  -a(t,a)  [x(t,s,n)  -J  K(t,s,o)x(t,o,u)do] 

b  "  <*6> 
+  nf(x,t,s,J  H(x,t,o,p)dc,u)  • 

’  a 

The  most  interesting  feature  of  this  equation  is  the  appearance  of  the 
factor  a(t,s)  ,  whose  dependence  on  s  destroys  the  symmetry  of  the 
operator  A  . 

Let  us  first  investigate  how  to  modify  the  construction  in  the  special 
case  when  equation  (46)  is  of  the  form 


13- 

=  -a(s)  [x(t,s,n)  -  J*  K(s,o)x(t,a,ti)da] 

a 

+  ^f(x,t,s,M>)  . 

(This  equation  has  been  studied  by  A.  Kasanov.)  For  xQ(t,s)  and  TTox(t,s) 
we  have  that 

_  k 

x  (t,s)  =  £  a.(t)co.  (s)  , 

u  i=l  1  x 

3Tf  x(T,s)  b 

- ^ -  =  -a(s)  [tT0x(t,s)  -  f  K(s,o)TT0x(T,o)dCT]  , 

a 

TrQx(0,s)  =  x°(s)  -  x0(0,s)  , 

and 

tTqX (  t  ,  s )  -  0  as  T  -  •  •  (*+7) 

Condition  (47)  reduces  to  the  condition  that  the  initial  value  Tt^xf  0,s)  be 
bounded,  that  is, 

b  TTnx(0,s)tp  (s) 

J  °a(3)  1  '  a.  -  0  (i 

a 

(cf.  (24));  whence,  0^.(0)  can  be  defined  by  the  formula  (cf.  (25)) 
b  x°(s)a  (s) 

ai^0^  =  I  a(s)  ds  (i  =  *■»  ■••>*)  * 

a 

We  must  of  course  choose  the  eigenfunctions  <p^(s)  (i  =  l,...,k)  to  be 
orthogonal  with  weight  l/a(s)  . 

The  equations  for  the  cr^(t)  are  obtained,  as  before,  from  a  solv¬ 
ability  condition  in  the  equation  for  x^(t,s) ,  that  is, 
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« 


0 


Sxn(t,s)  b 

— — -  =  -a(s) [x^(t,s)  -  J  K(s,o)x1(t,o)do] 

& 

+  f(xQ(t,s),t,s,0)  , 


and  they  are 


dor,  b 

«Kr 

a 


k 

f(  t  a .(t)cp.(s)  ,t,s,0)cp,(s) 

3-1  J  J _ _ _ 

a(s) 


ds 


2.  Equations  and  systems  of  the  form  p  —  =  A(t)x  +  pf  ,  where  A  is 

CL  v 

a  certain  integral  operator,  arise  in  kinetic  theory  (for  example,  Boltzmann's 
equation;  see  [6]) .  It  is  true  that  the  structure  of  the  operator  A  is 
more  complicated  in  such  problems  than  in  the  cases  considered  here;  moreover, 
the  problems  are  nonlinear.  However,  it  is  important  to  note  the  following: 
solutions  of  the  degenerate  equation  A(t)x  =  0  contain  an  indeterminacy. 

In  this  regard  we  note  that  seme  of  the  approximation  methods  in  kinetic 
theory  (cf.  [6])  lead  to  nonlinear  systems  of  the  type  considered  in  Chapter  2. 


* 
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